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SECTION  1 


SUMMARY 

A combined  experimental  and  analytic  program  was  conducted  to  obtain  required  data  and  develop 
improved  procedures  for  predicting  the  ablation  performance  of  graphitic  materials  in  rocket  nozzles. 
Particular  emphasis  was  placed  on  conditions  representative  of  the  throat  region  of  an  MX  nozzle  with 
graphitic  materials  which  are  representative  of  the  current  state  of  materials  development.  Signifi- 
cant results  of  this  program  are: 

• The  acquisition  of  a large  body  of  full  characterization  ablation  data  using  an  arc  plasma 
generator  which  simulates  rocket  nozzle  environments 

• The  harmonizing  of  this  data  with  available  roc'-et  motor  data  to  obtain  kinetic  correla- 
tion functions  which  may  be  used  to  predicc  ablation  response.  Correlation  functions  were 
obtained  for 

G-90  and  ATJ  bulk  graphite 

15%  silicon  carbide  modified  pyrolytic  graphite 

c plane  pyrolytic  graphite 

Pyrocarb  901  (p  ~ 1.83  gm/cc} 

• The  development  of  the  GASKET2  code.  This  code  is  a modification  of  the  GASKET  code  to 
include  the  new  correlation  functions,  new  JANNAF  species  data,  improved  generalized  in- 
put requirements,  and  kinetically-controlled  carbon  sublimation  models. 

• An  evaluation  of  rocket  nozzle  thermal  performance  procedures  currently  used  by  rocket 
motor  designers  to  predict  the  thermal  performance  of  a solid  propellant  rocket  nozzle. 

A significant  conclusion  of  this  evaluation  is  that  errors  in  the  mass  transfer  coeffi- 
cient have  a direct  and  significant  influence  on  the  predicted  ablation  rate.  With  the 
exception  of  the  kinetic  model,  all  other  variables  were  shown  to  have  a small  effect  on 
the  mass  loss  rate  in  the  throat  region  of  a typical  MX  nozzle.  Thus,  future  efforts 
should  concentrate  on  accurate  predictions  of  the  mass  transfer  coefficient  and  accurate 
kinetic  models. 
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• An  analytic  evaluation  of  the  relative  performance  of  different  graphitic  materials  in 
tiie  MX  throat  and  nose  cap  environments 

• Pretest  predictions  of  the  ablation  performance  of  several  AFRPL  advanced  technology  rocket 
nozzles 

• The  acquisition  of  limited  characterization  ablation  data  for  comparing  the  relative 
ablation  performance  of  materials  in  given  generic  classes.  These  classes  and  the  limited 
characterization  materials  were 

- Bulk  graphites  - P03 

- ATJS 

- Pyrolytic  graphites  - Hitco  a-b  plane 

- Pfizer  c plane 
5%  SiC/PG 

- 23%  SiC/PG 

- 65%  HfC/PG 

- Carbon/carbons  - Pyrocarb  903  (p  ~ 1.83  gm/cc) 

- Pyrocarb  903  HD  (p  ~ 1.90  gm/cc) 

- HRX  5125  (p  ~ 1.50  gm/cc) 

- HRX  5875  (p  - 1.80  gm/cc) 

- MDAC  3D(p  ~ 2.00  gm/cc) 
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SECTION  2 


INTRODUCTION 

The  rate  at  which  graphite  nozzle  components  are  consumed  by  hot  propellant  gases  depends  upon 
the  composition  of  these  gases.  In  addition,  however,  it  is  known  that  different  generic  classes 
(e.j.,  oulk  graphites  or  pyrolytic  graphites)  respond  differently  to  the  same  environments.  Although 
higher  density  materials  often  have  lower  ablation  rates,  the  contrasting  ablation  performance  of 
a-b  oriented  pyrolytic  graphite  compared  with  c oriented  pyrolytic  graphite  clearly  demonstrates  that 
density  is  not  the  only  important  factor.  Past  experience  (References  1-5)  has  shown  that  surface 
reaction  kinetics  have  an  important,  if  not  dominant,  role  in  controlling  carcon  consumption  rates. 

With  current  technology  it  is  not  possible  to  write  elementary  chemical  reactions  which  will 
define  the  carbon  consumption  process.  Yet,  some  basic  analytic  expressions  or  procedures  are  re- 
quired for  design  and  performance  analyses  of  rocket  motors.  The  Air  Force  Rocket  Propulsion  Labora- 
tory recognized  this  need  and  initiated  a program  (Reference  6)  to  develop  a semiemperical  procedure 
to  predict  ablation  response.  At  that  time  a-b  oriented  pyrolytic  graphite  was  considered  to  be  a 
viable  rocket  nozzle  material.  Thus  Reference  6 concentrated  on  the  development  of  a general  analytic 
model  and  its  validation  using  a-b  oriented  pyrolytic  graphite  data. 

The  analytic  model  was,  by  necessity,  based  on  an  engineering  approximation  of  the  overall 
ablation  reaction.  The  ablation  data  used  for  determining  the  correlation  coefficient  in  this  model 
were  obtained  by  exposing  a large  number  of  models  to  simulated  hot  propellant  environments.  The 
resultant  analytic  model  and  the  correlations  for  pyrolytic  graphite  were  assembled  into  the  GASKET 
code  (Reference  7).  The  GASKET  code  was  then  integrated  into  an  overall  rocket  nozzle  ablation  per- 
formance prediction  procedure. 

Since  the  completion  of  the  above  work,  thennostructural  problems  associated  with  a-b  oriented 
pyrolytic  graphite  have  cast  doubt  on  their  use  as  rocket  nozzle  liner  material.  Other  materials 
considered  and  applied  to  rocket  nozzles  include  bulk  graphites,  carbon/carbons,  c oriented  pyrolytic 
graphites  and  carbide  codeposited  pyrolytic  graphites.  A need  therefore  existed  to  generate  corre- 
lation coefficients  for  these  other  materials.  Because  of  the  many  varieties  of  materials  in  each 
of  these  generic  classes,  it  is  not  economically  practical  to  kinetically  characterize  all  possible 
graphitic  materials. 
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To  satisfy  the  requirement  for  predicting  ablation  response  of  other  graphitic  materials,  the 
program  described  in  this  report  was  conducted.  The  program  included  full  characterization  of  select 
materials  from  each  generic  class  and  limited  characterization  of  a larger  number  of  materials.  The 
behavior  of  the  limited  characterization  materials  was  deduced  by  comparing  their  ablation  response 
to  that  of  fully  characterized  materials  in  the  same  generic  class.  Appropriate  data  was  obtained 
from  controlled  ablation  tests  using  an  arc  plasma  generator  and  analyses  of  available  rocket  motor 
data.  The  data  analysis  procedure  was  very  similar  to  that  described  in  Reference  8. 

Data  for  full  characterizations  were  obtained  for: 

• ATJ  bulk  graphite 

• G-90  bulk  graphite 

• 15%  SIC  codeposited  pyrolytic  graphite 

• Pyrocarb  901  carbon/carbon  (p  = 1.83  lbm/ft3) 

Data  for  limited  characterizations  were  obtained  for: 

t Pyrocarb  903  carbon/carbon  (p  = 1.82) 

• High  Density  Pyrocarb  903  carbon/carbon  {p  1 1.90) 

• HRX  5125  carbon/carbon  (p  = 1.50) 

• HRX  5875  carbon/carbon  (p  = 1.86) 

• MDAC  3-D  carbon/carbon  (p  = 2.02) 

• ATJ-S  bulk  graphite 

• 5%  SiC  codeposited  pyrolytic  graphite 

• 23%  SIC  codepositied  pyrolytic  graphite 

• 65%  HfC  codeposited  pyrolytic  graphite 

• c oriented  pyrolytic  graphite  (Pfizer) 

• a-b  oriented  pyrolytic  graphite  (riitco) 

Full  characterization  data  was  analyzed  to  obtain  appropriate  correlation  coefficients  which  were 
compatible  with  the  GASKET  code.  The  GASKET  code  was  modified  to  include  this  new  data  base  and  des- 
ignated as  GASKET2  (Reference  9).  The  validity  or  accuracy  of  the  GASKET2  code  was  then  assessed  by 
tvmparing  predicted  and  post-test  measured  responses  of  rocket  motors. 


Each  of  the  test  materials  were  examined  for  pre-  and  post-test  microstructural  and  chemical 
characteristics.  Characterization  data  included  photomicrographs,  scanning  electron  microscopy, 
porosity,  chemadsorption,  thermal  expansion  coefficients,  and  lattice  parameters.  These  data  were 
examined  and  where  possible,  the  ablation  performance  was  related  to  the  observed  characteristics. 

The  GASKET2  code  was  also  used  in  a series  of  performance  calculations  (References  10-14)  to 
predict  the  response  of  various  graphitic  materials  in  propellant  environments  similar  to  those  ex- 
pected for  an  MX  rocket  motor.  These  analyses  were 

1.  Third  stage  Hercules  MX  nozzle  analysis  (Reference  10) 

2.  Rocketdyne  Condor  test  nozzle  analysis  (Reference  11) 

3.  Third  stage  Hercules  MX  nozzle  analysis  using  Pyrocarb  901  kinetics  (Reference  12) 

4.  Material /propellant  sensitivity  study  for  the  throat  location  of  the  standard  7-inch 
test  nozzle  (Reference  13) 

5.  Material  sensitivity  study  for  the  nose  cap  of  the  standard  7 inch  test  nozzle 
(Reference  14). 

In  addition  to  the  above  thermal  ablation  studies,  a thermostructural  analysis  was  performed 
on  the  BATES  motor  to  identify  probable  causes  for  nozzle  failures.  The  results  of  this  study  are 
presented  in  Reference  15. 
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SECTION  3 


PROCEDURE  FOR  MATERIAL  CHARACTERIZATION 

The  surface  recession  of  carbon  rocket  nozzles  due  to  chemical  erosion  by  propellant  gases 
is  a complex  phenomenon.  It  includes  events  such  as:  diffusion  of  the  reactive  species  to  the 
carbon  surface,  adsorption  and  desorption  of  the  reactive  species  and  reaction  products  at  the  sur- 
face, heterogeneous  reactions,  and  diffusion  of  the  reaction  products  into  the  bulk  stream.  With 
the  present  state  of  the  art,  an  exact  analytical  model  describing  these  events  is  not  attainable. 
Hence,  an  engineering  approximation  representing  the  overall  observed  phenomenon  is  used  as  an 
alternative. 

Under  U.S.  Air  Force  Contract  F0461 1 -69-C-0081  (Reference  6),  a combined  analytical  and 
empirical  procedure  was  developed  to  correlate  the  ablation  rate  of  pyrolytic  graphite.  This  pro- 
cedure included  ablation  tests  under  simulated  environments  in  the  arc  plasma  generator,  determina- 
tion of  the  correlation  function,  and  identification  of  the  kinetically  controlled  reactions.  This 
procedure  was  applied  to  the  consumption  rate  of  a-b  plane  oriented  pyrolytic  graphite.  The  result- 
ing correlation  has  been  applied  extensively  in  designing  rocket  nozzles,  and  satisfactoiy  predic- 
tions have  been  obtained. 

Under  the  current  U.S.  Air  Force  Contract  F0461 1 - 74-C-0023 , carbon  materials  were  ihararter- 
ized  using  two  procedures.  The  first  procedure  used  a full  characterization  procedure  similar  to 
Reference  6.  Five  commonly  used  carbon  materials  were  selected  from  generic  types  of  pyrolytic 
graphite,  modified  pyrolytic  graphite,  bulk  graphite,  and  carbon/carbon.  These  underwent  the  full 
procedure  to  determine  correlations  of  their  ablation  rates.  The  second  procedure,  which  only  de- 
termined the  relative  ablative  performances,  was  called  limited  characterization.  Considerably 
fewer  APG  data  were  taken  for  this  type  of  characterization  since  the  determination  of  the  empirical 
correlation  was  not  required.  The  carbon  materials  which  were  characterized  are  shown  in  Table  1. 

This  section  will  describe  the  experimental  apparatus  of  the  Aerotherm  1-megawatt  arc  plasma 
generator,  its  application  to  the  material  characterization  procedure,  the  selection  of  the  APG 
test  gases  and  test  conditions,  and  the  correlation  procedure.  Additional  details  are  available  in 
Reference  16. 
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TABLE  1.  CARBON  MATERIALS  CHARACTERIZED 


Manufacturer 

Material 

Test3 

Number 
- of 

Source 

F 

S 

Models 



ARC 

15%  SiC/PG 

X 

33 

HITCO 

Pyrocarb  901 

X 

28 

UC 

ATJ 

X 

28 

Carborandum 

G-90 

X 

28 

ARC 

5%  SiC/PG 

X 

3 

ARC 

23%  SiC/PG 

X 

10 

Raytheon 

HfC/PG 

X 

10 

Hitco 

A-B  PG 

X 

10 

Super  Temp 

PG  Plate 

X 

10 

Pfizer 

PG  Plate 

X 

10 

Hitco 

Pyrocarb  903 

X 

10 

Haveg 

HRX-5125 

X 

10 

Haveg 

HR X- 5 875 

X 

5 

McDac 

3-D  C/C 

X 

6 

UC 

ATJ-S 

X 

Pure  Carbon 

P03 

X 

Hitco 

High  Density 
Pyrocarb  903 

— 

■ 

5 

aF  - full  characterization 
S - limited  characterization 
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3.1  EXPERIMENTAL  APPARATUS 


The  experimental  apparatus  consisted  of  the  arc  plasma  generator  used  to  produce  the  high 
temperature  reactive  environments,  the  test  nozzles  which  were  exposed  to  these  environments,  the 
fume  collection,  cooling  and  scrubbing  system  used  to  remove  the  test  gases  from  the  facility,  and 
the  instrumentation  used  to  characterize  the  test  conditions  and  model  response.  The  arc  plasma 
generator  and  support  equipment  are  discussed  in  Section  3.1.1.  The  test  nozzles  are  described  in 
Section  3.1.2.  The  fume  collection  system  is  described  in  Section  3.1.3  and  the  instrumentation 
is  presented  in  Section  3.1.4. 

3.1.1  Arc  Plasma  Generator 

The  Aerotherm  1-megawatt  constricted  arc  plasma  generator  (AP6)  is  shown  schematically  in 
Figure  1 and  physically  in  Figure  ?.  The  AP6  is  a constant  mass  flowrate  device  with  a flowrate 
controlled  by  throttling  at  the  gas  injection  ports.  The  APG  uses  a segmented  constrictor  arc  with 
a tungsten  cathode  and  a water-cooled  copper  anode  to  transfer  energy  to  the  primary  test  gas. 

This  test  gas  i;  injected  tangentially  between  the  cathode  and  the  first  constrictor  segment  to  pro- 
vide a stable,  high  voltage  operation.  Additional  gases  to  simulate  propellant  gases  are  injected 
downstream  of  the  anode  and  mixed  with  the  primary  arc-heated  gas  in  a plenum  chamber.  Thermo- 
chemical equilibrium  is  achieved  in  this  plenum  and  the  resulting  simulation  gases  are  expanded 
through  a choked  converging-diverging  nozzle.  The  test  section  is  the  throat  region  of  this  nozzle. 

The  arc  unit  is  water-cooled  with  ambient  temperature,  high  pressure  deionized  water.  The 
APG  input  power  is  supplied  by  a 600  kW  continuous  rated,  saturable  core  reactor,  dc  rectifier 
power  supply.  A maximum  overload  power  level  of  1.2  MW  is  achievable  for  5 minutes.  The  power 
supply  has  1000,  2000,  or  4000  volts  open  circuit  voltage  modes  to  match  APG  operating  character- 
istics for  various  test  gases,  flowrates,  and  pressures.  Arc  starting  is  accomplished  by  imposing 
power  supply  open  circuit  voltages  across  the  APG  electrodes  while  an  argon  flow  is  maintained. 

Then  a momentary  RF  discharge  in  the  APG  column  provides  an  initial  ionization  path  for  the  arc. 

Once  the  arc  is  started,  test  gases  are  immediately  introduced  as  necessary  to  provide  the  required 
test  gas  composition. 

The  arc  unit  exhibits  very  low  contamination  levels.  Based  on  the  results  of  Reference  17, 
total  gas  stream  contamination  should  not  exceed  200  parts  per  million  (0.02  percent).  The  major 
sources  of  this  contamination  are  the  tungsten  cathode  and  copper  anode.  A third  potential  source 
of  contamination  is  the  boron  nitride  insulators  of  the  constrictor  section;  however,  their  contri- 
Dution  to  the  above  figure  is  felt  to  be  very  small. 
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Copper 

Tungsten  anode 


injection:  carbon 

monoxide  or  hydrogen 


Figure  1.  Aerotherm  1-MW  constrictor  arc  heater. 
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Figure  2.  Aerotherm  1-MW  constrictor  arc  heater. 


3.1.2  Test  Nozzle  Configurations 


The  nominal  test  configuration  was  an  axisymmetric  nozzle  as  shown  in  Figure  3.  The  test 
section  inserts  formed  the  throat  region  of  the  nozzle.  The  PG  washer  immediately  upstream  of  the 
test  section  insert  insured  a smooth  transition  into  the  insert  and  held  the  boundary  layer  trip. 
This  trip,  a thin  Grafoil  disk,  was  employed  to  promote  turbulent  flow,  and  therefore  high  transfer 
coefficients,  in  the  throat.  The  test  section  insert  was  retained  by  a crushable  high  temperature 
insulator  and  could  expand  thermally  without  suffering  excessive  compressive  stresses.  The  test 
section  insert  configuration  is  shown  in  Figure  4.  This  s the  nominal -dimension  configuration; 
the  details  of  the  actual  test  insert  depend  on  the  particular  requirements  of  the  test  material, 
e.g.,  backwall  insulation  in  the  throat  region. 

An  appropriate  ablation  sample  or  a water-cooled,  steady  state  calorimeter  was  placed  in  the 
test  section.  The  calorimeter  and  test  sample  both  have  the  nominal  interior  dimensions  shown  in 
Figure  4,  so  that  the  test  conditions  during  an  ablation  test  could  be  inferred  from  a correspond- 
ing calorimeter  test.  The  calorimeter  installation  is  shown  in  Figure  5 and  a view  of  the  assembly 
is  shown  in  Figure  6. 

3.1.3  Fume  Collection  System 

The  APG  for  these  tests  was  run  on  the  atmospheric  test  stand  with  the  test  gases  exiting 
directly  into  the  test  bay.  A fume  collection  system  was  employed  to  collect,  cool,  clean,  and  ex- 
haust the  gases  outside  the  test  area.  The  system  is  shown  schematically  in  Figure  7. 

The  first  component  of  the  system  is  the  heat  exchanger  section.  The  high  temperature  of 
the  test  gases  as  they  left  the  APG  required  a “cooldown11  to  less  than  250°F  before  they  entered 
the  remainder  of  the  system.  This  section  is  constructed  of  a high  temperature  alloy,  Hastelloy 
Alloy  C-276,  and  provided  a set  of  spray  nozzles  which  "quench"  or  cool  the  gases  with  a water 
spray.  Also  included  in  this  section  are  two  view  poi ts  to  allow  pyrometer  viewing  of  the  test 
section. 

The  gases  were  then  ducted  to  the  fume  scrubber  mounted  outside  the  test  bay.  This  scrubber 
is  of  the  packed  tower  type  and  is  designed  to  remove  all  toxic  fumes  (HCt,  HF)  from  the  gas  stream 
before  they  are  exhausted  to  the  atmosphere.  The  scrubbing  fluid  was  water  used  in  the  once-through 
mede. 

The  final  component  of  the  system  is  the  exhaust  fan,  mounted  on  the  roof  of  the  test  bay. 
This  provides  the  positive  draft  required  to  draw  the  gases  through  the  heat  exchanger  section  and 
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gure  6.  Axisymmetric  test  section.  Calorimeter  installed. 
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Figure  7.  Fume  collection  system. 


the  scrubber.  The  fan  has  been  sized  to  provide  a slightly  negative  pressure  in  the  system  when  it 
is  used  in  the  blanked-off  mode.  This  is  necessary  when  hazardous  or  toxic  test  gases  are  used  as 
it  prevents  the  release  of  such  gases  into  the  test  bay  and  insures  personal  safety. 

Due  to  the  corrosive  nature  of  certain  of  the  test  gases  (HC2.  or  HF),  the  fan,  scrubber  and 
all  ducting  exclusive  of  the  heat  exchanger  section  are  constructed  of  Rigedon  4837-AT-HF.  This  is 
a fire-retardant,  fiberglass-reinforced  polyester  plastic  resistant  to  corrosive  attack  by  both 
acids  and  alkali  and,  in  addition,  is  provided  with  a special  Dynel  veil  for  protection  against 
fluoride  attack. 

3.1.4  Instrumentation 

The  measurements  to  characterize  the  test  conditions  and  material  response  were: 

• Test  Condition 

- Gas  Total  Enthalpy,  hQ 

- Chamber  Pressure,  PQ 

- Cold  Wall  Heat  Flux,  q 

c « w • 

- Reactive  Gas  Composition, 

• Material  Response 

- Surface  Temperature  History,  Tw  and  6exper1mental 

- Surface  Recession,  S 

- Qualitative  Surface  Condition 

Tne  gas  total  enthalpy  was  defined  by  an  energy  balance  on  the  arc  heater  including  the 
plenum  chamber,  i.e. , 


h _ h = Ah  = Power  In-Cooling  Water  Losses  m\ 

no  " amb  " anarc  Total  Gas  Flowrate 

0.948  x 10'3  El  - m AT  C 

L C r_ 

♦ 

m„,„ 

gas 

where  hamp  is  the  enthalpy  of  the  test  gases  at  room  temperature.  Voltage  E and  current  I were  re- 
corded continuously  on  a digital  data  recording  system;  measurements  from  panel  meters  were  alsu 
taken  as  a check.  The  cooling  water  flowrate,  m , was  measured  continuously  during  each  test  with 
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a sharp-edged  orifice  and  differential  pressure  transducer  and  its  temperature  rise,  ATC<  was  mea- 
sured continuously  with  a differential  thermopile.  The  total  gas  flow  rate,  mgas,  was  the  sum  of 
all  gas  flowrates  delivered  to  the  APG.  All  gas  flowrates  were  measured  with  ASME  sharp-edged 
orifices  and  differential  pressure  gauges,  except  hydrogen-chloride  which  was  measured  with  a 
rotameter  with  a magnetic  float  follower. 

The  chamber  pressure  was  measured  continuously  with  strain  gauge  pressure  transducers.  The 
pressure  taps  are  located  at  the  downstream  end  of  the  plenum-mixing  chamber  (Figure  3).  The  cham- 
ber temperature  was  determined  from  the  calibrated  net  enthalpy  addition  due  to  arc  heating,  the 
measured  chamber  pressure,  the  test  gas  composition,  and  an  ACE  computer  code  computation  of  chamber 
conditions. 


Cold  wall  heat  flux  was  measured  at  the  throat  of  the  water-cooled  copper  calibration  nozzle 
with  a steady  state,  waler-cooled  calorimeter  section.  The  coolant  water  temperature  rise  ATc  was 
measured  with  a single-pair,  copper-constantan  differential  thermopile,  the  output  of  which  was  re- 
corded continuously.  The  calorimeter  water  flow,  mc,  was  measured  with  a standard  glass  tube 
rotameter  and  the  heat  flux  then  calculated  from  the  equation: 


if 


',cw 


(2) 


where  Ac  is  the  calorimeter  heated  area. 

Surface  temperature  history  was  measured  with  a Thermodot  TD-9CH  optical  pyrometer  which  is 
calibrated  with  a high  temperature  source.  For  each  nozzle  ablation  test,  this  pyrometer,  which 
has  a sensing  wavelength  of  0.8  microns,  viewed  the  nozzle  throat  at  an  angle  of  approximately  40° 
from  the  APG  centerline.  Output  data  was  recorded  both  visually  from  the  instrument  meter  and  in 
digital  form  from  the  data  acquisition  system.  In  some  tests,  a second  pyrometer  was  used  as  a 
check  on  the  primary  unit.  This  secondary  unit  was  a Thermodot  T0-9FH  optical  pyrometer  similar  to 
the  primary  instrument  except  calibrated  in  degrees  Fahrenhe,,.. 

The  test  sample  surface  recession  was  obtained  from  pre-  and  post-test  measurements  of  the 
throat  diameter.  Measurements  were  made  at  three  axial  stations  in  the  throat  region,  namely,  the 
entrance,  center,  and  exit.  In  addition,  at  each  station,  the  diameters  were  determined  at  two 
angular  positions  90°  apart.  The  measurement  accuracy  is  approximately  ±0.0005  inch. 
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3.2 


TEST  GASES  AND  TEST  CONDITIONS 


3.2.1  Test  Gas  Selection  Criteria 

The  selection  of  gases  for  APG  testing  is  very  important  since,  ideally,  one  would  like  to 
minimize  the  extent  to  which  experimental  results  must  be  extrapolated  in  order  to  predict  actual 
conditions.  Three  questions  must  therefore  be  addressed  in  selecting  appropriate  gases. 

1.  What  rocket  motor  environments  are  anticipated? 

2.  What  are  the  important  surface  reactions? 

3.  What  are  the  operating  limitations  of  the  APG? 

Test  gases  must  be  defined  for  two  different  kinds  of  tests.  First,  a comprehensive  set  of 
gas  mixtures  must  be  defined  to  allow  a full  kinetic  characterization  of  the  test  material.  Sec- 
ond, a gas  mixture  or  a set  of  gas  mixtures  must  be  defined  for  experimental  screening  or  ranking 
of  materials  similar  to  those  which  have  received  the  full  characterization  treatment.  Although 
it  is  likely  but  not  necessary,  the  screening  gases  and  their  test  conditions  will  be  a subset  of 
the  full  characterization  test  matrix. 

3. 2. 1.1  Rocket  Motor  Environments 

Rocket  motor  environments  were  based  on  three  advanced  MX  propellants,  namely, 

9 XLDB 

9 HTPB 

9 PEG/FEFO 

Representative  elemental  compositions  and  flame  temperatures  are  given  in  Table  2.  For  the 
purpose  of  studying  surface  kinetics  only,  the  elemental  composition  of  the  propellant  needs  to 
be  considered.  The  solid  ASIgOj  does  not  enter  into  the  surface  kinetics  problem  although  it  prob- 
ably contributes  to  surface  erosion  rates.  Table  3 gives  representative  compositions  of  the  propel- 
lant gases  with  all  the  At  and  an  appropriate  amount  of  oxygen  removed  as  At^CTj.  The  ACE/GASKET  pro- 
gram was  used  to  determine  the  concentration  of  gas  species  which  would  exist  at  the  carbon  surface 
for  three  conditions:  (1)  surface  equilibrium,  (2)  very  small  surface  ablation,  and  (3)  a nonreact- 

ing surface  at  typical  surface  temperatures  (2200°K  to  3300°K).  Those  species  with  significant  con- 
centrations would  then  be  candidates  for  reactants  and/or  poisons.  A typical  distribution  of  sur- 
face species  as  a function  of  temperature  for  an  HTPB  propellant  is  shown  in  Figure  8.  This  solution 
represents  the  kinetically  controlled  ablation  of  edge-oriented  pyrolytic  graphite  at  a thruat 


TABLE  2.  REPRESENTATIVE  COMPOSITION  AND  FLAME 
TEMPERATURE  OF  ADVANCED  MX  PROPELLANTS 


Propellant 

Flame  Temperatures  (°K) 
(°F) 

XLDB 

3880 

6524 

HTPB 

3690 

6182 

PEG/FEFO 

3787 

6360 

Mass  Fraction 

H 

2.5 

4.0 

2.6 

C 

13.5 

8.4 

12.5 

N 

24.0 

9.0 

23.0 

0 

39.5 

40.0 

37.9 

F 

- 

- 

1.5 

At 

18.5 

17.6 

18.5 

Ct 

2.0 

21.0 

4.0 

TABLE  3.  PROPELLANT  GAS  COMPOSITION 
(AigOj  REMOVED) 


Propellant 

XLDB 

HTPB 

PEG/FEFO 

Mass  Fraction 
II 

3.8 

6.1 

4.0 

C 

20.8 

’2.7 

19.2 

N 

36.9 

13.6 

35.4 

0 

35.5 

35.2 

32.9 

F 

- 

- 

2.3 

Ct 

3.0 

32.4 

6.2 

23 


2000 


Figure  8.  Typical  surface  gas  composition  at  throat  for  c plane  P6. 
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pressure  of  39.4  atm  from  the  ACE/GASKET  calculations.  Those  species  considered  as  possibly  signif- 
icant reactants  (molar  concentrations  greater  than  0.1  percent)  are: 

• CO 

• h2o 

• h2 

• n2 

• co2 

• HCt 

• HF  (HF  not  a specie  for  the  HTPB  solution) 

It  should  be  noted  that  other  species,  such  as  C£  and  H,  appear  in  representative  amount  and 
may  also  be  important.  Still  other  species,  such  as  0 and  OH,  though  present  only  in  small  quanti- 
ties, may  have  very  fast  reaction  rates.  The  concentrations  of  these  latter  species  decrease 
rapidly  as  the  surface  temperature  increases.  In  fact,  at  typical  rocket  motor  temperatures,  these 
concentrations  are  too  low  to  cause  any  significant  amount  of  carbon  removal.  Atomic  hydrogen  has 
been  shown,  at  least  in  one  case,*  to  react  slower  than  H2  and  since  the  concentrations  of  H2  are 
an  order  of  magnitude  greater  than  that  of  H,  ablation  due  to  the  latter  will  probably  be  insignif- 
icant. Cl,  a halogen,  is  a potential  poison;  however,  there  is  no  firm  evidence  for  this  behavior. 
Thus,  the  species  of  interest  are  those  previously  listed. 

3. 2. 1.2  Important  Surface  Reactions 

The  most  probable  surface  reactions  can  be  identified  by  considering  the  available  species, 
the  possible  reactions  with  carbon,  and  the  equilibrium  constant  for  each  reaction.  (The  equilib- 
rium rate  serves  as  an  upper  limit  to  the  surface  kinetic  rate.)  Of  the  reactions  considered,  only 
the  following  have  sufficiently  large  equilibrium  constants  in  the  temperature  range  of  interest: 

C*  + h2o  -*■  CO  + h2 

C*  + C02  2 CO 

2C*  + H2  - C2H2+ 

- M-r- 

Personal  communication.  Professor  0.  Rosner,  Yale  University. 

+ Although  C2H2  does  not  appear  in  the  list  of  gas  species,  the  reactions  should  not  be  ruled  out. 
Hydrogen  is  present  in  large  concentrations  (approximately  25  percent  by  mole)  and  the  C?H?  coming 
off  the  surface  may  be  eliminated  by  gas  phase  reactions.  £ i 
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The  reactions  of  carbon  with  CO,  ^ HCJt,  and  HF  are  not  considered  significant  since  the  equilib- 
rium formation  rates  are  too  small.  However,  they  may  have  inhibitor  properties. 

The  mechanism  of  poisoning,  or  inhibition  of  surface  reactions,  is  basically  one  of  active 
site  competition.  That  is,  a poison  specie  may  occupy  an  active  lattice  site  and  thus  prevent  a 
reactant  from  occupying  that  site.  Of  the  seven  species  listed  as  possible  poisons,  only  Ng  will 
not  be  considered  because  it  appears  to  be  inert  as  far  as  surface  kinetics  are  concerned.*  Al- 
though H2O  and  COg  readily  react  with  the  carbon  surface  they  also  occupy  lattice  sites  and  thus, 
in  that  sense,  are  poisons  for  each  other. 

3. 2. 1.3  Arc  Plasma  Generator  Limitations 

The  design  and  operation  of  the  arc  plasma  generator  imposes  restrictions  on  the  choice  and 
the  use  of  certain  test  gases.  There  are  two  basic  areas  of  concern;  first,  the  effect  of  a partic- 
ular gas  on  the  vital  components  of  the  APG  (cathode,  anode,  constrictor  segments,  etc.)  and  second, 
the  stability  of  the  arc  when  operating  with  a particular  gas  or  combination  of  gases.  The  situa- 
tion is  further  complicated  by  the  desire  to  produce  a test  gas  at  the  highest  possible  temperature. 
This  generally  requires  arc  heating  of  the  largest  possible  portion  of  the  total  test  gas  flow  to 
the  maximum  temperature  achievable,  i.e.,  maximum  energy  input,  while  minimizing  the  energy  losses 
to  the  cooled  walls  of  the  APG. 

With  the  design  of  the  APG  currently  being  used,  it  is  necessary  to  avoid  injecting  any 
oxidizing  species  into  the  arc  heater  as  the  primary  gas.  This  is  due  to  the  tungsten  material 
used  in  the  cathode,  which  when  rapidly  removed  through  oxidation  processes  can  both  limit  APG  run 
times  to  the  order  to  seconds  and  cause  catastrophic  failure  of  the  arc  heater.  The  normal  solu- 
tion employed  is  the  injection  of  such  gases  several  constrictor  duct  diameters  downstream  of  the 
cathode.  This  has  been  highly  successful  when  the  required  test  gas  is  simulated  air,  using  in- 
dividually injected  nitrogen  and  oxygen.  However,  in  the  case  of  propellant  simulation,  there  is 
an  additional  problem.  The  reactive  nature  of  the  base  species,  hydrogen,  which  for  reasons  of  arc 
efficiency  and  maximum  power  input  is  the  arc  heated  gas,  requires  the  injection  of  oxidizing  spe- 
cies downstream  of  the  arc  heater  portion  of  the  APG,  in  the  plenum  section  (see  Mgure  1).  fins 
is  primarily  due  to  the  combustion  induced  turbulence  which  adversely  affects  the  stability  of  the 
arc,  resulting  in  failure  of  the  constrictor  segments.  Therefore,  the  primary  or  arc-heated  gases 
must  be  either  inert  or  nonoxidizing;  the  remainder  of  the  test  gases  required  to  make  up  the 


it 

Kinetic  rate  data  in  Reference  6 substantiates  the  inert  behavior  of  Ng. 
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propellant  simulation  are  injected  in  the  plenum  section.  This  results  in  lower  overall  AFG  effi- 
ciencies due  to  the  portion  of  the  test  gas  which  is  not  directly  arc  heated  and  the  losses  to  the 
plenum  section  from  both  the  arc  heated  primary  gases  and  the  exothermic  reactions  which  take  place 
in  the  plenum.  The  net  effect  is  lower  test  gas  temperature  and  hence  lower  test  sample  surface 
temperature. 

The  "normal"  APG  limits  of  pressure,  current  and  power  input  must  also  be  considered.  These, 
in  general,  are  less  severe  than  those  discussed  above  and  typically  can  be  accommodated  through 
arc  heater  and  power  supply  configuration  changes.  It  should  be  noted  that  this  is  especially  true 
with  hydrogen,  which  is  very  sensitive  to  the  gas  injection  configuration  and  arc  heater  constric- 
tor length.  The  penalty  for  use  of  an  improper  configuration  is  usually  very  unstable  arc  opera- 
tion. 

3. 2.1. A Potential  Test  Gases 

As  described  in  Section  3. 2. 1.2,  the  potentially  important  reactants  are  H^,  H.,0,  and  CO,,. 

In  an  APG,  various  concentrations  of  these  gases  can  be  mixed  and  reacted  to  form  the  test  stream. 

By  judicious  selection,  various  reactants  and  poisons  could  be  isolated  in  a systematic  manner  so 
that  appropriate  reaction  rate  constants  could  be  determined.  Possible  test  gas  mixtures  are  shown 
in  Table  4.  These  gases  are  separated  into  three  groups,  reactions  which  include  H2,  H20,  C02 , 
and  CO,  reactions  with  these  gases  and  HC£,  and  reactions  with  HF  in  lieu  of  HCT  The  surface  reac- 
tion designations  are  shown  in  Table  5.  The  number  of  test  gases  to  be  used  in  a material  charac- 
terization test  matrix  would  be  selected  as  a subset  of  the  gases  tested  in  Table  4.  This  selection 
will  be  based  upon  a trade-off  between  the  degree  to  which  a particular  reactant  (or  poison)  can  be 
isolated  and  the  operating  limitations  of  the  APG.  Note  that  the  exhaust  gas  composition  is  only 
representative  and  that  all  gases  that  contain  CO  will  also  have  C02  ii  small  quantities.  At  high 
temperatures,  it  is  not  possible  to  have  large  concentrations  of  C02  in  the  presence  of  H2  since 
the  preferred  species  would  be  H.,0  and  CO. 

3.2.2  Test  Gas  Selection 

The  gases  shown  in  Table  4 that  contain  HCi  and  HF  require  special  toxic  gas  handling  sys- 
tems. The  current  Aerothenn  APG  facility  is  equipped  to  handle  HC)t,  although  a number  of  nontrivial 
additions  are  required  before  HF  can  be  used.  For  this  reason,  test  gases  that  contained  HF  were 
eliminated  during  the  test  gas  selection. 
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'Hass  ratio  of  plenura  injected  to  arc  heated  gases  (approxirate) 


TABLE  5.  SURFACE  REACTIONS 


Reaction  No. 

HHi 

2 

mgg/g 

■ffyipf 

! 6 
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3. 2. 2.1  Test  Gas  Evaluation 


Test  gases  1 through  8 have  been  evaluated  under  a wide  variety  of  APG  conditions  using  both 
water-cooled  calorimeters  and  carbon  test  sections.  These  tests  clearly  show  that  test  gases  2 and 
4 resulted  in  anomalous  heating  conditions.  The  probable  cause  can  be  defined  by  considering  the 
schematic  of  the  APG  shown  in  Figure  1.  In  normal  operations,  H 2 or  an  inert  gas  such  N,,,  Ar,  or 
He  is  used  as  the  arc  heated  column  and  all  other  gases  are  injected  between  the  arc  column  and  the 
plenum  chamber.  If  we  consider  test  gas  2 as  an  example,  the  ratio  of  injection  gas  (02)  to  arc 
column  gas  (H2)  is  4/9.  However,  mixing  of  the  two  gases  will  be  dependent  upon  their  relative  mass 
rates.  A simple  conversion  shows  that  the  relative  mass  of  injection  gas  to  arc  column  gas  is  ap- 
proximately 7/1.  It  was  originally  anticipated  that  combustion  induced  turbulence  could  result  in 
adequate  mixing  in  the  plenum,  however,  measured  data  suggested  a high  concentration  of  low  enthalpy 
injection  gases  near  the  walls  of  the  test  section.  This  rather  poor  mixing  of  arc  heateo  and  in- 
jection gases  made  the  test  data  impossible  to  adequately  analyze.  Subsequent  trial  and  error  exper- 
imentation showed  that  ratios  of  injected  gas  to  arc  heated  gas  of  less  than  5 (by  mass)  would  re- 
sult in  adequate  plenum  chamber  mixing.*  Thus  test  gas  number  5 was  also  eliminated. 

Experimentation  with  test  gases  4 and  10  revealed  a second  difficulty.  The  kinetic  reaction 
rates  of  C02  with  H2  are  much  slower  than  those  of  H2  with  02-  In  fact,  some  simple  kinetic  calcu- 
lations revealed  that  the/e  was  insufficient  residence  time  in  the  plenum  chamber  to  attain  thermo- 
chemical  equilibrium.  Thus,  all  gases  which  would  normally  inject  C0?  would  be  replamj  by  an  equi- 
valent combination  of  CO  and  02< 

From  the  above  discussion,  HCi  test  gases  12  and  13  can  be  eliminated  outright,  however,  oas 
number  11  can  be  made  acceptable  by  reducing  the  relative  moles  of  HCt  from  8 to  2. 

3. 2. 2. 2 Recommended  Test  Gases  for  Full  Characterization  Studies 

Based  upon  the  discussion  in  Section  3. 2. 2.1,  the  test  gases  for  material  characterization 
studies  were  reduced  to  the  subset  shown  in  Table  6.  Note  that  C02  was  not  used  as  an  injection 
gas  and  that  it  was  replaced  by  an  equivalent  quantity  of  02  and  CO,  Note  also  that  (IF  gases  were 
not  included  since  the  advisability  of  testing  with  HF  has  not  yet  been  assessed. 


★ 

It  is  assumed  that  at  least  one  of  the  injected  gases  will  be  02  so  that  there  will  be  ;ombustion 
induced  turbulence. 
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TABLE  6.  RECOMMENDED  MATERIAL  CHARACTERIZATION  TEST  GASES 


No. 

• APG  Input  Gases 
Relative  Moles 

R+ 

Equilibrium  Exhaust  Gases 
Relative  Moles 

Surface 
Reactions 
(see  Table  5) 

* 

H2 

°2 

CO 

. 

Ar 

HC£ 

H2 

h2o 

CO 

Ar 

HC£ 

i 

1 

0 

1 

) 

2 

2 

1 

4 

0.2 

2 

4 

2.4 

3 

8 

1 

2 

6 

2 

1.2.4 

4 

8 

1 

1 

4 

6 

2 

1 

1,2. 3, 4 

5 

2 

1 

2 

4 

0.5 

2 

2 

4 

2.4 

6 

10 

■ 

■ 

■ 

1 

1 

10 

■ 

n 

mm 

1,5 

7 

8 

n 

lill 

I 

■ 

6 

2 

■ 

■ 

n 

1.2,3 

8 

8 

n 

D 

■ 

1 

6 

2 

n 

■ 

H 

1.2, 3, 4,5 

9 

D 

5 

Inert  Gas 

* 

Arc  heated  gas 


Mass  ratio  of  plenum  Injected  to  ere  heated  gases  (approximate) 
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With  the  exception  of  the  HF  inhibitor,  all  other  surface  reactants  are  represented  by  this 
set  of  reactions.  It  is  clearly  not  possible  to  isolate  reactions  other  than  H2  since  oxygen- 
bearing species  (C02,  H20)  will  react  with  solid  carbon  to  form  CO  and  in  gas  phase  equilibrium,  a 
small  quantity  of  C02  will  also  be  present.  The  reactions  shown  in  Table  6 represent  a good  com- 
promise between  the  desire  to  isolate  reactants  and  still  stay  within  the  operating  limitations  of 
the  APG.  An  inert  gas  was  also  included  in  Table  6 to  test  for  shear  removal  affects.  Inert  gas 
tests  were  run  at  the  highest  heating  conditions  compatible  with  APG  limitations.  With  3 tests  for 
5 react'ig  gas  mixtures,  2 tests  for  3 mixtures,  and  1 for  the  inert  gas  test,  a minimum  of  22 
tests  were  required.  Six  additional  tests  ware  plannori  contingency  or  repeat  tests  and  were  per- 

formed as  required.  Thus,  a total  of  28  tests  were  planned  for  each  characterization  material. 

3. 2. 2. 3 Selection  of  Test  Gases  for  Limited  Characterization  Studies 

Since  fewer  samples  were  to  be  used  for  the  limited  characterization  materials  than  the  fully 
characterized  materials,  more  stringent  criteria  were  imposed  on  the  selection  of  these  test  gases. 
Test  gases  should  satisfy  the  following  thermal  performances,  without  exceeding  the  performance 
limitations  of  the  arc  plasma  generator: 

• High  condition  - q = 1600  Btu/ft2sec,  run  time  = 80  sec 

c • w • 

• Medium  condition  - q_  =1200  Btu/ft2sec,  run  time  = 100  sec 

Wi 

• Low  condition  - q^  w = 1000  Btu/ft2sec,  run  time  = 120  sec 

Furthermore,  these  test  gases  should  also  be  able  to  characterize  the  surface  reactions  as  shown  in 
Table  5. 

The  test  gases  were  selected  statistically  from  the  results  of  the  full  characterization 
studies.  The  probable  error  for  each  test  gas  with  respect  to  the  least  square  curve  fit  for  all 
test  gases  was  determined  ir.  the  process  of  correlating  kinetic  data.  The  test  gas  with  the  least 
probable  error  was  considered  to  be  the  potential  candidate  for  limited  characterization  studies. 

The  results  of  the  statistical  evaluation  of  test  gases  for  each  generic  material  are  shown  in 
Figures  9 through  11. 

Turther  analyses  are  required  to  reach  the  final  set  of  test  gases  for  limited  characteriza- 
tion studies.  The  reason  is  that  test  gases  which  were  selected  statistically  may  serve  a similui 
purpose  in  characterizing  surface  reacticns.  Such  redundancy  should  be  avoided  if  tie  number  uf 
data  to  be  taken  are  limited.  For  example,  in  determining  test  gases  for  carbon/carbon  material;,, 
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from  a minimum  error  standpoint,  test  gases  2,  3,  and  4 (see  Figure  11)  should  be  selected.  However, 
they  do  not  provide  enough  isolation  of  specific  reaction  to  be  a good  screening  gas  set.  To  arrive 
at  the  final  screening  gas  set,  gas  2,  which  is  similar  to  gas  3,  was  replaced  by  gas  1 in  order  to 


get  an  isolation  of  H2  reaction.  Similarly,  gas  4 was  replaced  by  gas  7 in  order  to  assess  HC2.  in- 
hibition. As  for  CO  inhibition,  a test  gas  with  an  arbitrary  amount  of  CO  is  not  required  for  all 
gas  systems  with  O2.  Thus,  the  final  test  gases  selected  for  limited  characterization  of  carbon/ 
carbon  are  1,  3,  and  7.  Similar  arguments  were  used  for  other  generic  materials.  The  resultant 
sets  of  limited  characterization  test  gases  are  shown  in  Table  7. 

3.3  CARBON  MATERIALS  SURFACE  KINETICS  CORRELATION  PROCEDURE 

The  ability  to  achieve  a successful  empirical  formula  which  describes  the  reactivity  of  a 
carbon  surface  with  propellant  gases  depends  strongly  on  the  selection  of  the  correlation  function 
and  the  kinetically-controlled  reactions.  Based  on  previous  Aerotherm  experience  (Reference  2), 
the  surface  kinetics  of  pyrolytic  graphite  were  correlated  by  applying  the  Langmuir-Hinshelwood 
model  and  assuming  that  the  following  reactions  were  kinetically  controlled: 

C*  + h2o  ->  h2  + CO 

C*  + C02  - 2C0 

2C*  H2  C2H2 

These  kinetically  controlled  reactions  can  be  inhibi tea  by  H20,  C02«  CO,  H2,  HCi,,  and  HF. 

Since  it  is  logical  to  conclude  that  all  carbon  materials  have  similar  kinetic  behavior,  the 
same  kinetic  model  was  adopted  to  correlate  all  the  fully  characterized  materials.  In  addition, 
the  mass  consumption  rate  due  lo  sublimation  is  included  when  surface  temperatures  exceeding  6000°F 
are  expected. 

Carbon  surface  kinetics  based  on  the  Langmuir-Hinshelwood  model  is  given  by: 


..TnYEim«/P. . _l 

1 W 1 Kn 
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pi 


prod. 
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(3) 


where  the  subscript  i denotes  each  of  the  reactants,  H20,  C02,  and  H2-  The  numerator  of  Equation 
(3)  describes  the  surface  reaction  with  gas  phase  species,  and  the  denominator  describes  the 
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TABLE  7.  RECOMMENDED  TEST  GASES  FOR  LIMITED  CHARACTERIZATION  STUDIES 


Test  Gas 
No. 

Modified 

Pyrolytic  Graphite 

Carbon/ 

Carbon 

Surface  Reactions 
(See  Table  5) 

1 

✓ 

/ 

/ 

- 

3 

/ 

/ 

/ 

1,  2,  4 

7 

/ 

/ 

1,  2,  3 

8 

/ 

1,  2,  3,  4,  5 
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surface  coverage  by  gas  phase  species.  As  can  be  seen,  the  surface  reaction  can  be  retarded  if  the 
surface  is  either  desorption  controlled  or  is  poisoned  by  species  such  HCji  or  HF. 

The  correlation  function  requires  further  manipulation  before  applying  an  Aerotherm  least 
squares  fitting  program  to  determine  the  coefficients.  Usually,  the  following  assumptions  are  made 
to  simplify  the  correlation  function:  the  reverse  rate  is  negligible  compared  to  the  forward  rate, 
and  H20  and  C02  surface  reactions  have  the  same  activation  energy  and  the  same  inhibiting  effect  by 
the  gas  phase  species.  These  two  assumptions  can  be  easily  removed  if  sufficient  kinetic  data  are 
available.  With  these  two  assumptions,  Equation  (3)  can  be  rewritten  as: 


MC  = B^e 
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where 


°1  = 


°2 


1 ♦ <AP)H20  + (AP)^  1 (AP)C0  4 (AP)Hz 


1 + <ap>h2o + (ap>co2  + <AP>co  + <Af\ 


H20,  co2 


(5) 


(6) 


The  coefficients  which  must  be  determined  from  the  data  analysis  are  A,  B,  n,  and  E. 

The  data  required  to  correlate  carbon  material  surface  kinetics  were  obtained  from  arc  plasma 
generator  (APG)  and  motor  firing  tests.  The  APG  data  provides  information  on  total  surface  reces- 
sion, surface  temperature,  experimental  time,  and  reactive  gas  chamber  conditions.  These  data  have 
a relatively  low  surface  temperature  (3000°R  - 55G0°R)  and  edge  pressure  (2-7  atm)  range,  but  can 
be  utilized  to  characterize  H20  and  C02  surface  reactions.  The  motor  firing  data  provides  similar 
information  except  that  the  surface  temperature  is  unknown.  However,  this  surface  temperature  can 
be  estimated  with  a semi-infinite  slab  analytical  solution.  The  motor  data  generally  have  higher 
surface  temperatures  (5500°R  - 6000°R)  and  edge  pressures  (30  - 100  atm)  than  the  APG,  and  can  he 
used  to  characterize  the  H2  surface  reaction. 

These  raw  data  required  further  reduction  before  use  to  determine  the  kinetic  coefficients. 
This  data  reduction  procedure  is  described  in  the  next  section.  The  reduced  data  contains  informa- 


tion on  mc,  T , Pg,  and  surface  composition. 
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Given  a sufficient  number  of  data  points  with  m , T , P„,  and  surface  composition  known,  the 
coefficients  in  Equations  (4),  (5),  and  (6)  can  be  determined  using  an  Aerotherm  least  squares  fit- 
ting program.  In  principle,  the  coefficients  are  adjusted  systematically  to  minimize  the  residual 
of  the  least  squares  curve  fit  through  a plot  of  log  (MCi ) versus  TWi-  The  systematic  adjustments 
of  the  coefficients  are  based  upon  the  method  of  steepest  descent  in  the  vector  space  of  R(A^ , , 

E^).  Because  of  the  large  number  of  coefficients  and  the  fact  that  there  is  no  assurance  that  the 
residual  R is  unimodal,  there  may  be  several  local  minimums.  Whether  or  not  the  true  minimum  is 
"discovered"  by  the  optimization  search  depends  largely  upon  the  initial  prescribed  values  of  the 
coefficients. 

Because  the  current  data  presumes  no  advance  knowledge  of  the  reaction  kinetics  of  the  fully 
characterized  graphite  materials,  the  starting  point  of  the  optimization  search  was  arbitrarily 
started  with  the  coefficients  determined  for  c plane  oriented  pyrolytic  graphite  (see  Reference  6). 

Since  the  surface  temperatures  for  the  motor  firings  were  approximated  with  a best  guess 
solution,  it  was  necessary  to  analyze  the  data  in  the  following  steps  to  arrive  at  the  final  kinetic 
coefficients. 

1.  A correlation  of  the  APG  data  and  the  motor  firing  data  (based  on  approximate  surface 
temperatures)  was  obtained  using  the  least  squares  optimum  seeking  code. 

2.  The  coefficients  were  input  into  GASKET  for  each  motor  firing  data  point. 

3.  CMA  solutions  were  obtained  using  the  GASKET  generated  surface  thermochemistry  tables. 

This  yielded  a predicted  surface  response  for  the  fully  characterized  graphite  material. 

4.  The  CMA  results  were  used  to  update  the  motor  firing  data  points.  An  arbitrary  choice 

was  made  to  use  tne  data  at  the  halfway  point  of  the  firing  duration. 

5.  A second  correlation  was  obtained  from  the  APG  data  and  the  updated  motor  firing  data. 
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SECTION  4 


ANALYSIS  OF  APG  TEST  RESULTS 

Arc  plasma  generator  test  data  are  shown  in  Table  8 for  the  full  characterization  materials. 
These  data  are  not  in  a form  compatible  with  the  surface  kinetics  correlation  procedure;  hence,  fur- 
ther data  reduction  was  required.  A procedure  was  set  up  to  reduce  the  APG  test  results  (PQ , H , , 

Tw>  AS,  experimental)  ^or  t*iese  kinetic  correlations.  This  procedure  is  as  follows: 

• Determine  the  boundary  layer  edge  condition 

• Evaluate  the  heat  and  mass  transfer  coefficients 

• Determine  the  open  system  surface  state  chemical  equilibrium  solution  using  measured 
carbon  consumption  rate  and  surface  temperature 

The  net  effect  is  a relationship  between  carbon  consumption  rate,  surface  temperature  and  partial 
pressures  of  gases  adjacent  to  the  ablating  surface.  The  reduced  data  before  correlation,  for  all 
the  full  characterization  carbon  materials,  are  presenteo  in  Table  9. 

APG  test  data  for  the  limited  characterized  materials  are  shown  in  Tables  10  through  12. 
Because  no  correlations  were  attempted,  no  further  reduction  of  this  data  was  required. 

4.1  CARBON  CONSUMPTION  RATE 

The  carbon  consumption  rate  can  be  calculated  from  the  expression: 

”c  ' 12  ” T <7> 

where  p is  the  carbon  density,  As  is  the  total  measured  recession  (in.),  and  0 is  the  actual  reaction 
time.  However,  the  determination  of  mc  is  not  straightforward  because  u is  inknown  without  a prior 
knowledge  of  the  surface  kinetics. 

One  way  to  estimate  the  reaction  time  is  to  assume  the  tested  carbon  material  has  a kinetic 
response  similar  to  edge  PG.  Substitution  of  the  measured  surface  temperature  history  into  the 
edge  PG  surface  kinetics  expression  allows  a recession  rate  history  to  be  obtained.  The  reaction 
time  is  approximately  the  time  interval  between  the  final  time  and  the  time  where  the  recession  rate 
suddenly  increases. 
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TABLE  8.  ARC  PLASMA  GENERATOR  DATA  FOR  FULL  CHARACTERIZATION  MATERIALS 


Test 

Number 

Model 

Number 

Test 

GdS 

Po 

(atm) 

Ho 

(Btu/lbm) 

G-90 

2809-1 

178  C 

1 

3.25 

35,700 

2810-1 

169  C 

1 

3.48 

46,230 

2811-1 

177  C 

1 

1 

3.55 

55,500 

2813-1 

196  C 

3 

5.20 

8,470 

2814-1 

176  C 

3 

5.18 

10,930 

2315-1 

168  C 

3 

4.825 

14,630 

2817-1 

167  C 

4 

5.13 

6,300 

2818-2 

194  C 

4 

5.10 

6,970 

2820-2 

195  C 

2 

5.00 

2,590 

2821-2 

187  C 

2 

4.95 

2,060 

2823-3 

163  C 

5 

6.25 

1 ,440 

2823-5 

172  C 

5 

6.25 

980 

2826-2 

186  C 

6 

5.60 

14,950 

2828-1 

171  C 

6 

5.70 

12,750 

2329-2 

185  C 

6 

5.40 

10,210 

2831-2 

162  C 

7 

4.32 

12,880 

2832-2 

170  C 

7 

3.98 

12,480 

2833-2 

161  C 

7 

4.00 

10,260 

2836-2 

190  C 

8 

5.30 

12,340 

2838-2 

184  C 

8 

5.05 

12,140 

ATJ 

2585-1 

003  C 

1 

4.48 

64,000 

2589-2 

004  C 

3 

5.50 

8,190 

2591-1 

005  C 

3 

5.50 

13,700 

2592-1 

006  C 

3 

5.70 

16,700 

2593-1 

007  C 

1 

4.20 

55,130 

2594-1 

009  C 

3 

5.00 

9,520 

2597-1 

010  C 

4 

5.85 

9,300 

2607-1 

008  C 

2 

7.70 

911 

2609-1 

Oil  C 
_ 

2 

4.80 

4,030 

AS 

e 

exp 

(mils) 

(sec) 

TABLE  8.  Continued 


Model 

Number 

Test 

Gas 

Po 

(atm) 

Ho 

(Btu/lbin) 

Tw 

(°R) 

m 

(lb/ft2sec) 

AS 

(mils) 

0 

exp 

(sec) 

mmmm 

2617-1 

2 

5.35 

2,?30 

3860 

0.00810 

39.20 

38.0 

2619-1 

014  C 

5 

6.40 

1,800 

4000 

0.00733 

34.70 

44.5 

2619-2 

015  C 

5 

6.40 

1,660 

4200 

0.00896 

43.40 

43.0 

2620-2 

016  C 

5 

7.15 

1,460 

3580 

0.00395 

10.00 

24.0 

262C-3 

017  C 

1 

905 

3720 

0.00540 

26.40 

44.5 

2815-3 

019  C 

B 

5.10 

14,700 

4950 

0.01325 

52.50 

40.0 

2817-3 

028  C 

5.30 

5,560 

4500 

0.01998 

82.50 

48.0 

2818-1 

029  C 

5.00 

6,430 

4450 

0.010991 

48.40 

45.0 

2818-4 

030  C 

H 

5.25 

7,000 

4800 

0.00903 

50.70 

54.0 

2820-1 

027  C 

2 

5.00 

2,550 

4900 

0.01075 

50.30 

55.0 

2823-1 

020  C 

5 

5.45 

1,470 

4750 

0.008819 

53.40 

59.0 

2826-1 

021  C 

6 

5.70 

14,710 

5300 

0.004192 

15.00 

37.0 

2827-1 

022  C 

6 

5.63 

13,250 

5050 

0.003292 

14.60 

57.0 

2829-3 

023  C 

6 

5.16 

[10,280] 

4600 

0.0008054 

2.66 

2831-3 

024  C 

1 

4.30 

13,280 

4850 

0.004814 

26.50 

59.0 

2832-3 

025  C 

B 

4.30 

11,810 

4850 

0.004723 

27.30 

60.0 

2834-1 

031  C 

B 

4.05 

9,290 

4500 

0.004704 

24.60 

38.5 

2837-1 

018  C 

H 

4.85 

12,940 

5000 

0.004872 

20.25 

48.5 

2838-3 

026  C 

8 

4.95 

12,640 

4850 

0.005450 

30.00 

63.0 

15%  SiC 

PG 

2668-1 

104  C 

I , 

6.60 

11,780 

4850 

0.01014 

26.00 

32.0 

2669-1 

105  C 

B 

6.20 

6,209 

4000 

0.00369 

9.45 

31.0 

2672-1 

106  C 

B 

3.35 

36,590 

4250 

-0.0002645 

-0  52 

33.0 

2674-1 

107  C 

i 

3.15 

31,370 

4050 

0.0000272 

0.05 

34.0 

2676-1 

108  C 

6.08 

2,650 

3720 

0.002111 

2 67 

32.0 

2678-1 

109  C 

5.62 

8,550 

4150 

0.004179 

11.25 

34.0 

2679-1 

110  C 

5.50 

11,880 

4775 

0.01493 

31.90 

31.0 

2680-1 

111  C 

5.65 

15,470 

4900 

0.01525 

40.40 

i 32.0 

i 

2705-1 

114  C 

7 

4.32 

9,890 

3760 

0.0002457 

1.05 

i 58.0 
1 

TABLE  8.  Continued 


Model 

Number 

Test 

Gas 

Po 

(atm) 

Ho 

(Btu/lbm) 

Tw 

(°R) 

m 

(lb/ft^sec) 

AS 

(mils) 

6exp 

(sec) 

15*  SiC  PG  (Concluded) 

2707-1 

115  C 

4.58 

12,180 

3900 

0.001577 

5.73 

55.5 

2709-1 

116  C 

H 

4.60 

14,140 

4350 

0.005113 

23.60 

59.0 

2711-1 

117  C 

V 

5.13 

7,860 

3960 

0.008093 

16.60 

30.0 

2712-1 

118  C 

V 

5.53 

8,430 

4000 

0.006700 

27.20 

57.5 

2715-1 

119  C 

8 

4.80 

10,040 

4200 

0.004235 

20.45 

60.0 

2716-1 

120  C 

8 

4.48 

10,190 

4150 

0.003978 

27.20 

90.0 

2719-1 

121  C 

3 

5.28 

11,680 

4300 

0.01086 

39.90 

44.0 

2721-1 

122  C 

3 

5.40 

9,730 

3920 

0.007406 

31.65 

55.0 

2723-1 

123  C 

1 

3.30 

31 ,580 

4400 

0.001009 

4.40 

57,0 

2725-1 

124  C 

1 

4.00 

45,450 

5050 

0.004020 

13.40 

42.5 

2727-1 

125  C 

4' 

5.85 

7,180 

3780 

0.002808 

9.00 

50.5 

Pyrocarb  901 

1 

2809-2 

140  C 

1 

3.25 

36,630 

4350 

0.0005241 

3.08 

55.0 

2810-2 

142  C 

1 

3.55 

49,320 

5000 

0.001997 

16.75 

52.5 

2811-2 

158  C 

1 

3.50 

55,430 

5358 

0.005077 

21.70 

59.0 

2813-2 

149  C 

3 

5.00 

8,880 

4150 

0.009770 

42.80 

36.5 

2814-2 

150  C 

3 

5.00 

10,300 

4700 

0.01116 

44.70 

37.5 

! 

2815-2 

159  C 

3 

5.00 

14,200 

5100 

0.01147 

51.50 

39.5  j 

2817-2 

151  C 

4 

5.00 

6,140 

4700 

0.009489 

58.80 

56.0  : 

1 

2818-3 

160  C 

4 

5.00 

7,060 

5100 

0.01329 

63.90 

45.0  | 

2820-3 

137  C 

2 

5.00 

2,540 

5200 

0.008035 

49.80 

58.0  j 

2821-1 

146  C 

2 

4.90 

2,020 

4950 

0.009405 

59.30 

59.0 

2823-2 

138  C 

5 

6.20 

1,530 

5000 

0.02257 

72.30 

60.0 

2823-4 

147  C 

5 

5.90 

950 

4700 

0.01856 

58.50 

59.0 

2825-1 

139  C 

6 

4.35 

14,850 

5750 

0.01402 

47.20 

38.0 

2826-3 

156  C 

6 

5.50 

15,910 

5500 

0.006820 

26.60 

50.0 

2828-2 

157  C 

6 

5.65 

13,250 

5500 

0.005640 

22.60 

50.0 

2829-1 

152  C 

6 

5.50 

10,420 

4800 

0.0006405 

1.54 

42.0 

44 


TABLE  8.  Concluded 


mm 

Model 

Number 

Test 

Gas 

Po 

(atm) 

Ho 

(Btu/lbm) 

Tw 

(°R) 

m 

(lb/ft2sec) 

AS 

(mils) 

e 

exp 

(sec) 

Pyrocarb  901  (Concluded) 

2831-1 

143  C 

7 

4.75 

14,010 

4900 

C. 008084 

40.60 

58  5 

2832-1 

134  C 

7 

4.15 

12,940 

5000 

0.008564 

33.40 

45.0 

2833-1 

153  C 

7 

3.82 

'9,820 

4950 

0.007745 

24.00 

48.0 

2836-1 

144  C 

8 

5.30 

11,600 

5000 

0.009186 

37.30 

55.0 

2838-1 

154  C 

8 

5.00 

12,170 

5050 

0.007721 

39.60 

60.5 

45 


46 


I 


48 


50 


TABLE  10.  ARC  PLASMA  GENERATOR  DATA  FOR  LIMITED  CHARACTERIZATION  MATERIAL  (CARBON/ CARBONS) 


Test 

Number 

Nodel 

Number 

Test 

Gas 

Po 

(atm) 

Ho 

(Btu/lbm) 

Tw 

(°R) 

m 

(lb/ft2sec) 

AS 

(mils) 

6 

(sec) 

903  Pyr< 

)carb 

2851-1 

218  S 

1 

3.60 

55,470 

5500 

0.003259 

16.875 

49.0 

2855-5 

216  S 

1 

3.50 

46,630 

5300 

0.002393 

22.00 

87.0 

2860-4 

219  S 

B 

3.33 

36,820 

4800 

0.00046 

4.76 

98.0 

2863-1 

220  S 

B 

5.15 

13,100 

5250 

0.00975 

44.30 

43.0 

2867-6 

214  S 

B 

4.93 

7,480 

4550 

0.00825 

40.54 

46.5 

2869-1 

223  S 

7 

4.05 

12,490 

4800 

0.00467 

32.78 

66.5 

2872-5 

217  S 

7 

4.10 

9,710 

4750 

0.00454 

35.55 

74.0 

903  HD 

2854-1 

250  S 

1 

34,700 

4700 

0.000303 

2.9C 

99.0 

2856-3 

243  S 

1 

44,280 

5200 

0.001258 

10.50 

82.0 

2858-2 

247  S 

1 

3.48 

53,980 

5600 

0.00309 

17.00 

67.0 

2864-5 

244  S 

3 

11,740 

5050 

0.01171 

50.90 

45.0 

2865-3 

246  S 

3 

6,200 

4200 

0.00604 

36.80 

63.0 

2871-2 

248  S 

7 

4.50 

12,660 

5100 

0.00533 

32.40 

63.0 

2874-1 

249  S 

7 

4.60 

12,840 

4750 

0.000802 

33.32 

43.0 

2876-3 

245  S 

7 

4.18 

9,280 

4800 

0.00398 

30.78 

60.0 

HRX-5125 

2852-1 

229  S 

■1 

3.50 

40,460 

4950 

0.001064 

13.38 

100.0 

2856-1 

232  S 

B 

3.45 

44,900 

5250 

0,001683 

18.40 

87.0 

2859-2 

228  S 

B 

3.65 

54,630 

5600 

0.001684 

12.38 

58.5 

2860-3 

231  S 

B 

3.60 

54,030 

5650 

0.001751 

12.98 

69.0 

2863-2 

230  S 

3 

4.80 

12,930 

5400 

0.01156 

52.30 

36.0 

2870-1 

225  S 

7 

4.12 

12,080 

4950 

0.00525 

37.60 

57.0 

2875-1 

226  S 

7 

4.35 

9,760 

4850 

0.00391 

29.70 

60.5 

HRX-5875 

2857-3 

239  S 

m 

3.50 

50,270 

5250 

0.002934 

17.89 

59.0 

2858-1 

238  S 

B 

3.48 

50,880 

5450 

0.003903 

23.80 

59.0 

2864-3 

241  S 

V J 

5.0 

11,740 

5150 

0.004624 

16.25 

34.0 

2870-2 

240  S 

7 

4.08 

11,830 

5000 

0.00420 

25.60 

59.0 

51 


TABLE  10.  Concluded 


Test 

Number 

Model 

Number 

Test 

Gas 

Po 

(atm) 

Ho 

(Btu/lbm) 

Tw 

(°R) 

m 

(lb/ft2sec) 

6s 

(mils) 

B 

MDAC  3-D  C/C 

2854-2 

199  S 

1 

3.40 

34,120 

4300 

0.0001004 

0.958 

99.0 

2861-1 

198  S 

1 

3.68 

53,470 

5450 

0.0007621 

4.625 

58.0 

2864-4 

200  S 

3 

5.15 

11 ,520 

5050 

0.01255 

53.71 

44.5 

2866-1 

202  S 

3 

5.12 

6,100 

4300 

0.002574 

14.85 

60.0 

2870-3 

197  S 

7 

4.35 

11,910 

5100 

0.00469 

26.62 

59.0 

52 


TABLE  11.  ARC  PLASMA  GENERATOR  DATA  FOR  LIMITED  CHARACTERIZATION  MATERIAL  (BULK  GRAPHITES) 


Test 

Number 

Model 

Number 

Test 

Gas 

Po 

(atm) 

Ho 

(Btu/lb) 

Tw 

(°R) 

m 

(lb/ft2sec) 

AS 

(mils) 

e 

(sec) 

PO  3 

2584-2 

034  S 

1 

3.68 

41,900 

4450 

0.0002138 

i .33 

59.0 

2585-2 

035  S 

1 

4.55 

61,900 

3200 

0.002931 

13.60 

44.0 

2585-3 

036  S 

1 

4.60 

52,200 

— 

0.001810 

9.54 

50.0 

2589-1 

037  S 

3 

5.35 

7,950 

3700 

0.00946 

46.90 

47.0 

2592 -i. 

038  S 

3 

5.95 

16,350 

5000 

0.01027 

39.20 

36.2 

2596-1 

039  S 

4 

5.55 

6,570 

4100 

0.00763 

37.80 

47.0 

2617-2 

068  S 

2 

5.40 

2,680 

4200 

0.00808 

39.20 

46.0 

2855-3 

070  S 

1 

3.53 

45,520 

5150 

0.05900 

7.64 

83.0 

2878-2 

071  S 

8 

4.95 

10,740 

4800 

0.00241 

17.16 

67.5 

2878-4 

069  S 

8 

4.75 

8,630 

4150 

0.00258 

21.77 

80.0 

ATJ-S 

2617-4 

076  S 

2 

5.70 

2,520 

3800 

0.00582 

1.864 

46.0 

2853-1 

084  S 

1 

3.38 

37,520 

4800 

0.00057 

5.88 

98.0 

2855-2 

077  S 

1 

3.53 

44,960 

5000 

0.0008184 

6.375 

74.0 

2855-4 

085  S 

1 

3.70 

55,100 

5350 

0.0016198 

10.06 

59.0 

2861-3 

082  S 

1 

3.48 

54,040 

5500 

0.001685 

10.64 

60.0 

2863-3 

078  S 

3 

5.15 

12,130 

5150 

0.009112 

46.04 

48.0 

2867-3 

081  S 

3 

4.92 

7,560 

4250 

0.00562 

35.50 

60.0 

2878-1 

079  S 

8 

5.30 

10,080 

4600 

0.00269 

17.29 

61.0 

2878-3 

080  S 

8 

4.75 

8,330 

4100 

0.00244 

20.81 

81.0 

53 


TABLE  12.  ARC  PLASMA  GENERATOR  DATA  FOR  LIMITED  CHARACTERIZATION  MATERIAL  (PYROLYTIC  GRAPHITES) 


Test 

Number 

Model 

Number 

Test 

Gas 

P 

0 

(atm) 

Ho 

(Btu/lb) 

Tw 

<°R) 

m 

(lb/ft2sec) 

AS 

(mils) 

0 

(sec) 

HITCO  a-t 

PG 

2853-2 

094  S 

1 

2.88 

38,980 

4800 

.000883 

7.525 

98.0 

2857-2 

093  S 

1 

3.20 

44,780 

5100 

.00111 

8.58 

89.0 

2859-4 

099  S 

1 

2.88 

56,260 

5250 

.001503 

7.71 

59.0 

2861-2 

098  S 

: 

3.25 

54,780 

5500 

.001755 

8.85 

58.0 

2865-1 

100  S 

3 

4.73 

11,640 

5000 

.005377 

17.30 

37.0 

2867-2 

096  S 

3 

5.60 

7,660 

4500 

.0006469 

3.38 

60.0 

2872-2 

092  S 

7 

4.23 

12,680 

5000 

.00309 

16.25 

60.5 

2872-4 

097  S 

7 

4.00 

13,710 

4850 

.003069 

15.48 

58.0 

2874-2 

091  S 

7 

4.00 

9,940 

4600 

.00104 

7.35 

81.0 

Supertemp 

PG 

2584-1 

040  S 

1 

3.75 

44,200 

4550 

.000647 

3.35 

59.5 

2586-1 

041  S 

1 

4.28 

65,400 

5050 

.00342 

16.79 

56.5 

2592-3 

043  S 

3 

6.10 

16,500 

4700 

.01292 

40.10 

35.7 

2596-2 

044  S 

4 

5.25 

7,260 

4650 

.01888 

46.80 

28.5 

2597-3 

051  S 

4 

6.45 

8,930 

4750 

.00557 

23.50 

48.5 

2617-5 

072  S 

2 

6.20 

2,630 

4500 

.00477 

18.67 

45.0 

2854-4 

074  S 

1 

3.45 

34,850 

4750 

.0002218 

1.89 

98.0 

2867-1 

073  S 

3 

3.30 

9,560 

4100 

.000104 

0.542 

60.0 

2870-4 

075  S 

7 

4.62 

11,940 

5050 

.00295 

15.50 

60.5 

Pfizer  PG 

2584-4 

045  S 

1 

3.80 

42,400 

4450 

.0001278 

0.667 

60.0 

2586-2 

046  S 

1 

4.30 

64,600 

4450 

.000341 

13.80 

46.5 

2590-1 

047  S 

3 

6.12 

9,130 

3650 

.00636 

16.60 

30.0 

2592-4 

048  S 

3 

5.78 

16,600 

4800 

.01624 

60.00 

42.5 

2596-3 

049  S 

4 

6.45 

6,770 

4550 

.01095 

55.70 

58.5 

2597-2 

050  S 

4 

5.90 

9,242 

5000 

.01126 

48.00 

49.0 

2584-5 

087  S 

1 

3.38 

34,330 

4250 

.000170 

1.45 

98.0 

2860-1 

086  S 

1 

3.75 

54,900 

5600 

.003127 

15.50 

57.0 

54 


TABLE  12.  Concluded 


Model  Test 


Pfizer  _PG  (Concluded) 


5%  SiC  PG 


23%  SiC  PG 

2617-3 

2854-3 

2856-4 

2859-3 

2864-2 

2867-5 

2872-3 

2876-1 

HCF  PG 


2853-3 

2865-2 

2867-7 

2871-1 

2874-3 


Ho 

(Btu/lb) 

T 

w 

(°R) 

m 

(lb/ft2sec) 

12,611 

5350 

.00939 

9,123 

5150 

.00657 

54,730 

5550 

.00287 

11,980 

5000 

.01013 

7,920 

4250 

.00447 

2,540 

3600 

.0000179 

33,270 

4850 

.001057 

45.620 

5250 

.00680 

55,080 

5450 

.010575 

12,400 

5100 

.02231 

7,500 

4550 

.01163 

12,810 

4900 

.00924 

9,850 

4600 

.00503 

38,820 

4850 

.001689 

11,050 

5250 

.02915 

6,960 

4500 

.00117 

13,280 

4850 

.01674 

9,080 

4950 

.00203 

AS 

(mils) 


4.2  BOUNDARY  LAYER  EDGE  SOLUTIONS 


The  free  stream  conditions  at  various  locations  within  the  nozzle  were  calculated  based  on 
the  assumption  of  one-dimensional  isentropic  flow.  This  calculation  can  be  easily  performed  with 
the  isentropic  expansion  option  in  the  GASKET  program.  In  order  to  determine  the  entropy  for  the 
isentropic  expansion  calculations,  the  thermodynamic  state  of  the  plenum  chamber  was  first  cal- 
culated based  on  the  measured  or  evaluated  parameters  (P0>  hQ,  ) . 

4.3  EVALUATION  OF  HEAT  AND  MASS  TRANSFER  COEFFICIENTS 


The  heat  transfer  coefficient  (Peue^)  was  evaluated  directly  from  experimental  measurements 
through  the  following  relationship: 


(p  u C, , ) 

' e e H'cw 


(8) 


Where  qcw  is  the  cold  wall  heat  flux  measured  from  a calibration  run  prior  to  each  ablation  test, 
hQ  is  the  total  enthalpy  which  is  known  from  the  APG  operating  conditions,  and  hcw  is  the  cold  wall 
enthalpy  whicn  is  merely  the  sum  of  the  heat  of  formation  multiplied  by  the  mole  fraction  of  each 
species  in  the  test  gas. 


In  order  to  account  for  hot  wall  effects,  the  Bartz  equation  was  used  as  a scaling  function. 
For  a given  geometry,  the  Bartz  equation  has  the  form: 

0-8  0*2  **  0*  6 


VeVlJ  ~ {prefue>  pref  Pr 


(9) 


Thus,  the  hot  wall  heat  transfer  coefficient  can  be  determined  from 


P 


u Cu 
e e H 


^peueCH^B,hw 

*peueCH*B,cw 


^peueCH^cw 


(10) 


Through  the  Reynold's  analogy,  the  mass  transfer  coefficient  can  also  be  determined.  The  correlation 
frequently  employed  is  (Reference  18). 


« -.  LeV3  , /P T)2/, 
H Le  'ScVef 
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For  the  above  computations,  all  properties  are  determined  from  the  GASKET  program.  The  subscript 
ref  denotes  the  values  at  the  reference  enthalpy  state,  defined  as: 


href  = °’36  he  + 0J9  ho  + °-45  hw 

The  calibration  results  and  the  evaluated  transfer  coefficients  are  shown  in  Tables  13  and  14. 
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4.4  OPEN  SYSTEM  SURFACE  STATE  SOLUTIONS 


ft'.'! 


The  final  step  of  the  data  reduction  procedure  is  to  calculate  the  mole  fraction  of  reactive 
species  at  the  carbon  surface.  By  specifying  the  edge  gas  thermodynamic  state,  Tw,  and  B',  the  sur- 
face state  solution  can  be  determined  using  the  Aerotherm  Chemical  Equilibrium  (ACE)  computer  pro- 
gram. B'  is  the  nondimensional  ablation  variable  and  is  defined  as 


B'  - 


peueCm 


(13) 


The  results  of  this  calculated  data  were  presented  in  Table  9. 


\ 
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SECTION  5 


ROCKET  MOTOR  TEST  RESULTS 

Previous  sections  have  described  how  botn  arc  plasma  generator  and  motor  firing  data  were  in- 
corporated in co  the  kinetics  model.  The  basic  data  used  was  in  the  form  of  mass  removal  rate,  wall 
temperature,  and  wall  partial  pressures.  This  data  is  readily  available  from  arc  plasma  generator 
tests  but  not  from  nozzle  test  firings.  For  a typical  test  firing  of  a full  or  subscale  nozzle, 
only  the  pressure  history  and  final  total  recession  are  reported.  To  deduce  the  appropriate  data, 
an  iterative  procedure  for  correlating  the  data  was  adopted.  This  was  described  in  a five-step 
process  at  the  end  of  Section  3.3.  This  section  will  summarize  the  nozzle  test  data  and  how  it  was 
derived  for  the  five  steps.  Step  1 of  the  correlating  procedure  requires  approximate  values  of  sur- 
face recession  rate,  surface  temperature,  and  species  partial  pressures  at  the  wall.  These  were  ob- 
tained by  performing  a complete  in-depth  conduction  solution  at  the  nozzle  location  of  interest 
using  the  standard  Aerotherm  procedure  (to  be  described  in  Section  7).  Since  the  kinetics  constants 
needed  for  GASKET  are  unknown  at  this  time,  a best  guess  is  used.  Previously  developed  constants 
for  materials  roughly  similar  to  the  material  of  interest  are  typically  used.  Figure  12  illustrates 
a typical  surface  response  prediction.  A complete  history  of  surface  temperature,  surfaie  recession 
rate,  and  cumulative  surface  recession  is  obtained  from  the  CMA  computer  code.  Notice  that  the  data 
needed  for  the  correlation  varies  with  time.  To  deduce  only  one  data  point  per  firing,  it  was  as- 
sumed that  the  conduction  prediction  procedure  accurately  predicted  the  relative  time  variation  of 
surface  temperature  and  recession  rate.  New  kinetics  constants  were  felt  to  improve  the  accuracy 
of  the  total  surface  recession  only.  For  this  reason,  the  point  in  time  used  for  data  in  all  cases 
was  chosen  to  be  the  midpoint  of  the  firing.  Surface  mass  removal  rate  was  deduced  from  the  surface 
recession  rate  as  follows: 

m = S*p  ( 14) 

where  m = mass  removal  rate 

a = surface  recession  rate 
p = material  density 

With  this  data  for  one  or  more  motor  firings,  Step  1 of  the  correlation  procedure  could  be  completed. 
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cal  surface  response  prediction  for  motor  fi 


A second  prediction  was  then  run  with  the  resulting  kinetics  model.  First,  a GASKET  predic- 
tion was  made  (Step  2)  and  then  a new  CMA  prediction  was  run  (Step  3).  This  second  prediction  was 
identical  to  the  first  except  for  the  more  refined  kinetics  constants.  Depending  on  the  accuracy 
of  the  first  guess,  the  total  predicted  recession  may  or  may  not  match  the  measured  recession.  If 
not,  the  newly  predicted  surface  temperature  and  a modified  surface  recession  rate  was  used  to  up- 
date the  data  used  in  the  correlation  (Step  4).  The  modified  surface  recession  rate  was  defined  as: 


S = § 

^Modified  ^predicted*  S 


measured 

predicted 


(15) 


With  this  information  the  final  kinetics  correlation  model  was  determined  (Step  5).  The  final  motor 
firing  data  used  for  all  of  the  kinetics  models  are  summarized  in  Table  15. 

A final  complete  prediction  was  made  for  each  motor  firing  to  verify  the  kinetics  model.  In 
most  cases  the  model  was  verified  by  predicting  the  measured  total  surface  recession  within  25  per- 
cent. These  verification  predictions  are  referred  to  as  correlation  studies  and  are  presented 
briefly  in  Section  7 and  in  full  detail  in  a second  final  report  (AFRPL-TR-76-71 ). 
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SECTION  6 


EVALUATION  OF  KINETIC  CONSUMPTION  RATES 

Two  procedures  were  applied  to  evaluate  kinetic  consumption  rates  of  carbon  materials.  The 
first  procedure,  namely  full  characterization  studies,  correlated  the  carbon  consumption  rate  data 
from  both  APG  and  motor  firing  testings  as  a function  of  the  following  parameters: 

• Surface  temperature 

• Boundary  layer  edge  pressure 

• The  chemical  composition  of  the  propellant  gas 

The  correlation  function  used  was  based  on  the  Langmuir-Hinshelwood  model.  This  model  takes  into 
account  such  detailed  mechanisms  as  chemadsorption , desorption,  and  heterogeneous  reactions.  Thus, 
not  only  is  the  obvious  ablation  performance  evaluated,  but  significant  physical  insights  may  be 
obtained  on  the  competition  of  active  sites  by  reactive  species  and  poisoning  effects  by  halogen 
species  This  correlation  procedure  was  described  in  Sections  3 and  4. 

The  second  procedure,  namely  limited  characterization  studies,  used  the  APG  to  evaluate  the 
relative  ablation  performance  of  the  screening  carbon  materials.  This  procedure  included  the  fol- 
lowing items: 

• Selection  of  test  gases  that  characterized  the  actual  motor  firing  surface  kinetics 

• Extrapolation  of  the  APG  results  to  determine  the  relative  ablation  performance  at 
actual  rnotor-firing  conditions  (T  = 5500°R  - 6000°R) 

Fully  characterized  materials  were  used  as  a baseline  material  for  the  screening  materials  fur  the 
same  generic  class  in  order  to  determine  their  relative  ablation  performance. 

6.1  RESULTS  OF  FULL  CHARACTERIZATION  STUDIES 

Kinetic  constants  and  inhibition  coefficients  of  the  following  materials  were  determined  from 
the  APG  and  motor  firing  data  using  an  Aerotherm  least  squares  data  optimization  program: 

• Supertemp  Edge  P^ 

• 15*  SiC/PG 
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ATJ  Bulk  Graphite 


• G-90  Bulk  Graphite 

• Pyrocarb  901  Carbon/Carbon 

Carbitex  700  data  was  available  from  previous  work  and  its  correlation  functions  were  reassessed  in 
this  study.  The  quality  of  the  data  optimization  are  shown  in  Figures  13  through  17  and  numerical 
results  are  presented  in  Table  16. 

The  maximum  standard  deviation  obtained  in  these  correlations  is  o = 0.16.  From  a statistical 
point  of  view,  these  correlations  are  considered  to  be  satisfactory.  The  accuracy  of  these  correla- 
tions, however,  must  be  further  tested  by  comparing  predictions  with  the  motor  firing  data.  The 
results  of  these  comparisons  will  be  discussed  in  Section  7. 

The  following  conclusions  were  reached  from  the  full  characterization  studies: 

• The  - C*  reaction  has  a high  activation  energy  and  a relatively  large  number  of  active 
sites.  This  kinetic  reaction  becomes  significant  at  surface  temperatures  above  5500°R  and 
is  a major  contributor  to  carbon  surface  recession  at  motor  firing  conditions. 

• The  H20  and  COg  - C*  reactions  have  relatively  low  activation  energies  and  relatively 
small  numbers  of  active  sites  compared  with  the  - C*  reaction.  These  reactions  play 
a dominant  role  in  carbon  surface  recession  at  temperatures  between  4000°R  and  5500°R. 
Their  effects,  however,  diminish  due  to  the  dissociation  cf  H^O  and  C02  as  the  surface 
temperature  increases. 

• HCJl  retards  the  CO^  and  H^O  - C*  reactions,  but  has  no  effect  on  the  - C*  reaction. 

• G-90  graphite  is  slightly  more  active  with  respect  to  reactants  of  H-,0  and  CO^  than  the 
other  three  full  characterization  materials, 

i Carbon  sublimation  results  in  a nonnegligible  rate  of  carbon  consumption  at  anticipated 
rocket  motor  temperatures.  The  significance  of  sublimation  becomes  increasingly  im- 
portant as  propellant  and  surface  temperature  increase,  especially  for  materials  with 
low  carbon  ablation  rates. 

Because  of  the  complexity  of  the  kinetic  formulations,  it  is  not  possible  to  make  general 
comments  about  the  magnitudes  of  each  reaction  for  a given  material.  However,  one  can  compare  the 
reaction  rate  coefficients,  as  presented  in  Table  16,  as  a function  of  temperature.  These  compari- 
sons are  shown  in  Figures  18  through  22  for  the  materials  characterized  in  this  study.  In  addition, 
results  from  earlier  studies  for  Carbitex  700  and  Atlantic  Research  Corporation  layer  PG  are  shown 
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Figure  13.  Supertemp  PG,  c plane  surface  kinetics. 
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Figure  16.  Results  of  G-90  kinetic  correlation. 
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TABLE  16.  CARBON  KINETICS  COEFFICIENTS  FOR  FULL  CHARACTERIZATION  MATERIALS 
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in  Figures  23  and  24.  For  the  Carbltex  700  material,  data  from  Reference  19  (Thiokol)  were  re- 
correlated to  be  consistent  with  first-order  kinetics.  For  layer  PG,  the  results  were  taken  directly 
from  Reference  6 (Schaefer). 


Before  the  reader  draws  any  conclusions  from  these  figures,  he  is  cautioned  to  qualitatively 
account  for  the  facts  that: 

1.  Figures  18  through  24  do  not  include  the  effects  of  inhibitors  and  poisons 

2.  The  actual  carbon  consumption  rate  is  the  product  of  the  reaction  rate  coefficient 
(modified  by  inhibitor  terms)  and  the  local  partial  pressure  of  the  reactant  at  the 
reacting  surface 

Because  of  these  two  constraints,  valid  comparisons  of  the  reaction  rates  can  only  be  made  for 
specific  propellants.  Conclusions  based  only  on  Figures  18  through  24  can  be  very  misleading.  For 
instance,  the  edge  PG  results  of  Figure  18  indicate  that  the  reaction  rate  for  H.,0  at  6000°R 
(1/TW  = 3 x 1 0~ 4 K~ 1 ) is  about  two  orders  of  magnitude  greater  than  that  for  H^.  However,  GASKET 
code  calculations  show  that  the  mole  fraction  of  HgO  is  about  four  orders  of  magnitude  less  than 
that  of  H2  (Figure  8).  Hence  the  carbon  consumption  rate  for  the  H.,0  - C*  reaction  will  about  two 
orders  of  magnitude  less  than  that  for  the  H2  - C*  reaction. 

6.2  RESULTS  OF  LIMITED  CHARACTERIZATION  STUDIES 
H2  Test  Gas 

From  the  results  of  the  full  characterization  studies,  it  was  concluded  that  H2  is  the  major 
contributor  for  carbon  consumption  in  the  temperature  range  corresponding  to  typical  motor  firing 
conditions.  H20  and  C02  are  somewhat  less  reactive  in  this  same  temperature  range.  In  fact,  their 
significance  diminishes  due  to  the  disappearance  of  both  species  by  dissociation.  These  findings 
suggest  that  H2  is  the  most  appropriate  test  gas  for  material  screening  in  the  APG. 

An  H2  test  gas  simplifies  the  interpretation  of  data.  The  reason  Is  that  the  H2  - C*  reac- 
tion is  not  strongly  inhibited  by  other  gas  species,  halogens  included.  Thus,  the  kinetic  rate  of 
this  reaction  closely  obeys  the  Arrhenius  expression,  i.e.,: 

-E/RTw 

mc  = A e P^ 

By  plotting  log  (m  /P  ) versus  1/T  with  the  APG  data,  a straight  line  should  be  obtained.  This 
straight  line  can  tin.ii  be  extrapolated  to  the  motor  firing  temperature  range  to  estimate  the  ab- 
lation performance  of  the  screening  materials.  Some  engineering  judgment  must  then  be  used  to  in- 


clude the  effect  of  H20  and  C02  reaction. 


>/T„  X 10-  (°K- 1 ) 

Figure  24.  Reaction  rate  coefficients  for  Atlantic 

Research  Corporation  layer  pyrolytic  graphite 


Figures  25  and  26  show  the  relative  ablation  performance  of  the  following  carbon/carbon 
materials  in  the  APG. 

• HRX  5125 

• HRX  5875 

• MDAC  3-D  C/C 

• Pyrocarb  903 

• Pyrocarb  903  HD 

These  are  compared  to  the  full  characterization  material,  Pyrocarb  901.  Based  on  these  results, 
the  ranking  of  ablation  performance  at  motor  firing  conditions  was  determined  and  is  presented  in 
Table  17. 


TABLE  17.  RANKING  OF  CARBON/CARBONS  BASED 
ON  MASS  CONSUMPTION  PERFORMANCE 


Ranking 

Material 

Designation 

Source 

Measured 

Density 

(gm/cc) 

Reinforcement  Precursor 

1 

3-D  C/C 

MDAC 

2.0 

Rayon  yarn  3D  orthogonal  weave 

2 

HRX  5125 

Haveg 

1.55 

Rayon  fabric  0°-22°-45°  ply  orientation 

3 

Pyrocarb  903  HD 

Hitco 

1.90 

PAN  fabric  0°-45o-90o  ply  orientation 

3 

Pyrocarb  903 

Hi  tco 

1.84 

Same  as  Pyrocarb  903  HD 

4 

HRX  5875 

Haveg 

1.82 

PAN  fabric  0°-22o-45°  ply  orientation 

5 

Pyrocarb  901 

Hitco 

1.65 

Rayon  fabric  0o-45°-90<>  ply  orientation 

From  the  carbon/carbon  data,  an  attempt  was  made  to  relate  the  measured  relative  performance 
to  the  construction  of  the  composite.  In  particular,  the  relative  orientation  of  the  adjacent  plies, 
the  reinforcement  precursor  and  the  composite  density  were  considered.  No  ply  orientation  relation- 
ship was  observed;  however,  some  general  relat.onships  were  found  for  the  reinforcement  precursor  and 
the  comoos ite  density. 

Generally,  it  was  found  that: 

• Rayon  precursor  composites  are  superior  to  PAN  precursor  composites  of  equal  density 

t Within  classes  of  materials  that  have  the  same  precursors,  high  density  materials  per- 
form better  than  low  density  materials 


83 


m/P  ( Ibm/sec-f t2atm) 


A.-ISI87 


nj/P  ( lbro/sec-ft2atm) 


Figure  26.  Carbon/carbons,  APG  test  gas 


These  conclusions  were  made  based  on  the  results  presented  in  Table  17  and  Figure  27.  Figure  27 
was  constructed  from  the  results  of  Figures  25  and  26  by  extrapolation  of  a mean  line  through  the 
APG  data.  For  reference,  these  lines  are  shown  in  Figure  28. 

Figure  29  shows  the  performance  variations  of  pyrolytic  graphite  materials  in  a hydrogen 
environment.  These  materials  were  fabricated  by: 

t Hitco  and  Atlantic  Research  Corp.  (ARC)  for  a-b  plane 

• Supertemp  and  Pfizer  for  c plane 

The  ARC  and  Supertemp  materials  were  used  as  references  and  are  shown  in  Figure  29  as  lines.  It  is 
apparent  from  this  figure  that  the  material  source  Is  not  Important  for  the  ablation  performance  of 
pyrolytic  graphites. 

The  relative  performance  of  modified  pyrolytic  graphite  materials  in  a hydrogen  environment 
are  shown  in  Figure  30.  These  materials  included: 

• 5,  15  and  23  percent  (by  weight)  silicon  carbide  codeposited  pyrolytic  graphite 

• 65  percent  (by  weight)  hafnium  carbide  codeposited  pyrolytic  graphite 

Although  it  is  presumptuous  to  reach  any  conclusions  from  the  limited  amount  of  APG  data,  it  appears 
that: 

• 23  percent  SiC/PG  performs  similarly  to  edge  oriented  PG 

• 5 percent  SiC/PG  performs  marginally  better  than  15  percent  SiC/PG,  which  performs 
noticeably  better  than  23  percent  SiC/PG.  It  should  be  noted  that  the  23  percent  SiC/PG 
was  highly  nonuniform  so  the  23  percent  is  merely  a nominal  indication  of  the  SiC  content. 

• 65  percent  HFC/PG  material  performs  similarly  to  edge-oriented  PG;  however,  when  the  mass 
loss  is  transformed  to  surface  recession,  the  performance  of  this  material  is  quite  good 
since  it  has  a density  of  about  9 gm/cc  compared  to  2.2  gm/cc  for  PG 

Figure  31  shows  the  relative  ablation  performance  of  P03  and  ATJ-S  graphite  relative  to  the 
full  characterization  material,  G-90  qraphlte.  The  ablation  performance  of  these  three  materials  is 
in  the  order  ATJ-S,  P03,  and  G-90,  although  the  differences  between  these  materials  is  not  signifi- 
cant. 

H.,  - O2  Test  Gas 

APG  results  for  Hg  - Og  gas  mixtures  are  presented  without  further  discussions  in  Figures 


32  through  35. 
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Figure  33.  Ablation  data,  pyrolytic  graphites 
APG  test  gas  H9/O0. 
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Figure  34.  Ablation  data,  modified  pyrolytic 
graphites,  APG  test  gas  H2/02> 
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Figure  35,  Ablation  data,  buU  graphites,  APG  test  gas  - 
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SECTION  7 


PREDICTED  ROCKET  NOZZLE  PERFORMANCE 

The  GASKET2  code  (Reference  9),  which  includes  the  kinetic  models  presented  in  Table  16,  was 
used  to  predict  many  actual  motor  firings.  Specifically,  15  correlation  studies  and  6 performance 
studies  were  performed.  The  purpose  of  the  correlation  studies  was  to  provide  basic  motor  firing 
data  to  be  incorporated  in  the  kinetics  correlation.  Once  this  data  was  provided,  the  correlation 
studies  were  used  to  verify  the  resulting  model.  These  verification  calculations,  along  with  per- 
formance studies  in  support  of  AFRPL  technology  programs,  have  been  performed  and  will  be  summarized 
in  this  section. 

7.1  ANALYSIS  PROCEDURE 

The  analysis  procedure  used  to  predict  motor  firings  is  described  in  this  section.  The  events 
which  occur  near  the  internal  surface  of  an  ablative  rocket  nozzle  are  illustrated  in  Figure  36.  An 
inviscid  flow  field  comprises  the  bulk  of  the  flow  of  propellant  gases  through  the  nozzle.  Typically, 
the  propellant  contains  less  than  21  percent  (by  weight)  of  aluminum,  which  forms  liquid  alumina 
particles  that  flow  with  the  other  products  of  combustion  through  the  nozzle.  Near  the  surface, 
the  flow  field  is  represented  by  a boundary  layer  in  which  alumina  particles  are  not  considered  to 
be  present.  Chemically  reactive  species  diffuse  through  this  boundary  layer  and  cause  surface  ab- 
lation. The  nozzle  thermal  protection  material  must  respond  to  three  sources  of  energy  transfer 
the  convection  and  diffusion  of  energy  across  the  boundary  layer  and  the  radiation  or  energy  from 
high  temperature  alumina  particles  in  the  inviscid  case  flow.  All  of  these  energy  events  result  in 
a given  amount  of  energy  being  conducted  into  the  material  to  cause  internal  component  heating. 

The  procedure  used  for  ablation  predictions  treats  the  inviscid  flow  field,  the  boundary 
layer  and  the  radiant  energy  transfer  separately.  Although  the  analyses  are  separate,  the  final 
solution  is  properly  coupled  as  illustrated  in  Figure  37.  Each  area  will  be  discussed  individually 
in  the  following  paragraphs. 

Flow  Field 

The  flow  field  was  analyzed  with  the  "Aerotherm  Chemical  Equilibrium  (ACE)"  computer  code 
(Reference  20).  This  is  a chemical  equilibrium  code  which  is  used  to  compute  the  local  thermodynamic 
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Figure  37.  Thermal  analysis  flow  chart. 


state  throughout  the  nozzle.  The  primary  variables  of  concern  are  the  local  pressures,  temperatures 
and  chemical  species.  These  are  determined  by  performing  many  isentropic  closed-system  equilibrium 
solutions  throughout  the  nozzle.  The  ACE  code  can  handle  any  arbitrary  real  gas  by  knowi-ng  the 
chamber  state  (temperature  and  pressure)  and  the  propellant  elemental  composition.  Both  gaseous 
and  liquid  phases  are  computed,  but  no  thermal  or  velocity  lags  are  accounted  for.  The  closed  sys- 
tem solutions  are  associated  with  given  positions  in  the  nozzle  by  assuming  one-dimensional  flow 
and  using  conservation  of  mass  relations. 

Boundary  Layer 

The  analysis  of  the  boundary  layer  Is  performed  by  using  the  "Aerotherm  Real  Gas  Energy 
Integral  Boundary  Layer  (ARGEIBL)"  computer  code  (Reference  21).  This  is  an  energy  integral  tech- 
nique which  can  handle  any  arbitrary  real  gas  with  the  input  of  general  Mollier-type  tables  gener- 
ated by  the  ACE  computer  code.  Since  no  alumina  particles  are  in  the  boundary  layer,  these  proper- 
ties are  for  the  gas  phase  only.  Other  input  consists  of  edge  state  variables  (P  and  T ) generated 
by  the  flow  field  analysis.  Arbitrary  wall  temperatures  are  also  handled  by  this  code.  Since  these 
are  not  predicted  until  a conduction  analysis  Is  performed,  an  estimation  of  the  wall  temperature  is 
required.  If  this  estimation  proves  to  be  far  in  error,  an  iteration  back  through  the  boundary  layer 
analysis  is  necessary. 

Turbulent  flow  is  assumed  and  the  boundary  layer  is  started  (with  zero  energy  thickness)  at 
the  nose  of  submerged  nozzles  and  at  the  base  of  the  aft  closure  for  conventional  nozzles.  The 
nos.,  of  a submerged  nozzle  is  defined  as  that  point  furthest  axially  upstream  of  the  throat  and  the 
base  of  an  aft  closure  is  considered  as  the  point  where  the  nozzle  Insulation  mates  to  the  insula- 
tion in  the  propellant  case. 

To  compute  the  heat  transfer  coefficient  for  ablating  surfaces,  a blowing  correction  is  made 
to  the  ARGEIBL  results.  This  is  a multiplying  factor  which  is  a function  of  the  mass  injected  into 
the  boundary  layer.  The  correction  is  made  internal  to  the  CMA  code  and  a detailed  discussion  can 
be  found  in  Reference  22. 

Previous  experience  has  shown  that  the  heat  transfer  coefficient  predicted  by  the  ARGEIBL 
code  is  high  and  should  be  multiplied  by  0.75;  that  is: 

peueCH  = 0 75  (peueCH*ARGEIBL  (16) 

where  (peue^ARGEIBL  '5  the  convective  heat  transfei  coefficient  predicted  by  the  ARGEIBL  code. 

This  factor  has  been  included  in  all  of  the  predictions  reported.  For  the  transient  in-depth  thermal 
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analysis  presented  in  subsequent  sections.  Equation  (16)  is  modified  further  to  account  for  the  ac- 
tual chamber  pressure  history  by  using  the  relationship: 


peueCH  ’ °-75  (p-2-*’'1 


e e H'ARGEIBL 


ave 


07) 


where  PQ  is  the  instantaneous  measured  (or  predicted)  chamber  pressure  and  PQ 
chamber  pressure  assumed  for  the  flaw  field  analysis. 


ave 


is 


the  average 


Surface  Thermochemistry 

The  surface  thermochemistry  analysis  is  performed  by  using  the  Aerotherm  "Graphite  Surface 
Kinetics  (GASKET2)"  computer  code  (Reference  9).  It  computes  the  surface  state  (complete  wall  gas 
thermodynamic  and  chemical  make-up)  of  many  graphitic  materials  exposed  to  a corrosive  rocket  nozzle 
environment.  The  reactions  which  occur  between  the  graphite  material  and  the  propellant  gases  are 
considered  kinetically-controlled  and  are  modeled  by  Arrhenius-Langmuir  type  reaction  rate  equations. 
This  code  can  handle  arbitrary  propellant  gases  and  arbitrary  surface  reaction  kinetic  constants. 
Kinetics  constancs  used  for  the  predictions  discussed  in  this  section  are  currently  included  in  the 
GASKET2  code  for  the  following  graphitic  materials: 

• a-b  plane  PG 

• c plane  PG 

• 15%  SiC/PG 

• Pyrocarb  901  carbon/carbon 

• Carbitex  700  carbon/carbon 

• ATJ  bulk  graphite 

• G-90  bulk  graphite 


As  input,  the  GASKET2  code  requires  the  local  thermodynamic  edoe  state  (P„  and  T ) and  the 

6 6 

elemental  composition  of  the  gas,  which  are  obtained  directly  from  the  ACE  generated  flow  field 
analysis.  The  mass  transfer  coefficient  is  also  required  and  is  obtained  directly  from  the  heat 
transfer  coefficient  computed  by  ARGEIBl  from  the  following  relationship: 


" ©,/i 


peueCH 


(18) 


where 

Pr  = Prandtl  number 
Sc  = Schmidt  number 
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The  Prandtl  and  Schmidt  numbers  are  both  obtained  from  the  ACE  generated  Mol  Tier  input  to  the 
ARGEIBL  code.  For  nozzles  whose  pressure  traces  are  far  from  constant  throughout  a firing,  Equation 
(17)  shows  that  the  heat  transfer  coefficient  (and  thus  mass  transfer  coefficient)  will  vary  with 
time.  For  these  cases,  up  to  three  surface  thermochemistry  analyses  are  performed  for  a given  firing. 
The  conduction  solution  procedure  then  interpolates  between  these  solutions  to  obtain  instantaneous 
firing  and  boundary  conditions. 

To  model  the  radiative  boundary  conditions  during  nozzle  firings,  a parallel  plate  model  was 
used.  This  model  applies  to  aluminized  propellants  and  assumes  that  the  particle  laden  stream  of 
combustion  products  Is  optically  thick  and  that  it  exchanges  radiant  energy  with  the  surface  as  if 
the  stream  and  wall  were  parallel  plates.  In  this  way,  multiple  reflections  between  the  wall  and 
stream  were  taken  into  account.  In  addition,  the  assumption  was  made  that  both  the  stream  and  wall 
behave  as  gray  bodies  and  that  they  emit  and  reflect  radiant  energy  diffusely.  Based  on  the  above 
assumptions,  the  net  radiant  heat  flux  relation  is  given  as 


%et  rad  * eeff^s  ' cTw^ 

where 

ceff  " Effective  emlssivity  = 


(19) 


ew  -Wall  material  emissivity 
es  - Particle  laden  stream  emissivity 
o - Stefan-Bol tzmann  constant 

T$  - Free  stream  (edge)  temperature 

Tw  - Wall  temperature 

To  determine  tne  effective  emissivity  using  Equation  (19),  the  stream  emissivity  was  defined  as 

es  « 1 - exp  (-C^  pD)  (20) 

where 

C - Empirical  constant  (0.808) 

n - Percentage  of  aluminum  loading 

p -local  density  of  propellant  combustion  species  (lb/ft3) 

D -Local  beam  length,  usually  taken  as  the  diameter  (in.) 
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In-Depth  Conduction 


Both  the  Aerotherm  "Axi symmetric  Transient  Heating  and  Material  Ablation  (ASTHMA)"  (Reference 
23)  and  the  Aerotherm  "Charring  Material  Ablation  (CMA)"  (Reference  22)  computer  codes  were  used  for 
the  in-depth  conduction  analyses.  Bascially,  the  ASTHMA  code  is  two-dimensional  and  CMA  is  one- 
dimensional. Both  have  the  capability  of  handling  temperature  dependent  material  properties.  CMA 
also  has  the  capability  of  modeling  materials  which  internally  decompose.  This  capability  was  not 
exercised  in  the  motor  firing  predictions  since  only  graphitic  surface  materials  were  considered. 

Figure  37  shows  that  the  flow  field  (H  ),  boundary  layer  (pQuQCu),  surface  thermochemistry 
(B1  map),  and  radiation  (net  radiation  flux)  analyses  are  all  used  as  input  to  CMA  and  ASTHMA.  Also 
shown  as  input  are  the  component  geometry  and  material  properties.  Material  properties  are  re- 
ferenced for  each  nozzle  firing  analyzed,  but  a majority  of  the  properties  were  obtained  from  the 
"Aerotherm  Graphitic  Material  Handbook  of  Thermophysical  Properties"  (Appendix  B). 

7.2  CORRELATION  STUDIES 

This  section  summarizes  the  15  correlation  studies  which  were  performed.  The  following  five 
kinetics  models  were  verified  by  these  studies: 

• c plane  PG 

t 15%  SiC/PG 

• G-90  bulk  graphite 

• ATJ  bulk  graphite 

t Pyrocarb  901  carbon/carbon 

Predictions  were  made  for  test  nozzles  fired  by  several  companies,  using  various  propellants. 
The  firings  were  both  fullscale  and  subscale,  with  both  submerged  and  conventional  nozzles.  Those 
chosen  in  most  cases  were  typical  of  advanced  ICBM  conditions. 

The  motors  used  for  the  15  correlation  studies  are  as  follows: 

C Plane  PG 

1.  TCC,  MMIII,  HTPB  Demo 

2.  CSD,  C-4,  3rd  Stage  Demo 

3.  ASPC,  C-4,  1st  Stage  Demo 

4.  TCC,  C-4,  1st  Stage  Demo  (FST-002) 
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5.  ASPC,  Nimrod  14,  Subscale 


6.  ASPC,  Nimrod  15,  Subscale 

7.  CSD,  MX,  Lower  Stage,  Subscale 

8.  TCC,  C-4,  30-Inch  Material  Evaluation  (3SF-24) 

15%  S1C/PG 

9.  ARC,  7-inch 

10.  ARC,  3.5-inch 

11.  ASPC,  Nimrod  6,  Subscale 

G-90 

12.  Rocketdyne  Condor 

13.  CSD,  FW-5 

ATJ 

14.  Hercules,  X259  Antares  II 
Pyrocarb  90? 

15.  TCC,  C-4,  30-inch  Material  Evaluation  (3SF-24) 

Complete  details  of  all  of  the  predictions  are  given  in  a second  final  report  ( AFRPL-TR- 76- 7 1 ). 
Since  much  of  the  work  is  classified,  only  a brief  summary  of  the  results  will  be  presented  here. 

Figure  38  shows  that  12  of  the  studies  predicted  total  surface  recession  within  25  percent. 

The  three  that  fell  beyond  the  ±25  percent  band  were  studies  1,  5,  and  6.  The  second  final  report 
discusses  in  detail  the  comparisons  between  the  predictions  and  the  data. 

7.3  PERFORMANCE  STUDIES 

As  mentioned  previously,  the  performance  studies  were  done  in  support  of  AFRPL  technology 
programs.  A total  of  six  were  performed  and  are  summarized  as  follows: 

1.  Hercules  3rd  Stage  MX  (Carbitex  700}  (Reference  10) 

2.  Rocketdyne  Condor  (Reference  11) 

3.  Hercules  3rd  Stage  MX  (Pyrocarb  901)  (Reference  12) 

<4.  Standard  7-inch  Nozzle  (Throat  Sensitivity)  (Reference  13) 
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larison  of  measured  and  predicted  ablation 
correlation  studies. 


105 


5.  Standard  7-inch  Nozzle  (Nosecap  Sensitivity)  (Reference  14) 

6.  BATES  Motor  (Thermostructural ) (Reference  15) 

Each  study  was  reported  in  detail  as  individual  nozzle  bulletins  published  under  this  contract.  In 
the  following  paragraphs  a brief  overview  is  given  of  each  study. 

Study  1 - Hercules  3rd  Stage  MX  (Carbitex  700)  (Reference  10) 

The  purpose  of  this  study  was  to  compare  the  predicted  performance  of  a Hercules  3rd  Stage 
MX  test  nozzle  fired  in  an  HTPB  propellant  environment  (AFRPL,  January  1976)  versus  the  original 
design  environment  of  an  18.5  percent  A1  XLOB  propellant.  The  results  of  the  Aerotherm  performance 
study  were  used  to  select  a firing  time  for  the  nozzle  test. 

A complete  two-dimensional  thermal  analysis,  including  surface  recession,  was  performed  for 
two  carbon/carbon  components  upstream  and  downstream  of  the  throat  and  five  pyrolytic  graphite 
washers  which  formed  the  throat  pack.  Another  complete  thermal  analysis  was  performed  assuming  a 
21  percent  At,  90  percent  solids  HTPB  propellant.  The  chemical  compositions  and  ideal  chamber  con- 
ditions for  the  propellants  are  listed  in  Table  18.  The  chamber  pressure  history  used  is  presented 
in  Figure  39.  The  nozzle  geometry,  obtained  from  Hercules,  is  presented  in  Figure  40. 

Typical  Internal  temperatures  predicted  by  the  ASTHMA  code  are  shown  In  Figure  41.  This  par- 
ticular prediction  is  for  the  HTPB  propellant  after  60  seconds.  Conclusions  reached  from  this  study 
were  that  the  HTPB  propellant  was  more  corrosive  than  the  XLDB  propellant  and  that  the  nozzle  was  ex 
pected  to  withstand  a 60-second  firing  in  the  HTPB  environment.  Figure  41  shows  that  approximately 
50  percent  of  the  leading  PG  washer  has  receded  by  60  seconds.  Although  it  was  felt  that  this  nozzl 
could  sustain  a 60-second  firing,  it  was  recommended  that  it  be  fired  in  the  HTPB  propellant  for  40 
seconds. 

This  test  occurred  in  March  1976  and  was  successful.  A detailed  comparison  of  measured  and 
predicted  recession  has  not  been  performed  to  date. 

Study  2 - Rocketdyne  Condor  (Reference  11) 

For  this  performance  study,  a one-dimensional  ablation  analysis  of  the  threat  was  used  to 
compare  the  accuracy  of  the  newly  developed  edge  PG  kinetics  model  to  Condor  firing  data.  This 
motor  was  chosen  to  demonstrate  the  applicability  of  the  kinetics  modeling  to  conditions  other  than 
that  of  typical  ICBM  type  motors.  The  Condor  firing  durations  were  on  the  order  of  200  seconds  in 
a low  aluminum,  low  flame  temperature  propellant.  The  nozzle  is  a blast  tube  type  with  a throat 
diameter  of  1,1  inches.  The  ICBM  type  test  motors  (which  were  primarily  used  for  developing  the 
current  kinetics  model)  typically  have  firing  durations  of  less  than  60  seconds  in  high  performance 
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TABLE  18.  PROPELLANT  DATA 


Designation  - XLDB  (18.5*  A1 ) 

Formulation  - Element 

Gm  Atoms/100  Gms 

Hydrogen 

2.3475 

Carbon 

1.2745 

Nitrogen 

1.8271 

Oxygen 

2.3897 

Aluminum 

0.6857 

Ideal  Chamber  Conditions  - Pressure 

= 800  psia 

Temperature  = 6700°R 

Designation  - HTPB  (21%  A1 , 90%  solids) 

Formulation  - Element 

Gm  Atoms/ 100  Gms 

Hydrogen 

3.4864 

Carbon 

0.6933 

Nitrogen 

0.5933 

Oxygen 

2.3766 

Aluminum 

0.7784 

Chlorine 

0.5873 

Ideal  Chamber  Conditions  - Pressure 

= 800  psia 

Temperature  = 6649°R 

Pressure 


*18* 


j 


Figure  40.  Nozzle  geometry,  Hercules  test  nozzle. 
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Figure  41.  Predicted  nozzle  response  to  HTPB  propellant,  60.0  seconds. 


propellants  (18  - 21  percent  aluminum  and  6000°R  - 6900°R  flame  temperatures).  These  motors  are 
usually  submerged  and  have  throat  diameters  from  4 to  15  inches. 

The  firing  data  used  was  obtained  from  Rocketdyne  Solid  Rocket  Division  in  MacGregor,  Texas. 

A series  of  seven  firings  were  given,  each  of  which  differed  mainly  in  the  prefire  temperature  of 
the  motor.  The  resulting  throat  erosion  rate  was  more  than  an  order  of  magnitude  less  than  that 
exhibited  by  typical  ICBM  type  motors  and  varied  considerably  among  the  seven  firings.  Both  the 
propellant  and  firing  data  are  confidential  and  are  not  included  in  this  report.  The  geometry  is 
given  in  Figure  42. 

Figure  43  shows  the  predicted  surface  recession  as  a function  of  time.  The  predictions  were 
found  to  be  approximately  80  percent  higher  than  the  average  measured  recession.  This  is  an  accept- 
able prediction  for  the  following  reasons: 

• Firing  data  indicates  that  alumina  condensed  out  on  PG  washers 

o Data  scatter  was  large 

• Measured  recession  was  very  small  (making  a percentage  comparison  invalid  as  an  indica- 
tion of  accuracy) 

Study  3 - Hercules  3rd  Stage  MX  (Pyrocarb  901)  (Reference  12) 

This  performance  study  is  almost  identical  to  Study  1.  Recall  that  Study  1 incorporated  PG 
washers  in  the  throat  and  carbon/carbon  components  up  and  downstream  of  the  washers.  The  carbon/ 
carbon  material  has  a specific  gravity  of  1.60  and  was  made  by  Kaiser.  At  the  time  of  the  predic- 
tion, the  only  carbon/carbon  kinetics  model  available  was  one  developed  by  Aerotherm  for  Thiokol 
Corporation  (P.0.  414011).  This  kinetic  node!  was  specifically  developed  for  Carbitex  700  (S.G.  = 
1.5)  and  was  used  for  this  prediction  by  adjusting  the  material  to  a 1.6  specific  gravity. 

After  Study  1 was  performed  a kinetics  model  for  Pyrocarb  901  (S.G.  = 1.83)  was  developed 
under  this  contract.  Preliminary  indications  showed  that  the  kinetics  model  was  markedly  different 
from  that  developed  for  Carbitex  700.  With  the  indication  that  some  carbon/carbon  materials  behave 
differently  from  others,  it  was  of  interest  to  rerun  the  prediction  of  Study  1 with  the  newly  devel- 
oped Pyrocarb  901  model  (S.G.  = 1.83). 

As  in  performance  Study  1,  two  complete  ASTHMA  analyses  were  performed  for  the  carbon/carbon 
components  and  PG  washers  in  the  region  of  the  throat.  One  was  for  an  18.5  percent  AX.  XLDB  propel- 
lant and  the  other  was  for  a 21  percent  A&,  90  percent  solids  HTPB  propellant. 
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Figure  44  shows  the  predicted  isothermal  profiles  after  60  seconds  of  exposure  to  the  HTPB 
environment.  As  expected  from  the  kinetics  data,  the  predictions  indicated  more  recession  for  the 
Pyrocarb  901  than  the  Carbitex  700.  This  can  be  seen  by  a comparison  between  Figures  41  and  44. 

This  study  showed  that  the  Pyrocarb  901  could  not  sustain  a 60-second  firing  without  severe  under- 
cutting of  the  PG  washers. 

Table  19  is  included  as  a further  comparison  of  Studies  1 and  3.  Basically,  this  shows  that 
the  XL0B  propellant  is  not  as  severe  as  the  HTPB  and  that  predicted  recession  for  Pyrocarb  901  is 
more  than  Carbitex  700  for  comparable  conditions. 

Study  4 - Standard  7-Inch  Nozzle  (Throat  Sensitivity}  (Reference  13) 

The  purpose  of  this  study  was  to  compare  the  predicted  performance  of  four  possible  throat 
materials  in  the  standard  7-1nch  test  nozzle  configuration,  namely: 

• Pyrocarb  901  carbon/carbon  (density  =1.83  g/cc} 

• c plane  pyrolytic  graphite  (edge) 

• a-b  plane  pyrolytic  graphite  (layer) 

i 15*  SiC/PG 

For  each  material,  three  propellant  environments  were  considered. 

• HTPB 

• XLDB 

• PEG/FEF0 

For  each  propellant/material  combination,  a complete  one-dimensional  thermal  analysis  including  sur- 
face recession  was  carried  out  at  the  throat.  A total  of  12  nozzle  analyses  were  performed. 

The  chemical  composition  and  actual  chamber  conditions  for  the  propellants  considered  are 
summarized  in  Table  20.  The  firing  duration  was  60  seconds  at  a cnamber  pressure  of  1000  psia.  The 
nozzle  geometry  is  shown  in  Figure  45(a).  Figures  45(b)  and  45(c)  show  the  geometries  of  the  throat 
insert  depending  on  whether  a coated  or  solid  component  is  used. 

Of  most  Interest  in  this  study  is  a comparison  of  the  predicted  surface  recession  for  all  of 
the  material/propellant  combinations.  The  predicted  total  recessions  and  average  recession  rates 
are  tabulated  in  Table  21.  The  average  recession  rates  are  illustrated  in  Figure  46.  This  shows 
that  the  Pyrocarb  901  recedes  at  least  a factor  of  two  or  more  than  the  other  materials  in  a given 
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60  seconds. 


TABLE  19.  RECESSION  RATE  SUMMARY,  HERCULES  3RD  STAGE  MX  NOZZLE 


Location 

Propellant 

Material 

Average  Recession 
Rate  O 60  secs 
(mils/sec) 

Upstream  Throat 

XLDE 

Carbitex  700 

6.6 

M 

ll 

Pyrocarb  901 

8.9 

II 

HTPB 

Carbitex  700 

9.5 

ii 

II 

Pyrocarb  901 

12.2 

Throat 

XLDB 

PG  Washers 

3.8 

II 

i<  ‘PB 

II 

5.9 

Downstream  Throat 

XLDB 

Carbitex  700 

2.1 

ll 

II 

Pyrocarb  901 

3.3 

<1 

HTPB 

Carbitex  700 

1.8 

ll 

II 

Pyrocarb  901 

5.5 
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TABLE  20.  PROPELLANT  DATA,  STUDY  4 


Propellant 

XLDB 
(19%  Al) 

HTPB 
(21%  Al, 
90%  Solids) 

PEG/FEFO 
(20%  Al) 

„ , ,GMAT0MS, 

Hydrogen  (1Q0  GMS) 

2.3885 

3.4864 

2.5416 

Carbon 

1.2348 

0.6933 

1.0967 

Nitrogen 

1.6779 

0.5933 

1.5528 

Oxygen 

2.4224 

2.3766 

2.2753 

Fluorine 

- 

- 

0.0750 

Aluminum 

0.7042 

0.7784 

0.7417 

Chlorine 

0.0423 

0.5873 

0.1320 

Pressure  (psia) 

1000. 

1000. 

1000. 

Temperature  (°R) 

6808. 

6699. 

6772. 

0.41 


i 


Layer  PC  or  SiC/PG 

(b)  Layerpd  throat  insert 


Pvrocarb  901 
or  edge  PG 


(c)  Solid  throat  insert 
Figure  45.  Concluded. 


9 


96trtrt  -V 


TABLE  21.  MATERIAL  RESPONSE  SUMMARY,  STUDY  4 


Propellant 

Throat 

Material 

Total 

Recession3 

(in) 

Average 

Recession  Rate 
(mils/sec) 

XLDB 

C/C  901 

0.444 

7.4 

Edge  PG 

0.203 

3.4 

Layer  PG 

0.119 

2.0 

15%  SiC  PG 

0.145 

2.4 

HTPB 

C/C  901 

0.588 

9.8 

Edge  PG 

0.289 

4.8 

Layer  PG 

0.104 

1.7 

15%  SiC  PG 

0 . 32  5b 

5.4 

PEG/FEFO 

C/C  901 

0.383 

6.4 

Edge  PG 

0-156 

2.6 

Layer  PG 

0.112 

1.9 

15%  SiC  PG 

0.095 

1.6 

a60-second  firing 

^Extrapolated  (coating  burned  through  at  59  secs.) 


propellant  group.  It  Is  also  evident  that  HTPB,  XLDB,  and  PE6/FEF0,  respectively  decrease  in  cor- 
rosivity for  all  materials  except  layer  PG. 

Study  5 - Standard  7-Inch  Muzzle  (Nosecap  Sensitivity)  (Reference  14) 

The  purpose  of  this  study  was  to  compare  the  predicted  performance  of  five  possible  nosecap 
materials  in  the  standard  7-inch  test  nozzle.  The  five  materials  are: 

• ATO  graphite 

• G-90  graphite 

e Pyrocarb  901  carbon/carbon  (density  = 1.83  gr/cc) 

• Carbitex  700  carbon/carbon  (density  1 1.50  gr/cc) 

• 15%  S1C/PG 

A one-dimensional  thermal  analysis  including  surface  recession  was  performed  for  each  mate- 
rial in  an  HTPB  propellant  environment.  This  propellant  is  identical  to  the  HTPB  used  in  Study  4, 
as  are  the  firing  duration  (60  seconds)  and  average  chamber  pressure  (1000  psia).  The  nozzle  geom- 
etry is  shown  in  Figure  47. 

Table  22  summarizes  the  predicted  total  recession  and  average  recession  rate  for  the  five 
materials.  A more  graphic  presentation  of  the  average  recession  rates  is  given  in  Figure  48.  This 
shows  that  the  materials  can  be  ranked  by  their  resistance  to  ablation  as  follows: 

1.  15%  S1C/PG 

2.  Carbitex  700  carbon/carbon 

3.  G-90  bulk  graphite 

4.  Pyrocarb  901  carbon/carbon 

5.  ATJ  bulk  graphite 

Study  6 - BATES  Nozzle  (Reference  15) 

Previous  studies  were  to  predict  tne  ablation  performance  of  rocket  nozzle  materials.  Since 
the  Air  Force  Rocket  Propulsion  Laboratory  (AFRPL)  has  experienced  numerous  failures  of  the  one 
piece  graphite  nozzle  utilized  in  various  versions  of  the  BATES  rocket  motor,  the  purpose  of  this 
study  was  to  perform  a preliminary  structural  analysis  of  the  BATES  nozzle  to  understand  the  reasons 
for  failure  and  to  recommend  solutions  to  the  problem. 
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Insulator 


TABLE  22.  MATERIAL  RESPONSE  SUWARY 


Propellant 

Material 

Total  Recession3 
(in) 

Average 
Recession  Rate 
(mils/sec) 

HTPB 

ATJ 

0.645 

10.8 

6-90 

0.505 

8.4 

Pyrocarb  901 

0.558 

9.5 

Carbitex  700 

0.456 

7.6 

15%  SiC/PS 

0.309 

5.2 

a60-second  firing 


Average 
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Three  basic  classes  of  BATES  motors  are  in  use  at  AFRPL:  a 15-pound  motor,  a 70-pound  motor, 
and  a high  pressure  (approximately  80-pound)  motor.  A drawing  of  tne  high  pressure  BATES  configura- 
tion is  shown  in  Figure  49.  The  nozzle  is  of  single  piece  construction,  specified  as  HLM  85  or 
equivalent  graphite.  The  nozzle  is  held  in  place  by  16  steel  bolts  torqued  to  600  inch-pound.  This 
preloads  the  graphite  between  a retaining  ring  and  the  motor  aft  closure.  An  0-ring  seal  prevents 
gas  leakage  at  the  entrance  region  of  the  nozzle. 

The  BATES  motors  exhibit  an  essentially  flat  pressure-time  characteristic,  reaching  full 
chamber  pressure  quite  rapidly.  Typical  run  times  are  only  a few  seconds,  with  9 or  10  seconds 
being  a maximum.  A worst  case  propellant  assumption  would  be  a high  solids  loaded,  high  aluminum 
content  formulation  such  as  90  percent  solids  - 20  percent  aluminum  formulation. 

Nozzle  failures  have  occurred  with  all  three  motor  types,  particularly  when  operating  at 
pressures  greater  than  1500  psia.  The  mode  of  failure  is  consistent,  involving  radial  fracture 
planes  as  a result  of  hoop  stresses.  Generally  the  nozzle  fractures  into  two  nearly  equal  pieces. 
However,  three  radial  fracture  planes  producing  three  essentially  120°  segments  are  sometimes  ob- 
served. Occasionally,  a transverse  fracture  is  also  observed  at  the  juncture  of  the  forward  nozzle 
flange  and  the  motor  aft  closure.  It  is  felt  that  this  fracture  may  occur  subsequent  to  the  pre- 
viously described  hoop  failure.  All  failures  are  felt  to  occur  quite  early,  possibly  upon  reaching 
full  chamber  pressure. 

Many  different  graphites  including  Great  Lakes  HLM  85,  Stackpole  2020,  Airco  Speer  8882  and 
873,  Union  Carbide  ATJ,  and  Carborundum  G83  have  failed.  Although  a grain  direction  is  not  speci- 
fied on  the  nozzle  insert  drawing,  it  is  safe  to  assume  that  the  grain  direction  was  perpendicular 
to  the  nozzle  centerline  for  molded  graphites  and  along  the  nozzle  centerline  for  extruded  graphites 

Per  AFRPL  request,  the  analysis  was  initially  concentrated  on  the  high  pressure  BATES  design. 
For  approximately  a dozen  firings  up  to  the  onset  of  this  study,  the  high  pressure  nozzle  had  ex- 
hibited approximately  a 75  percent  success  rate  for  pressures  up  to  2000  psia  and  essentially  no  sue 
cess  beyond  2000  psi. 

In  addition  to  the  nozzle  failures  occurring  during  actual  motor  firings,  a hydrostatic  test 
has  also  resulted  in  nozzle  failure.  A special  hydrostatic  test  nozzle  was  fabricated  of  HLM  85 
with  the  internal  features  of  the  exit  cone  left  unmachined  (solid).  The  nozzle  was  mounted  in  a 
fixture  identical  to  the  aft  closure  of  the  motor  and  pressurized  hydrostatically  to  failure  at  am- 
bient temperature.  Failure  in  two  nearly  equal  halves  identical  to  actual  firing  failures  occurred 
at  1700  psia. 
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The  general  approach  to  the  problem  consisted  of  first  conducting  a structural  analysis  of 
the  existing  BATES  high  pressure  configuration  in  both  the  nydrostatic  test  and  motor  firing  modes 
to  see  if  the  observed  failures  could  be  correlated.  Once  a correlation  was  established,  these 
same  analytical  techniques  were  used  to  examine  various  potential  fixes.  The  DOASIS  finite  element 
computer  code  was  used  (Reference  24)  to  solve  the  displacements,  strains  and  stresses  in  the  nozzle 
using  orthotropic  properties  ard  axi symmetric  pressure  and  thermal  loading.  All  motor  calculations 
were  performed  for  a chamber  pressure  of  2000  psia,  which  appears  to  be  the  threshold  level  for  re- 
latively consistent  nozzle  failure. 

As  a result  of  the  finite  element  calculations,  it  was  concluded  that  the  failure  of  the 
BATES  nozzle  probably  occurs  as  a result  of  circumferential  tensile  stresses  in  a region  just  for- 
ward of  the  nozzle  throat.  It  was  further  concluded  that  the  failure  occurs  very  early  in  the  burn 
and  is  almost  solely  the  result  of  pressure  loading  on  the  inner  contour  at  a time  when  the  inter- 
face leading  to  the  0-ring  is  essentially  sealed. 

Assuming  the  above  failure  mode  is  indeed  correct,  the  most  promising  fix  to  the  problem  in- 
volves intentional  design  features  to  assure  that  chamber  pressure  reaches  the  0-ring  rapidly.  In 
this  way  the  nozzle  is  in  effect  "self-healing".  A reduction  in  maximum  hoop  tensile  stress  of 
greater  than  a factor  of  two  is  realized  from  this  change  alone. 

Another  fix  with  the  potential  of  increasing  nozzle  capability  by  more  than  a factor  of  two 
is  the  use  of  ATJ  or  equivalent  strength  graphite  rather  than  HLM  85  quality  graphHe.  Several  cf 
the  qraphites  that  have  failed  in  the  BATCS  testinq  are  quite  inferior  to  ATJ. 

A slip-fit  sleeve  does  not  appear  promising  as  a fix  to  the  problem  for  two  reasons.  First, 
the  radial  expansion  of  the  baseline  design  with  no  sleeve  is  only  about  2 mils,  which  is  probably 
the  same  order-of-magni tude  as  the  gaps  present  in  a slip-fit  sleeve  of  this  large  diameter.  Sec- 
ond, the  calculations  for  a zero  gap  indicate  only  about  a 30  percent  reduction  in  hoop  tensile 
stress  due  to  the  sleeve. 

Pressurization  of  the  nozzle  00  appears  to  be  quite  beneficial  in  reducing  hoop  tensile 
stresses.  The  benefits  are  nearly  as  great  in  the  region  of  major  interest,  the  entrance/throat 
region,  if  pressurization  is  applied  only  to  the  forward  portion  of  the  0D. 

There  appears  to  be  a tendency  for  the  couple  produced  by  the  offset  action  lines  of  the 
net  pressure  force  and  the  reaction  force  at  the  aft  end  of  the  nozzle  to  splay  the  aft  end,  pro- 
ducing hoop  tensile  stresses.  Although  the  failures  are  not  felt  to  originate  in  this  region, 
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angling  the  aft  end  of  the  nozzle  somewhat  will  alleviate  tnis  situation.  Although  it  was  not  ex- 
amined, it  might  be  possible  to  design  the  forward  end  of  the  nozzle  in  a similar  fashion  so  that 
the  bolting  force  would  produce  a compressive  hoop  prestress  in  the  entrance  region. 

All  of  these  conclusions  were  presented  to  AFRPL  as  recommendations  for  solving  the  failure 
problems  of  the  high  pressure  BATES  nozzle.  Further  recommendations  were  also  made  to  solve  similar 
problems  on  the  lower  pressure  BATES  motors. 


SECTION  8 


CONCLUSIONS  AND  RECOMMENDATIONS 

Kinetic  expressions  for  ATJ,  G-90,  c oriented  PG,  15%  SiC/PG,  and  Pyrocarb  901  were  obtained 
from  correlations  of  APG  data  and  rocket  motor  firing  data.  These  expressions,  along  with  expres- 
sions from  previous  results  for  a-b  oriented  PG  and  a recorrelation  of  previous  Carbitex  700  data, 
were  incorporated  into  the  former  GASKET  code  for  thermochemical  ablation  analysis.  The  new  code, 
GASKE12,  also  includes  an  improved  accounting  of  sublimation  kinetics,  additional  gas  phase  species 
(Duff-Bauer),  and  a generalized  procedure  for  specification  of  kinetic  reactions. 

Using  the  GASKET2  code  and  the  Aerotherm  rocket  nozzle  prediction  procedure,  the  ablation 
performance  of  15  motor  firings  was  predicted.  All  but  three  of  these  predictions  were  within  +25 
percent  of  measured  results.  From  these  predictions,  which  included  a wide  range  of  propellants 

and  nozzle  designs,  it  was  concluded  that  the  kinetic  expressions  were  acceptably  accurate.  How- 

ever, additional  verifications  are  required  since  in  several  cases  only  one  motor  tiring  was  used 
to  validate  the  kinetic  expressions  for  a particular  material. 

In  addition  to  the  data  used  for  kinetic  correlations,  ablation  data  for  13  other  materials 
were  obtained  in  the  Aerotherm  APG.  These  data  were  compared  within  a generic  class  of  materials 
to  determine  qualitative  performances.  Data  obtained  to  date  show  that  some  materials  are  sensitive 
to  manufacturing  source  and  others  are  not. 

Significant  conclusions  from  this  investigation,  including  discussions  from  the  appendix, 

are: 

1.  The  GASKET2  code  and  the  Aerotherm  rocket  nozzle  prediction  procedure  can  accurately  pre- 
dict the  ablation  performance  of  a rocket  nozzle  which  uses  materials  that  have  been 

characterized. 

2.  If  propellant  temperatures  increase  beyond  current  values  of  about  6800°R,  carbon  sub- 
limation kinetics  will  become  important.  In  anticipation  of  this,  the  new  GASKET2  code 
includes  appropriate  sublimation  terms. 
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3.  At  surface  temperatures  representative  of  MX  nozzles,  the  ^ - C*  reaction  is  the  domi- 
nant carbon  removal  mechanism.  Although  the  reaction  rate  constants  are  large  for  1^0 
and  CC^,  their  reaction  rates  at  motor  temperatures  are  usually  small  because  of  low  re- 
actant partial  pressures. 

4.  Pyrolytic  graphite  ablation  rates  are  not  sensitive  to  the  material  suppliers. 

5.  Bulk  graphite  and  carbon/carbon  composite  ablation  rates  are  sensitive  to  the  material 
suppliers.  In  general,  density  alone  is  not  a performance  scaling  parameter;  however, 
for  carbon/carbon  composites  it  was  observed  that  a nonlinear  density  dependence  could 
exist,  if  the  material  is  separated  according  to  the  reinforcement  precursor.  That  is, 
rayon  precursor  composites  and  PAN  precursor  composites  each  exhibit  a density  dependence. 
It  is  highly  p.obable  that  the  matrix  precursor  should  also  oe  a correlating  parameter. 

6.  For  equal  densities,  rayon  precursor  composites  performed  better  in  the  APG  than  PAN 
precursor  composites.  As  shown  in  Appendix  A,  PAN  precursor  composites  show  evidence  of 
fiber  cleavage  whereas  rayon  precursor  composites  do  not. 

7.  At  very  high  temperatures  all  materials  seem  to  approach  an  asymptotic  mass  removal  rate. 

8.  No  correlations  of  ablation  performance  were  observed  for  material  porosity. 

9.  Chemadsorption  measurements  show  no  surface  adsorption  of  reaction  species,  namely, 

CO.,,  t^O,  and  H.,. 

10.  Significantly  more  data  and  data  analysis  will  be  required  to  relate  ablation  performance 
to  processing  or  material  microstructure. 

The  following  are  recommended  in  order  to  improve  rocket  nozzle  ablation  performance  predic- 
tions and  to  obtain  an  understanding  of  the  relationship  between  ablation  performance  and  material 
processing. 

1.  APG  data  should  be  obtained  for  additional  carbon/carbon  composites  which  include  a wide 
/•ange  or  processing  variables.  Ideally,  the  composites  should  be  processed  with  a sys- 
tematic variation  in  these  variables. 

2.  The  microstructure  of  each  material  should  be  carefully  examined  and  recorded.  This  ex- 
amination should  include,  as  a minimum,  metal lographs,  SEM's,  porosity,  and  chemoadsorp- 
tion. 
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3.  Thermophysical  properties  should  also  be  measured  and  include,  as  a minimum,  local 
density,  thermal  expansion  coefficients  and  thermal  conductivity. 

4.  Additional  motor  firings  should  be  examined  and  predicted  to  improve  the  prediction  ac 
curacy  of  the  GASKET2  code. 

5.  Alternate  procedures  should  be  examined  for  obtaining  accurate  high  temperature  kineti 
ablation  data. 
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APPENDIX  A 


MICROSTRUCTURAL  EXAMINATIONS 


A. I INTRODUCTION 

Surface  kinetic  constants  depend  on  the  microstructure  of  the  material.  This  was  shown  by 
the  large  difference  in  measured  ablation  rates  of  a-b  plane  and  c plane  pyrolytic  graphite.  Be- 
cause of  this,  parameters  such  as  relative  fiber/matrix  content,  composite  density  and  porosity, 
and  the  degree  of  graphitization  are  potentially  important  variables.  In  addition,  the  orientation 
of  carbon  platelets,  which  depend  on  the  precursor  materials  and  their  processing,  is  important. 

For  nstance,  PAN  and  rayon  precursors  yield  fibers  with  platelet  orientations  as  shown  in  Figure 
A-l.  For  ablation  along  the  cylindrical  face,  PAN  fibers  appear  like  layered  PG  and  rayon  fibers 
appear  like  random  carbon  edges.  Thus  the  ablation  kinetics  will  also  depend  on  the  relative  fabric 
orientation  and  the  orientation  of  the  fabric  relative  to  the  ablation  surface. 

Processing  of  carbon/carbon  rocket  nozzle  materials  is  quite  proprietary;  hence,  only  general 
information  is  available.  Some  of  the  information  for  materials  used  in  this  investigation  is  de- 
scribed in  Section  A. 2.  However,  greater  details  are  required  before  a meaningful  relationship  can 
be  established  between  materials  processing  and  ablative  performance.  Figure  A-2  shows  the  inter- 
dependence between: 

• Fabrication  and  processing 

• Pretest;  microstructure 

• Ablation  performance 

• Post-test  microstructure 

The  ideal  situation  is  to  be  able  to  obtain  a desired  ablation  performance  by  prescribing  fabrica- 
tion and  processing  variables.  Unfortunately,  the  current  state-of-knowledye  is  only  sufficient 
to  provide  simple  guidelines  which  may  have  many  qualifiers  or  exceptions. 

With  the  current  engineering  approach  to  determining  reaction  kinetics,  it  is  not  possible 
to  define  the  fundamental  mechanisms  for  carbon  consumption.  Thus,  it  would  not  be  unexpected  for 
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Desired  relationship 


Easily  obtained  relationships 
which  can  be  varied 


Links  to  be  determined  for 
material  gu'dance 


Figure  A-2.  Relationship  between  ablation  and  fabrication. 
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each  form  of  carbon  (e.g.,  bulk  graphite,  pyrolytic  graphite,  etc.)  to  be  represented  by  a different 
reaction  rate  expression.  This  behavior  has  in  fact  been  verified  by  this  investigation.  For  in- 
stance Figure  A-3  is  a theoretical  prediction  based  on  the  empirical  correlations  described  in  Sec- 
tion 6.  This  figure  shows  that  the  mass  removal  rate  is  dependent  upon  the  type  of  carbon  under 
consideration.  It  also  implies  that  preferential  attack  will  occur  in  carbon  composites  which  are 
heterogeneous  mixtures  of  different  carbons.  This  implication  has  been  verified  by  tl.e  ablation 
tests  described  in  Section  4.  Figure  A- 4 shows  the  results  of  some  of  these  measurements  for 
carbon/carbon  materials  in  a hydrogen  environment.  Different  precursor  materials  are  apparently 
attacked  more  or  less  vigorously  by  hot  hydrogen  gases.  Preferential  ablation  can  also  be  seen  in 
SEM's  of  post-test  surfaces.  For  example,  Figures  A-5  and  A-6,  respectively,  show  the  ablated  sur- 
faces of  rayon  and  PAN,  precursor  composites.  Rayon  is  apparently  attached  preferentially,  leaving 
a visible  surface  which  is  virtually  all  matrix  material.  PAN  composites,  however,  show  surfaces 
which  are  mixtures  of  fibers  and  matrix.  Significant  micromechanical  breakage  of  fibers  is  observed 
and  apparently  enhances  the  ablation  rates.  This  fiber  cleavage  is  more  evident  in  Figure  A-7. 

Two  tentative  conclusions  can  be  reached  from  a preliminary  examination  of  the  carbon/carbon 
pre-  and  post-tcsT 'microstructures.  First,  for  a resin  (or  pitch)  impregnated  c/c  at  a given  dens- 
ity, the  recession  rate  performance  of  rayon  precursor  reinforcements  is  superior  to  that  of  PAN 
precursor  reinforcements  (Figure  A-4).  Second,  post- test  SEM's  show  that  the  reinforcement  in  a 
rayon  precursor  c/c  recedes  below  the  surface  of  the  matrix  (Figure  A-5).  This  second  conclusion 
suggests  that  the  matrix  is  a better  ablator  than  the  rayon  reinforcement.  Thus,  a potentially  good 
material  would  use  a low  volume  content  of  rayon  precursor  reinforcement  and  a high  density  graph- 
itized  matrix. 

It  is  recognized  that  continuous  filament  rayon  material  will  soon  be  unavailable  so  that  a 
first  impression  would  be  to  avoid  the  use  of  rayon.  However,  there  is  a mounting  effort  to  develop 
low  density  PAN  and  pitch  precursor  reinforcements  with  physical  properties  similar  to  that  of  rayon. 
Hence,  basic  information  obtained  with  rayon  precursor  c/c  materials  will  be  valuabTe  for  providing 
ablation  and  microstructure  relationships  to  guide  material  development. 

Fiber  volumes  for  two-dimensional  c/c  materials  are  typically  45  to  55  percent.  Fiber  yolumes 
for  3-D  materials  are  around  30  to  35  percent.  Thus,  3-D  materials  offer  two  advantages.  First, 
they  are  structurally  better  than  2-D  materials  and  second,  they  have  inherently  lower  fiber  volumes. 
The  disadvantage  of  3-D  materials  is  their  higher  fabrication  cost  and  potentially  greater  scale-up 
problems.  If  3-D  materials  prove  to  be  ablation-wise  superior,  then  these  disadvantages  can  be 
tolerated. 
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Figure  A-3.  Relationship  between  mass  transfer  coefficient  and  ablation 
rate  for  various  carbon  materials. 
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Figure  A-4.  Measured  ablation  rates  in  hydrogen  environment,  T = 5500°R. 
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Figure  A-5 

Exposed  matrix  at  surface  of  rayon  precursor  composite. 
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Exposed  matrix  and  fibers  of  PAN  precursor  composite 
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Evidence  of  fiber  cleavage  in  PAN  precursor  composite  . 


As  part  of  this  investigation,  pre-  and  post-test  microstructural  characteristics  were  ex- 
amined for  two  purposes.  First,  it  was  hoped  that  some  obvious  relationships  would  be  observed  to 
relate  microstructure  to  ablation  performance  and  material  fabrication.  Believing,  however,  that 
there  would  not  be  sufficient  data  to  establish  firm  relationships,  the  second  purpose  was  to  col- 
lect this  microstructure  data  as  a "data  bank"  for  future  reference.  As  additional  data  is  added 
to  this  data  bank,  the  relationships  between  fabrication  and  ablation  performance  should  become 
more  obvious.  Eventually,  the  kinetic  behavior  of  a graphitic  material  can  be  determined  by  a de- 
tailed knowledge  of  its  fabrication  or  by  examination  of  its  pretest  microstructure.  This,  of 
course,  is  a very  ambitious  objective  and  may  take  a decade  to  realize. 

The  microstructural  characterization  measurements  performed  in  this  study  are  tabulated  in 
Table  A-l.  All  characterizations,  with  the  exception  of  ordinary  photography  and  scanning  electron 
microscopy,  were  performed  by  Mr.  Jay  Baetz  of  tne  Aerospace  Corporation,  El  Segundo,  California. 

A. 2 MICROSTRUCTURAL  EXAMINATION 

Microstructural  examination  was  conducted  to  characterize  both  virgin  and  tested  material 
conditions.  The  primary  areas  of  interest  in  this  examination  included  in-depth  assessment  of  un- 
tested material  structure  and  the  surfaces  resulting  from  exposure  to  arc  test  conditions.  The 
principal  microstructural  features  of  interest  for  tested  specimens  included  response  of  reinforce- 
ments ana  matrices  to  the  imposed  test  conditions. 

Selected  specimens  were  examined  by  photomicrographic  and  scanning  electron  microscopy  (SIM) 
techniques.  Representative  selected  materials  were  examined  in  tne  tested  and  uraosted  states  by 
both  techniques.  However,  primary  emphasis  in  tne  photomicrographic  studies  was  directed  toward 
untested  material  examination.  This  emphasis  was  maintained  so  that  high  resolution  definition  of 
polished  material  samples  would  be  obtained.  The  SEM  was  primarily  used  for  arc  tested  material 
analysis,  which  permitted  utilization  of  the  SF.M  large  depth  of  focus  to  assess  tested  surfaces 
without  surface  modifications. 

Materials  Summary 

Three  classes  of  materials  were  examined.  These  classes  included  carbon/carbon  composites, 
bulk  graphites,  and  vapor-deposited  carbons.  Within  the  carbon/carbons  both  two-  and  three- 
dimensional  ly  reinforced  materials  were  analyzed.  The  bulk  graphites  included  aerospace  grade, 
fine  grained  systems  and  commercial  grade  graphites.  Vapor  deposited  carbon  material  systems  ex- 
amined included  conventional  pyrolytic  graphites,  a-b  oriented  graphites  and  graphites  codeposited 
with  metal  carbides.  A summary  of  the  materials  examined  is  provided  in  Table  A-2. 
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TABLE  A-l . MICROSTRUCTURAL  CHARACTERIZATION  MEASUREMENTS 
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Chemadsorption 


TABLE  A-2.  SUMMARY  OF  MATERIALS  EXAMINED 


Material  Class 

Manufacturer 

Identi ty 

Carbon/carbon 

Hitco/Defense  Products  Division 

Pyrocarb  901 

Hitco/Defense  Products  Division 

Pyrocarb  903 

Haveg 

HRX  5125 

Haveg 

MRX  5875 

McDonnell -Douglas 

MDAC  3-D  c/c 

Carborundum 

Carbitex  700 

Bulk  Graphite 

Union  Carbide 

ATJ 

pure  Carbon 

PO-3 

Carborundum 

G-90 

Union  Carbide 

ATJ-S 

Vapor  Deposited  Carbons 

Pfizer 

Standard  pyrolytic  graphite 

Supertemp 

Pyrolytic  graphite  (101-B) 

Hitco/Defense  Products  Division 

a-b  oriented  pyrolytic  graphite 

Atlantic  Research  Corporation 

5%  SiC/pyrolytic  graphite 

Atlantic  Research  Corporation 

15%  SiC/pyrolytic  giaphite 

Atlantic  Research  Corporation 

23%  SiC/pyrolytic  graphite 

Raytheon 

HfC  coated  pyrolytic  graphite 
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These  materials  constitute  a representative  cross-section  of  all  classes  of  materials  cur- 
rently considered  for  nozzle  application.  Within  each  materials  class,  a number  of  constituents  or 
processing  techniques  are  considered  proprietary  by  the  respective  manaufocturers.  A summary  of 
each  material  examined  is  provided  within  the  above  cited  proprietary  restrictions  as  follows. 

Pyrocarb  901 

Pyrocarb  is  a two-dimensional  reinforced  composite  with  a rayon  precursor  square  weave  carbon 
fabric.  The  laminate  is  initially  formed  by  fabric  impregnation  with  a phenolic  resin.  Individual 
plys  are  rotated  within  the  billet  plane  as  specified  by  the  customer.  Details  of  processing  and 
densification  are  considered  proprietary  by  the  manufacturer. 

Pyrocarb  903 

This  material  is  also  a two-dimensional  reinforced  carbon/carbon.  The  reinforcement  is  an 
eight  harness  satin  weave  with  a PAN  precursor.  Laminates  are  fabricated  with  the  fabric  impregnated 
with  a phenolic  matrix.  Multiple  density  levels  are  available  depending  upon  the  customer  require- 
ments. In  this  program  two  density  levels  were  tested.  Ply  rotation  is  as  specified  by  the  customer. 

HRX  5125 

HRX  5125  is  a two-dimensional  carbon/ carbon  fabricated  from  a rayon  precursor  fabric  such  as 
WCA.  The  initial  densi fication  is  accomplished  with  H-resin.  This  matrix  is  carbonized  at  800°C. 
Redensification  is  accomplished  with  the  H-resin  for  an  unspecified  number  of  cycles.  Subsequent 
densification  is  accomplished  with  Allied  Chemical's  new  coal  tar  pitch.  The  final  processing  step 
performed  is  at  2800°C  graphitization.  Adjacent  plys  are  rotated  to  obtain  quasi-isotropic,  in-plane 
properties. 

HRX  5875 

HRX  5875  is  a two-dimensional  reinforced  carbon/carbon  composite.  The  reinforcement  consists 
of  a PAN  precursor  carbon  fabric  with  an  eight  harness  satin  construction.  A typical  yarn  used  is 
Thomel  300.  Prior  to  use  in  laminates,  the  fabric  is  heat  treated  to  the  graphitization  tempera- 
ture. K-resin  is  used  as  the  initial  impregnating  matrix  and  is  carbonized  at  800"C.  At  the  mid- 
point of  densification  processing,  tne  laminate  is  graphitized.  Continued  densification  is  accom- 
plished with  15  V coal  tar  pitch.  The  final  processing  step  is  graphitization  at  2800°C.  Ply  rota- 
tion is  performed. 
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Carbitex  700 


Carbitex  700  is  a two-dimensional  reinforced  carbon/ carbon.  The  reinforcement  is  a square 
weave  rayon  precursor  graphitized  in  fabric  form.  The  matrix  utilizes  a s i n precursor.  Graphiti- 
zation  is  conducted  at  a minimum  temperature  of  2700°C. 

MDAC  3-0  C/C 

The  single  three-dimensional 1y  reinforced  material  was  MDAC  3-D  c/c.  This  composite  utilized 
orthogonal  reinforcements  of  Thornel  50  yarns.  Weave  construction  was  224.  The  prefonn  was  woven 
by  FMI  and  carried  the  identity  of  173B.  Densification  was  accomplished  with  15  V pitch  and  a 
thermosetting  resin.  Graphitization  between  densification  cycles  was  conducted  at  2800“C.  The 
final  thermal  treatment  was  a 1000°C  pyrolysis.  No  definition  of  the  pitch/resin  impregnating  se- 
quence is  available. 

ATJ-S 

ATJ-S  was  the  initial  bulk  graphite  considered.  This  product,  the  best  characterized  of 
the  aerospace  grade  graphites,  is  fine  grained  (0.006  inch  maximum)  with  a density  of  1.83  gm/cm3 
and  an  ash  level  of  193  ppm. 

ATJ 

ATJ  is  a molded  graphite  with  a grain  size  equivalent  to  an  ATJ-S.  However,  its  density 
(1.74  gm/cm3)  is  somewhat  lower  and  its  ash  level  (1200  ppm)  is  somewaht  higher  than  ATJ-S.  The 
specific  material  utilized  ir.  this  program  was  obtained  from  the  center  of  a 15-inch  diameter  billet. 

£0-3 

P0-3  is  a commercial  grade  molded  graphite,  stated  by  the  manufacturer  to  be  somewhat  porous. 
Specific  information  regarding  manufacturing  processes  has  not  been  obtained.  A number  of  grades 
are  produced  including  some  with  carbide  additions. 

G-90 

G-90  is  an  aerospace  grade  graphite.  It  has  a 1.9  gm/cm3  minimum  density  and  is  manufactured 
as  an  extruded  product.  Ash  content  is  0.06  percent.  The  maximum  grain  size  is  0.037  inch. 
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Supertemp  PG 


Pyrolytic  graphite  was  obtained  from  Supertemp  in  accordance  with  their  specification  101-B. 
No  information  as  to  source  gases,  deposition  temperatures,  or  other  processing  conditions  were 
available. 

HfC/PG 

This  product  consists  of  pyrolytic  graphite  codeposited  with  hafnium  carbide.  Codeposition 
•'■’formed  at  approximately  10  Torr  pressure  and  a temperature  between  1800°C  to  1850°C.  The  re- 
sulting deposition  rate  is  approximately  10  mils/hour  on  a graphite  mandrel  substance.  The  result- 
ing product  is  characterized  as  reasonably  fine  grain  renucleated. 

A-B  PG 

This  product  was  obtained  by  deposition  of  pyrolytic  graphite  on  premachined  ATJ  substrates. 
The  material  was  not  machined  prior  to  testing  so  that  the  initial  ablation  surface  was  the  as- 
deposited  surface. 

Pfizer  PG 

Pfizer  pyrolytic  graphite  is  deposited  as  pure  carbon  to  obtain  a continuously  nucleated 
structure.  Further  details  are  not  known. 

5,  15,  and  23  Percent  SiC  PG 

Silicon  carbide  is  codeposited  with  graphite  to  obtain  this  product.  Fabrication  information 
is  contained  in  Reference  A-l. 

A. 2.1  Metal  lographs 

Representative  samples  of  pre-  and  post-test  materials  were  photographed  using  a metal lograph 
These  metal lographs  yielded  qualitative  information  on  the  structure  and  uniformity  of  each  material 
Select  post- test  samples  were  also  photographed  to  reveal  the  roughness  of  the  ablation  surface. 
Magnifications  of  40  to  1000X  were  used  although  the  preponderance  of  photographs  were  for  magnifica 
tions  less  than  200.  A list  of  the  materials  and  magnification  levels  is  shown  in  Table  A-3. 

No  noticeable  changes  were  observed  for  the  in-depth  microstructure  between  pre-  and  post- 
test states.  Some  differences  were  noted,  however,  in  the  characteristic  surface  roughness  of  dif- 
ferent materials.  Since  surface  roughness  metallographs  also  show  the  in-depth  structure,  there 
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TABLE  A- 3.  PRE-  AND  POST-TEST  METALLOGRAPHS 
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would  be  little  value  in  showing  pretest  microstructures.  Thus,  only  representative  post-test 
metal lographs  will  be  shown. 

Metallographs  for  the  bulk  graphites,  PO-3  and  G-90  are  shown  in  Figures  A-8  and  A-9,  re- 
spectively. These  photographs  show  the  effect  of  grain  size  on  the  texture  or  roughness  of  the 
ablated  surface.  PO-3,  being  a very  fine  grain  material,  has  a virtually  smooth  ablation  surface, 
while  G-90,  which  has  a large  grain  size,  has  a rough  surface. 

ATJ  and  ATJ-S  have  grain  sizes  between  PO-3  and  G-90  and  although  not  shown  have  an  ablation 
surface  with  a roughness  between  that  of  PO-3  and  G-90. 

Metallographs  for  Pyrocarb  903  and  HRX  5875  are  shown  in  Figures  A-iO  and  A- 1 1 , respectively. 
These  are  both  PAN  precursor  reinforcements;  however,  there  are  some  obvious  differences  in  fabrica- 
tion. The  HRX  5875  apparently  uses  a heavier  carbon  yarn  or  is  processed  at  lower  pressures.  This 
is  evidenced  by  the  uniformity  of  the  plys  compared  with  those  of  Pyrocarb  903,  The  latter  shows 
significantly  more  yarn  distortion.  Surface  roughness  of  these  two  materials  are  about  the  same 
with  roughness  heights  o^  about  3 to  5 mils. 

Figures  A-12  and  A-13  show  the  microstructure  of  Pyrocarb  901  and  HRX  5125,  respectively. 

The  differences  in  ablation  surface  roughness  ara  apparent  although  both  are  fabricated  from  rayon 
precursor  fabric.  The  "evenness"  of  the  Pyrocarb  901  surface  is  unusual  when  compared  with  the 
other  carbon/carbon  composites. 

A metal lograph  for  the  3-D  carbon/carbon  is  shown  in  Figure  A- 14 . The  carbon  yarns  are 
clearly  thicker  than  those  of  2-D  composites  and  because  of  the  orthogonal  weave,  there  are  distinct 
and  systematic  matrix  pockets. 

A. 2. 2 SEM  Analysis 

As  stated  above,  primary  emphasis  in  the  SEM  studies  was  placed  on  post-test  examination  of 
the  eroded  nozzle  surfaces.  Sections  were  removed  from  the  tested  models  and  viewed  as  shown  in 
Figure  A-15.  In  all  SEM  photos  of  tested  models  taken,  the  entrance  of  the  throat  is  to  the  left 
of  the  picture.  All  specimens  were  examined  without  surface  preparation.  A number  of  magnification 
levels  were  taken;  however,  only  50  and  500  magnification  levels  will  be  shown.  As  with  the  metal lu- 
graphs,  only  representative  SEM  photographs  will  be  shown. 

A summary  of  the  models  and  test  conditions  for  which  SEM  analysis  was  conducted  is  provided 
in  Table  A-4. 
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f 'i]ure  A- 10.  Post-test  microphotoqraphs  for  Pyrocarb 
903  1 00X . 
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Fiyure  A - 1 1 . Post-test  microphoto<jraphs  for 
HRX  58/S  100/ 
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Figure  A-15.  SEM  sample  orientation. 
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TABLE  A-4.  SUMMARY  OF  BULK  GRAPHITES  EXAMINED  BY  SEM 
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TABLE  A-4.  Concluded 


Vapor  Deposited  Carbons 


SEM  photographs  of  Supertemp  pyrolytic  graphite  specimens  tested  in  hydrogen,  hydrogen/ 
oxygen,  and  hydrogen/oxygen/hydrogen  chloride  mixtures  are  shown  in  Figures  A-16,  A-17,  and  A-18. 

The  low  magnification  hydrogen  tested  model  exhibited  three  primary  features.  These  features  in- 
cluded the  formation  of  cavities,  microcracks  perpendicular  to  the  throat  axis  and  apparent  waviness 
of  structure  as  shown  in  Figure  A-16(a).  The  waviness  is  also  evident  at  higher  magnification 
(Figure  A- 16(b))  as  is  the  presence  of  a nodular  structure  on  the  cavity  walls. 

Exposure  of  Supertemp  pyrolytic  graphite  to  the  hydrogen/oxygen  test  gas  at  somewhat  lower 
temperatures  resulted  in  the  grainy,  near  fibril  appearance  shown  in  Figure  A-17(a).  At  high 
magnification,  these  features  are  plainly  shown.  As  indicated  in  Table  A-4,  a lower  recession  rate 
was  found  with  the  hydrogen/oxygen  mixture  than  the  pure  hydrogen  environment.  This  apparently 
anomalous  condition  may  be  related  to  the  lower  test  temperatures  employed  in  the  hydrogen/oxygen 
environment. 

The  structure  resulting  from  hydrogen,  oxygen,  and  hydrogen  chloride  closely  represents  the 
more  severe  recession  conditions  normally  seen  in  oxygen-containing  environments  on  other  materials. 
This  is  shown  in  Figure  A-18.  It  should  be  noted  that  the  temperature  employed  in  this  test  approx- 
imated that  used  on  other  oxygen-containing  environments.  Small  cracks  are  again  evident  in  a 
generally  scalloped  macrostructure.  At  higher  magnification,  the  generally  layered  structure  is 
evident  outside  the  crack  region. 

The  a-b  pyrolytic  graphite  in  a hydrogen  environment  exhibited  the  presence  of  partially  de- 
laminated layers  parallel  to  the  throat  and  the  presence  of  a nodular  structure  at  low  magnifica- 
tions. This  is  shown  in  Figure  19(a).  At  higher  magnification,  the  detail  nodular  structure  is 
apparent.  Some  alignment  of  the  nodule  peaks  parallel  to  the  test  gas  flow  is  also  observed  in 
Figure  A-19(b). 

Compared  to  the  Supertemp  material,  the  Pfizer  pyrolytic  graphite  material  exhibited  a 
similar  but  more  extreme  case  of  cavity  formation  in  hydrogen  environments.  Surface  pitting  and 
a-b  plane  delaminations  were  also  more  severe. 

The  codeposited  pyrolytic  graphites  are  typified  by  SiC/PG.  Figures  A-20  and  A-21  illustrate 
the  structures  observed  for  the  23 % and  5%  silicon  carbide  materials.  A definite  nodular  structure 
is  evident  at  the  low  magnification  in  both  the  23%  and  5%  SiC  specimens  tested  in  hydrogen  (Figure 
A-20(a)  and  A-21 (a) ) . Additionally,  these  figures  show  the  cavity  structure  observed  on  other  pyrol- 
ytic graphite.  A higher  magnifications,  small  cracks  at  the  base  of  the  nodules  are  evident  in  the 
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Figure  A- 17.  SEM  of  Supertemp  PG  tested  in 
H2/02,  Tw  = 4il2°R. 
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500X 


Figure  A-18.  SEM  of  Supertemp  PG  tested 
in  H2/02/HCJt,  Tw  = 5095°R. 
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Figure  20.  SEM  of  l t>X  Sio/PG  tested  in 
H2,  Tw  = 5478°R. 
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23%  silicon  carbide  material  (Figure  A-20(b)).  In  both  the  5%  and  23%  SiC  materials,  the  high 
magnification  shows  the  nodules  to  have  a fine  grained  appearance  (Figures  A-20(b)  and  A-21(b)). 

For  both  the  5%  and  23%  SiC  specimens,  introduction  of  oxygen  produced  a significantly  dif- 
ferent appearance.  At  low  magnification,  the  23%  silicon  carbide  model  formed  a white  coating 
(probably  S^02)  interrupted  by  holes.  At  high  magnification,  this  coating  is  observed  to  consist 
of  ceramic  type  globules.  The  5%  SiC  model  exhibited  a similar  coated  nodul-. r structure;  however, 
coating  density  was  less  complete.  Although  tested  at  almost  identical  temperatures,  high  magnifi- 
cation views  of  the  5%  silicon  carbide  coating  showed  a fiberous,  needle-like  structure  compared  to 
the  globular  structure  for  the  23%  material. 

Hafnium  carbide  pyrolytic  graphite  tested  in  hydrogen  gas  showed  a relatively  dense  uniform 
coating.  Presence  of  Hafnium  in  this  coating  is  indicated  by  darkened  areas  whicn  changed  while  in 
the  SEM  unit.  At  high  magnification,  minor  evidence  of  a substrate  structure  was  evident  below  the 
coating  cracks. 

Bulk  Graphites 

ATJ-S  is  representative  of  the  bulk  graphites  considered  in  the  SEM  analysis.  Introduction 
of  oxygen  in  the  test  gas  results  in  progression  of  the  microstructure  from  a comparatively  smooth, 
fine-grained  appearance  to  a roughened,  coarse  structure.  This  is  shown  in  Figures  A-22  and  A-23. 
Roughening  corresponded  to  a large  recession  rate  increase  as  indicated  in  Table  A-4.  This  change 
was  not  due  to  temperature  differences  used  in  the  test.  This  is  uased  upon  the  increased  erosion 
rate  and  coarsened  50X  microstructure  which  was  obtained  at  the  approximately  600°R  lower  H,,/02 
temperature.  The  coarsened  microstructure  is  retained  if.  the  H2/02  test  environment  with  an  increase 
of  1000°R  but  with  less  recession  rate  increase  than  the  change  in  test  gas  from  H2  to  H2/09.  No 
significant  change  in  500X  microstructure  is  evident  with  introduction  of  oxygen  into  the  environ- 
ment as  shown  comparing  Figures  A-22  and  A-23.  No  significant  changes  occurre  i either  in  macro  or 
microstructure  with  C0/HC5,  introduction. 

The  effect  of  oxygen  was  also  observed  in  somewhat  larger  grain-sized  graphites.  For  example, 
ATJ  progressed  from  an  essentially  uniform,  relatively  smooth  macrostrrcture  to  the  coarse,  rough 
structure  with  the  introduction  of  oxygen.  Again,  no  significant  change  in  500X  microstructure  was 
observed. 

The  PO-3  graphite  presented  a unique  feature  among  the  molded  or  extruded  graphites.  This 
feature  was  the  formation  of  macro  level  cavities  within  the  hydrcgen  environment.  However,  at  500X 
no  significant  difference  in  macrostructure  is  noted  compared  tc  the  ATJ-S  on  ATJ  systems. 
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Figure  A-23.  SEM  of  ATJ-S  in  Ho/0o, 
Fw  * '1258°R.  2 2 
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Carbon/Carbons 


The  progressive  microstructure  change  in  carbon/carbons  from  untested  material  through  expo- 
sures in  typical  environments  was  demonstrated  by  Pyrocarb  903.  This  progression  is  shown  in 
Figures  A-24  through  A-27  for  Pyrocarb  903.  Progression  from  untested  material  to  a hydrogen  ex- 
posed material  results  in  a moderate  amount  of  recession.  This  recession  is  evidenced  by  the  initial 
observation  of  peaks  and  valleys  around  the  PAN  based  yarns  parallel  to  the  film  plane  (Figures 
A-24(a)  aid  A-25(a)).  At  higher  magnifications,  the  structure  changes  from  thoroughly  intermixed 
fibers  and  matrix  to  exposure  of  some  filaments  with  a majority  of  the  matrix  retained. 

Introduction  of  oxygen  results  in  much  higher  yarn  definition  as  shown  in  Figure  A-26.  This 
definition  resulted  from  the  approximately  threefold  increase  in  ablation  rate  (see  Table  A-4).  It 
should  be  noted  that  this  change  with  H2/ 02  introduction  occurred  at  a slightly  lower  temperature 
than  H2  only.  The  ply  rotation  characteristics  of  this  material  is  shown  by  filaments  approximately 
perpendicular  to  the  film  plane  both  in  the  low  and  high  magnifications  of  Figure  A-26.  As  observed 
in  Figure  A-26(a),  introduction  of  oxygen  resulted  in  removal  of  significant  amounts  of  matrix  from 
the  filament  groups. 

HRX  5125  was  examined  as  another  representative  two-dimensionally  reinforced  carbon/carbon. 
This  rayon  precursor  reinforced,  H-resin/15  V pitch  densified  material  exhibited  the  same  dependence 
on  oxygen  presence  in  appearance  as  did  the  preceding  material.  All  of  the  features  previously  cited 
including  fiber  exposure,  ply  rotation,  and  matrix  loss  with  oxygen,  were  observed  to  some  degree. 

HRX  5785  presented  an  anomaly  to  the  above  cited  microstructure  trend.  Although  tested  at 
similar  or  lower  temperatures,  the  specimen  exposed  to  hydrogen  showed  an  apparently  more  severe 
loss  in  both  low  and  high  magnifications. 

The  single  three-dimensionally  reinforced  material  examined  was  MDAC  3-D  c/c.  This  material 
was  examined  at  three  locations  for  the  specimen  tested  in  hydrogen.  Two  of  these  positions  were 
orthogonal  while  the  third  was  at  45°  to  the  previous  two.  No  relationship  to  the  throat  circum- 
ferential position  was  noted.  In  all  cases,  the  surface  exhibited  distinct  fiber  exposure  and  re- 
tained matrix  in  the  low  and  high  magnifications.  Progressive  fiber  exposure  and  smoothness  were 
evident  in  the  H^/O 2 specimen. 

Although  tested  at  approximately  the  same  temperature  as  the  l'2/02  environment,  introduction 
of  CO  and  HCJt  resulted  in  better  retention  of  the  fibers  parallel  to  the  film  plane.  This  is  ac- 
counted for  by  the  CO  and  HC£  producing  a less  severe  environment  than  C2/02  and  was  confirmed  by 
the  reported  recession  rates. 
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Figure  A-24.  SCM  of  Pyrocarb  903,  as 
machined. 
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From  the  standpoint  of  microstructure,  changes  in  density  for  the  same  material  (over  the 
density  levels  tested)  provide  no  significant  alterations.  This  is  determined  by  comparing  Pyrocarb 
903  HD  (high  density)  with  the  Pyrocarb  903  previously  discussed.  Both  were  tested  in  Hg/O^  environ 
merits  at  similar  temperatures.  At  low  magnifications  the  same  roughened,  exposed  fiber  appearance 
is  evident.  At  high  magnifications  loss  of  matrix  around  filaments  is  evident  for  both  materials. 
Somewhat  more  matrix  retention  may  be  present  in  the  “HD"  version,  and  this  may  be  when  measured 
recession  rates  for  the  HD  version  were  slightly  lower  than  the  standard  Pyrocarb  903  version. 

The  roughened  surface,  fiber  exposure,  ply  rotations,  and  matrix  removal  features  were  re- 
tained with  the  introduction  of  HCt  into  the  H2/O2  test  gas.  It  should  be  noted  that  the  depth  of 
recession  between  yams  is  somewhat  less  for  this  case  than  for  H2/O2  only,  which  is  in  accordance 
with  the  lower  recession  rate  reported. 

Pyrocarb  901  showed  the  same  general  trends.  The  reported  recession  rates  in  all  test  en- 
vironments were  higher  for  this  material  than  for  Pyrocarb  903.  It  should  be  noted,  however,  for 
Pyrocarb  901,  definite  fiber  exposure  and  ply  exposure  was  observed  in  the  hydrogen  test  environment 
Introduction  of  oxygen  resulted  in  severe  removal  of  the  exposed  yarns.  This  removal  is  evidenced 
by  a "smoothing"  of  the  rayon  based  yarns. 

Examination  of  the  high  magnification  views  of  Pyrocarb  901  in  both  H2  and  l^/Og  is  somewhat 
less  informative.  Here,  matrix  was  retained  intermixed  with  the  fibers  in  both  environments. 

A. 2.  3 Thermal  Expansion  Coefficients 

Materials  and  specimens  on  which  thermal  expansion  coefficients  were  determined  am  summar- 
ized in  Table  A-5.  The  specimens  have  been  divided  into  the  following  three  material  classes: 

• Carbon/carbon  (c/c) 

• Pyrolytic  graphite  (PG) 

• Bulk  graphite 

Except  as  noted  below,  all  the  specimens  listed  in  Table  A-5  were  machined  out  of  discs  of 
approximately  2-1/4-inch  diameter  and  1/4  inch  thickness.  The  blackened  portion  of  sketches  drawn 
in  Column  7 of  Table  A-5  shows  the  location  of  specimen  within  the  disc  and  the  specimen  shape  as 
it  was  submitted  for  testing.  Specimen  L-l  was  machined  out  of  a G-90  bille*  and  its  location  in 
the  billet  is  shown  in  Figure  A-28. 

In  the  following  discussion,  c direction  refers  to  the  cross-ply  direction  and  a-b  to  the 
with-ply  direction. 
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TABLE  A-5.  TEST  MATERIALS 


No. 


Sample 

ID 


R-l 


Material 


No. 


26 


Name 


HRX  5125 
(Haveg) 


Class 


c/c 


Density 

(g/c/c) 


1.86 


Specimen 

Made 

From 

+ 


Disc 


© 


L-2 


34 


MDAC  3-D 


c/c 


Disc 


L-3 


D-8 


26 


HRX  5875 


c/c 


1.86 


28 


Pyrocarb  901 
(Hitco) 


c/c 


1.80 


Disc 


Disc 


D-9 


28 


Pyrocarb  903 


c/c 


1.84 


Disc 


® 


D-10 


D-3 


J-7 


26 


HRX  5125 
(Haveg) 


c/c 


12 


Supertemp  PG 


PG 


28 


Hitco  PG 


PG 


Disc 


Disc 


Disc 


® 


L-l 


41 


G-90 


Bulk  Graphite 


Billet  G-90 
(Figure  1 ) 
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Carbon/Carbon  Materials 


Thermal  expansion  data  as  a function  of  temperature  for  two-dimensional ly  reinforced  carbon/ 
carbon  materials  tested  is  summarized  in  Figure  A-29.  Thermal  expansion  is  lower  in  the  a-b  direc- 
tion, which  is  a result  of  carbon  fibers  not  expanding  as  much  as  the  matrix  material.  In  the  c- 
direction  at  a temperature  of  1800°F,  the  highest  value  of  thermal  expansion  is  9.4  x 10'3  inch/inch 
for  the  Pyrocarb  903  material.  The  corresponding  lowest  value  is  7.1  x 10'3  inch/inch  for  HRX  5125. 

In  the  a-b  direction  at  a temperature  of  1800°F,  the  highest  value  is  3.2  x 10"3  inch/inch 
for  HRX  5125  and  the  corresponding  lowest  value  is  0.02  x 10'3  inch/inch  for  HRX  5875.  The  increase 
in  the  thermal  expansion  values  is  approximately  linear  with  temperatures  over  the  temperature 
range  shown. 

Thermal  expansion  data  for  the  single  3-0  c/c  material  is  shown  in  Figure  A- 30 . As  expected, 
the  values  of  thermal  expansion  in  the  z-direction  are  considerably  lower  than  those  in  the  x-  or 
y-direction  since  the  z-direction  has  a higher  volume  of  c/c  fibers.  There  is  a nonlinear  relation- 
ship between  the  thermal  expansion  and  temperature  at  higher  temperatures  for  this  material. 

Vapor  Deposited  Carbons 

For  pyrolytic  graphite,  the  thermal  expansion  values  are  much  higher  in  the  c direction  than 
in  the  a-b  direction.  This  is  shown  in  Figure  A-31.  Both  specimens  (Hitco  a-b  PG  and  Supertemp  PG) 
exhibit  similar  thermal  expansion  characteristics  in  their  respective  c and  a-b  directions.  However, 
the  increase  in  the  thermal  expansion  values  with  temperature  in  the  c direction  is  much  greater 
than  that  in  the  a-b  direction.  The  increase  in  the  a-b  direction  expansion  values  becomes  nonlinear 
with  temperatures  above  about  3600°F. 

Bulk  Graphites 

The  only  thermal  expansion  measurements  of  bulk  graphites  were  for  G-90.  These  results  are 
shown  in  Figure  A- 32.  Some  anisotropy  in  this  material  is  reflected  in  the  two  different  curves. 

This  is  a result  of  the  preferential  grain  orientation  induced  by  processing  forces. 

A . 2 . 4 Ion  Microprobe  Mass  Analysis  (IMMA) 

IMffA  data  was  obtained  for  representative  materials  in  pre-test  and  post-test  states.  Post- 
test samples  were  selected  on  a basis  of  the  measured  surface  temperature  and  exposure  gas.  Ideally 
these  conditions  would  be  representative  of  those  found  in  rocket  nozzles,  however,  these  conditions 
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Figure  A-29.  Thermal  expansion  for  2-0  carbon/carbons. 
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Figure  A-26.  SEM  of  Pyrocarb  903  tested 
H9/0o,  T , - 5309°R. 
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could  not  be  duplicated  in  the  APG.  The  IMMA  post-test  samples  were  therefore  selected  as  a compro- 
mise between  test  gas  and  high  surface  temperatures.  Identification  information  on  the  selected 
IMMA  samples  is  shown  in  Table  A-6. 

IMMA  scans  were  made  for  both  positive  and  negative  ions.  Typical  surface  scans  are  shown 
in  Figures  A-33  and  A- 34.  Since  only  the  molecular  weight  of  molecules  are  shown,  it  is  not  always 
possible  to  unequivically  identify  the  molecule.  For  instance,  a molecular  weight  at  28  may  be  a 
molecule  of  CO,  or  Si.  In  addition,  the  concentration  of  each  surface  species  can  not  be  de- 
termined from  the  IMMA  data.  Thus,  the  scans  are  largely  qualitative. 

Rather  than  presenting  all  of  the  IMMA  scans,  the  dominant  lines  from  each  scan  were  identi- 
fied and  are  presented  in  tabular  form  in  Table  A-7.  Also  shown  in  this  table  are  possible  chemical 
species  that  could  be  represented  by  each  molecular  weight.  Only  molecular  weights  up  to  50  are 
shown. 

Two  samples  of  one  of  the  materials,  namely  Pyrocarb  901,  were  scanned  to  gain  an  apprecia- 
tion of  the  repeatability  of  the  material.  These  samples  were  at  slightly  different  reported  densi- 
ties but  were  presumably  otherwise  identical.  The  positive  ion  scans  showed  a noticeable  difference, 
which  may  be  attributed  to  Impurities  from  handling  or  to  material  processing.  The  former  appears 
more  likely  since  most  of  the  differences  are  from  elements  of  salts  which  are  contained  in  human 
perspiration. 

Comparisons  of  post-test  surface  scans  with  (1)  pre-test  scans  and  (2)  post-test  subsurface 
scans  show  that  the  surface  may  contain  a large  number  of  hydrocarbon  species.  These  would  be 
natural  by-products  of  the  reaction  between  APG  gases  and  the  surface.  A known  poison  species, 

HCt,  is  observed  only  in  one  case,  whereas  Cl  is  observed  in  almost  all  materials.  This  again 
may  possibly  be  attributed  to  sample  handling. 

From  the  amount  and  quality  of  data  taken,  no  firm  conclusions  can  be  reached  abojt  the  value 
of  the  IMMA  scans,  either  in  terms  of  guidance  for  material  development  or  as  aids  in  interpreting 
kinetic  reaction  rates.  However,  it  can  be  concluded  that  future  samples  must  be  handled  much  more 
carefully  and  the  IMMA  scans  should  be  conducted  very  shortly  after  testing.  New  procedures  should 
also  be  used  to  identify  adsorbed  species.  In  addition,  future  samples  should  be  scanned  at  a num- 
ber of  positions  to  reduce  the  uncertainties  of  local  impurities  or  contamination. 

A. 2 5 Porosity 

A mercury  porosimetry  was  used  to  determine  pore  size,  pore  size  distrioution  and  density 
jf  selected  materials.  Data  obtained  by  the  mercury  penetration  technique  includes  pore  size  and 
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TABLE  A-6.  IMMA  SAMPLES 


Pre- 

test 

9 

Post- 

test 

Exposure  Conditions 

Material 

Identification 

Tw 

(°R) 

Gas 

Time 

(sec) 

Recession 

(in.) 

Pyrocarb  901  (J2) 

Hitco  919946-2 

X 

1.72 

Mi 

Pyrocarb  901  (Jl) 

Hitco  919950-6 

X 

1.84 

■ 

Pyrocarb  901  (PA4) 

Hitco  919951-1 

1.83 

X 

5000 

1 

36.5 

0.0334 

Pyrocarb  903  (PA4) 

Hitco  919957-3 
Aerotherm  7113-134 

X 

1.84 

MDAC  3-0  (C2) 

X 

2.02 

G-90 

Aerotherm  41 
Layer  3 , 547  108° 

X 

SG-90 

Aerotherm  7113-170 

X 

4850 

7 

48.5 

0.024 

15%  SiC/PG 

ARC  8328-10 
003-23 

15%  SiC/PG 

Aerotherm  7113-116 

X 

4350 

1 

54.0 

0.024 

TABLE  A-7.  RESULTS  OF  IMMA  SCANS 
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volume.  Density  is  calculated  from  the  results  obtained.  These  data  are  determined  by  measuring 
the  quantity  of  mercury  which  can  be  forced  into  the  pores  of  the  test  material  at  various  increas- 
ing pressures.  Evidence  of  pore  shape  is  obtained  by  determining  the  amount  of  mercury  expelled 
from  the  pores  at  various  decreasing  pressures. 

The  theoretical  basis  for  the  mercury  penetration  method  is  negative  capillary  action:  this 

effect  results  from  the  nonwetting  nature  of  mercury.  Since  mercury  exhibits  a contact  angle  of 
greater  than  90°  with  most  materials,  it  will  not  penetrate  pore  openings  in  a material  unless 
forced  by  an  applied  pressure.  Mercury  will  then  penetrate  the  pores  in  accordance  with  the  force 
applied  and  the  size  o*  the  openings. 

The  relationship  describing  penetration  of  mercury  into  circular  openings  under  pressure  is: 

PD  3 -4o  cos  0 (A-l) 

where: 

P = applied  pressure 

D = diameter  of  smallest  pore,  filled  at  pressure  P 

a = surface  tension  of  mercury 

0 = contact  angle  or  wetting  angle 

Reported  measurements  of  contact  angle  between  mercury  and  a large  number  of  materials  range 
from  about  112°  to  142°,  with  a contact  angle  of  130°  as  the  most  frequent  value.  Accordingly, 
caking  0 = 130°  and  surface  tensin  of  mercury  o = 474  dynes/cm  (vacuum,  25°C),  the  following  rela- 
tionship is  obtained  for  cylindrical  pores: 


F0  = 177  (A-2) 

where  pressure  P is  in  psia  and  diameter  D is  in  microns. 

As  the  pressure  is  increased,  the  amount  of  liquid  mercury  forced  into  the  pores  increases. 
From  Equation  (A-2),  the  diameter  of  the  pore  is  obtained  for  that  particular  pressure.  The  volume 
of  mercury  forced  into  a pore  is  a direct  function  of  the  volume  of  the  pore.  Thus  a penetration- 
volume  versus  applied-pressure  curve  can  be  driwn  and  analyzed  for  particle  size.  Material  density 
can  be  calculated  with  pores  larger  than  any  limiting  size  excluded. 
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Table  A-8  lists  the  measured  physical  properties  of  different  materials  by  material  state, 

i.e.,  virgin,  backface  or  fired.  The  virgin  state  refers  to  the  as-received  material,  the  fired 
state  refers  to  the  fired  surface,  and  the  backface  state  refers  to  the  material  a small  distance 
below  the  fired  surface.  Within  a state,  the  materials  are  divided  into  three  groups  by  material 
class.  These  three  groups  are: 

1.  Bulk  graphite 

2.  Carbon/carbon  (c/c) 

3.  Pyrolytic  graphite  (PG) 

In  order  to  make  a comparison  between  different  material  classes  at  various  material  states, 
average  material  properties  are  plotted  in  Figures  A- 35  through  A-37. 

In  preparing  these  figures,  only  those  materials  for  which  data  was  available  in  all  three 
states  were  utilized. 

Density  values  of  the  virgin  state  are  greater  than  those  at  the  fired  state  as  seen  in 
Figure  A-35.  As  expected,  the  porosity  values  show  just  the  reverse  behavior  as  seen  in  Figure 
A- 36.  As  observed  in  Figure  A-2,  while  the  total  porosity  of  c/c  materials  is  greater  than  the 
others,  its  open  porosity  is  smaller  than  that  for  bulk  graphite. 

Figure  A-37  presents  the  average  specific  open  pore  volume  by  material  class  for  various 
materials  conditions.  The  shapes  of  the  curves  are  similar  to  tnose  of  the  open  porosity  curves. 

Data  for  pyrolytic  graphite  is  insufficient  to  draw  any  convincing  conclusions  regarding  the  pro- 
perties behavior. 

Figures  A- 38  through  A-40  show  typical  pore  size  spectrums  for  selected  materials.  There 
was  no  consistent  pore  size  occurrence  within  a given  material  class  for  a given  material  state. 

For  example,  the  pore  size  distribution  for  G-90  shows  that  about  22  percent  of  the  pores  fall  in 
the  diameter  range  from  0.069  to  0.035  microns.  The  remaining  porosity  is  scattered  around  it  from 
3.43  to  0.013  microns.  On  the  other  hand,  for  ATJ-S  about  15  percent  of  the  pores  range  from  4.66 
to  3.42  microns  in  diameter,  about  13  percent  range  from  0.023  to  0.018  micron  in  diameter  and 
others  are  scattered  around  these  intervals  from  17.0  to  0.013  microns  in  diameter. 

Except  for  a few  exceptions,  it  cu.i  be  said  that  the  majority  of  the  pores  lie  in  the  range 
of  about  17  to  0.013  microns  in  diameter  regardless  of  material  class  or  state.  A shift,  however, 
towards  smaller  pore  diameters  with  progression  from  virgin  to  backface/ fired  conditions  was  observed. 
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TABLE  A-8.  POROSITY  MEASUREMENT  RESULTS 
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Average  density  tg/cc) 


Figure  A-35.  Average  density  of  specimen  materials  at 
different  material  states. 
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Virgin 
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Material  state 


Figure  A-37.  Average  specific  open  pore  volume  by  material 
class  at  different  material  states. 
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Figure  A-40.  Typical  pore  size  distribution,  15%  SiC/PG 


Limited  data  was  developed  for  vapor  deposited  carbons,  but  it  was  observed  that  the  pore 


structure  was  radically  different  compared  to  bulk  graphites  or  carbon/carbons. 

REFERENCE  FOR  APPENDIX  A 

A-l . Hughes,  M.  C.,  et  al.,  "Codeposited  PG/SiC  Nozzle  Liners  for  Advanced  ICBM  Systems  Vol.  I 
Deposition  Process  Development,"  AFRPL-TR-74-15,  April  1974. 
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APPENDIX  B 


GRAPHITIC  MATERIALS  THERMOPHYSICAL  PROPERTIES 


B.l  INTRODUCTION 

There  are  many  graphitic  type  materials  typically  used  in  the  construction  of  solid  propel- 
lant rocket  motor  nozzles.  It  is  important  that  the  thermal  properties  of  these  various  materials 
be  adequately  established  and  generally  accepted  for  meaningful  and  consistent  predictions  of  ther- 
mal response  of  rocket  motor  nozzles.  The  purpose  of  this  appendix  is  to  provide  a collection  or 
recommended  property  values  for  selected  types  of  yraphitic  materials.  The  properties  collected 
are  those  that  are  relevant  to  the  prediction  of  thermal  response,  i.e.,  specific  heat,  thermal  con- 
ductivity, density,  heat  of  formation  and  emissivity.  For  each  type  of  material  the  following  in- 
formation is  provided. 

• List  of  manufacturers 

• Manufacturers  suggested  property  values 

• Uncertainties  in  property  values 

• Sensitivity  of  thermal  response  to  uncertainties  in  property  values 

• Recommended  property  values 

This  appendix  is  organized  so  that  additional  types  of  materials  or  new  materials  within  an  already 
included  type  may  be  easily  added. 

The  recommended  property  values  result  from  a consideration  of  the  sensitivity  of  the  ther- 
mal response  to  an  assumed  uncertainty  in  material  properties.  This  sensitivity  analysis  was  per- 
formed for  the  nozzle  throat  insert  of  an  MX  nozzle  using  ablation  rate  as  the  principal  criteria. 
Since  flow  conditions  in  other  regions  of  the  nozzle  are  significantly  different,  the  conclusions  on 
throat  response  sensitivity  do  not  necessarily  apply  elsewhere. 

Material  property  uncertainties  may  arise  from: 

• Lack  of  property  data 
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• Errors  in  property  measurements 

• Process  variations  (billet  to  billet) 

• Within-billet  variation 

The  thermal  response  is  characterized  by  the  time  variation  of: 

• Surface  temperature 

• Surface  recession 

• Backwal 1 temperature 

For  the  purposes  of  the  sensitivity  analysis  the  environmental  conditions  were  assumed  to  be: 

• HTPB  propellant  (Pc  1 1000  psla,  Tc  * 6642°R) 

a M « 1 

• P»u  Cu  = 0.7  - 14  lbm/ftJsec 

e e m 

• CM/CH  = 0.6  - 0.8 

The  propellant  gas  properties  and  B‘  versus  T maps  were  calculated  by  the  GASKET  code  (Reference  B-l) 
using  the  most  appropriate  surface  kinetic  model  Included  In  the  GASKET  code.  The  thermal  response 
was  calculated  by  the  CMA  code  (Reference  B-2)  for  the  conditions  given  below. 

• Material  thickness  equal  to  2 inches  (unless  otherwise  noted) 

• Insulated  backwal  1 

• 40-second  exposure 

Only  the  principal  results  of  the  sensitivity  study  are  presented  in  each  Section  to  illustrate  the 
sensitivity  of  the  thermal  response  to  the  various  properties.  General  conclusions  regarding  the 
thermal  response  are: 

• The  thermal  response  is  Insensitive  to  surface  emissivlty  and  heat  of  formation  of  the 
particular  graphitic  material 

• The  surface  temperature  is  insensitive  to  values  for  the  specific  heat,  thermal  conduc- 
tivity and  density 

• For  surface  temperatures  less  than  approximately  6500°R* 

n 

This  requirement  arises  from  limitation  of  the  slope  of  the  B'  versus  T curve.  In  the  region  of 
large  slope,  near  the  sublimation  region,  the  conclusions  may  not  be  valid. 
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- The  mass  loss  rate  and  surface  recession  are  insensitive  to  specific  heat  and  thermal 
conductivity  variations 

- The  surface  recession  is  directly  proportional  to  density  variation 

• The  backwall  temperature  is  sensitive  to  variations  in  specific  heat  and  thermal  con- 
ductivity 

These  conclusions  are  valid  for  the  types  of  graphitic  materials  examined  and  for  typical  MX  nozzle 
throat  conditions.  They  are  not  dependent  upon  the  particular  surface  kinetic  model  or  propellant 
(for  XLDB,  PEG/FEFO  or  HTPB). 

Properties  to  be  presented  are  categorized  into  four  generic  classes,  namely: 

t Bulk  graphite 

• Carbon/carbon 

• Pyrolytic  graphite 

• Modified  pyrolytic  graphite 

B.2  BULK  GRAPHITES 

Bulk  graphites  are  in  widespread  use  in  nozzle  thermal  protection  systems  which  operate  at 
high  temperatures  and  pressures.  Presented  herein  are  a brief  description  of  the  manufacturing  pro- 
cess, a list  of  manufacturers  of  aerospace  grades  of  bulk  graphite,  representative  thermal  proper- 
ties, and  the  results  of  the  analysis  of  the  sensitivity  of  the  thermal  response  to  property  varia- 
tion. 

B.2.1  Manufacturing  Sumnary 

Manufactured  graphite,  in  aerospace  grades,  is  produced  in  a similar  fashion  by  all  manufac- 
turers. A filler,  usually  petroleum  coke,  and  a binder,  coal  tar  pitch,  are  mixed  and  then  formed 
to  a shape.  This  is  baked  to  form  amorphous  carbon  which  is  graphitized  at  temperatures  ranging 
from  40C0°F  to  5500°F.  Variations  in  each  step  have  a significant  effect  on  the  properties  of  the 
finished  product  and  contribute  to  the  differences  observed  among  the  various  conmercial  grades  pro- 
duced. In  addition  to  these  standard  processes,  certain  grades  of  graphite  are  further  densit.ed, 
either  by  hot-working  finished  billets  or  by  multipitch  impregnations  and  regraphitization  of  the 
finished  billets.  These  processes  produce  a material  that  is  quite  anisotropic  with  respect  to  ther- 
mal and  structural  properties.  The  highest  degree  of  anisotropy  results  from  extrusion  of  a coarse 
grain  material;  wherease,  an  isostatically  molded,  fine  grain  material  can  be  practically  isotropic. 
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8.2.2  Materials  Summary 


Manufacturers  and  materials  for  which  thermal  property  data  are  available  are  shown  in  Table 
b-l.  Detailed  values  for  thermal  conductivity  and  specific  heat  are  presented  in  Table  B-2  through 
Table  B-ll.  Figures  B-l  through  B-3  show  these  same  thermal  properties  in  graphical  form.  It 
should  be  noted  that,  in  many  cases,  material  properties  are  obtained  from  tests  conducted  by  in- 
dependent organizations,  such  as  SoRI,  ARC,  LMSC,  6E,  and  Aerotherm,  rather  than  the  material  manu- 
facturer. 

The  heat  of  formation  for  the  various  bulk  graphites  is  assumed  to  be  zero  (the  value  for 
elemental  carbon).  The  surface  emissivity  (total  hemispherical)  is  initially  a function  of  the 
surface  finish;  however,  once  the  surface  starts  ablating  all  the  graphites  have  essentially  the 
same  value,  approximately  0.90. 

For  all  the  materials  investigated  it  was  necessary  to  estimate  probable  thermal  property 
variations  within  one  grade  (billet  to  billet  variations)  and  probable  variations  of  properties 
v 'thin  one  billet  of  material.  There  is  general  agreement  that  these  variations  do  exist  and  that 
they  may  have  an  impact  on  design  considerations.  No  large  body  of  data  was  found  treating  this 
area;  however,  References  B-l 0 and  B-ll  offer  some  Insight  to  the  problem.  Reference  B-ll  is  a 
study  of  ATO-S  graphite  and  its  property  variations  and  is  used  as  a source  for  the  magnitude  of 
variations.  Variations  In  thermal  conductivity  of  ±14  percent  with  respect  to  the  nominal  values 
are  reported  for  billet  to  billet  variations;  variation?  of  specific  heat,  density,  and  emissivity 
were  found  to  be  negligible.  This  reference  also  considered  within-billet  variations  and  found 
that  they  do  exist,  but,  from  a thermal  p”operty  standpoint,  these  were  negligible. 


B . 2 . 3 Sensitivity  Analysis  and  Recommended  Values 

The  results  of  the  sensitivity  analysis  for  bulk  graphite  are  shown  in  Table  B-12  and  are 
consistent  with  the  general  conclusions  given  in  Section  B.l.  Nominal  property  values  for  this  sen- 
sitivity study  were  taken  as  those  of  ATJ-S.  The  variations  in  surface  emissivity  and  specific  heat 
wers  arbitrarily  selected  to  determine  their  impact  on  the  thermal  response.  The  variation  in  den- 
sity covered  the  range  of  values  for  all  the  bulk  graphites  about  a nominal  value  of  1.83  g/cc. 

The  thermal  conductivity  was  allowed  to  vary  by  ±14  percent  to  reflect  the  reported  billet  to  billet 
variations.  (The  error  for  measurement  of  thermal  conductivity  was  reported  as  ±7  percent  in  Refer- 
ence B-l 5. ) The  nominal  values  of  Tg^,  T^,  and  recession  for  the  with-grain  case  are  shown  in  Table 
B-12.  The  results  for  the  sensitivity,  i.e.,  the  percent  fluctuation  in  Tgy,  Tg,  and  recession, 
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TABLE  B-l.  BULK  GRAPHITE  MANUFACTURERS 


Manufacturer 

Product  (s) 
(Density,  G/CC) 

Carborundum  Co. 

G-90  (1.90),  Graphitite-G  (1.88) 

Great  Lakes  Carbon  Corp. 

H-205-85  (1.81) 

Poco  Graphite,  Inc. 

AXF-5Q  (1.81) 

Pure  Carbon 

P-03  (1.83) 

Speer  Carbon 

8882-E  (1.76) 

Union  Carbide 

ATJ  (1.73),  ATJ-S  (1.83), 
AGSR  (1.55),  CS  (1.72) 

TABLE  B-2.  THERMAL  PROPERTIES  ATJ  (Reference  B-12) 


Temperature 

°R 

Specific 

Btu/lbm- 

460 

0.283 

960 

.340 

1460 

.390 

1960 

.430 

2460 

.462 

3450 

.505 

4460 

.521 

5460 

.525 

6460 

.525 

Thermal  Conductivity 
Btu/ft-sec-cR  x 10" 3 

W/G 

A/G 

20.52 

15.6 

16.08 

11.76 

12.36 

9.37 

9.67 

7.96 

7.86 

6.29 

6.46 

5.09 

5.78 

4.74 

5.39 

4.56 

5.32 

4.44 

4460 

5460 

6460 


TA8LE  B-3,  THERMAL  PROPERTIES  ATJ-S  (Reference  B-ll) 


Temperature 

°R 

Specific  Heat 
Btu/lbm-°R 

Thermal  Co 
Btu/ft-sec 

nductivity 
-°R  x 10"3 

v;/g 

A/G 

460 

0.283 

25.80 

19.56 

960 

.340 

18.72 

14.88 

1460 

.390 

14.64 

11.55 

1960 

.430 

11.74 

9.16 

2460 

.462 

9.70 

7.43 

3460 

.505 

7.32 

5.00 

4460 

.521 

6.34 

4.66 

5460 

.526 

5.3b 

4.42 

6460 

.526 

5.44 

4.16 

Density  = 114.2  lbm/ft3  = 1.83  g/cc 

zu 


1 
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TABLE  B-5.  THERMAL  PROPERTIES  A6SR  (Reference  B-12) 


Temperature 

°R 

Specific  Heat 
Btu/lbm-°R 

1 

Thermal  Conductivity 
Btu/ft-sec-°R  x 10'3 

W/G 

A/G 

460 

| 

0.283 

24.48 

19.44 

960 

.340 

18.72 

15.00 

1460 

.390 

14.28 

11.34 

1960 

.430 

11.27 

8.90 

2460 

.462 

9.26 

7.36 

3460 

.505 

6.95 

5.51 

4460 

.521 

6.25 

4.86 

5460 

.525 

6.20 

4.86 

6460 

.525 

6.20 

4.86 

■ ■ i 

Density  = 96.77  lbm/ft3  = 1.55  g/cc 


TABLE  B-6.  THERMAL  PROPERTIES  CS  (Reference  B-12) 


Temperature 

°R 

Specific  Heat 
Btu/lbm-°R 

Thermal  Conductivity 
Btu/ft-sec-°R  x 10~3 

W/fr 

A/G 

460 

0.283 

25.68 

20.40 

960 

.340 

20.88 

16.32 

1460 

.390 

16.68 

12.96 

I960 

.430 

13.56 

10.66 

2460 

.462 

11.17 

8.80 

3460 

.505 

8.80 

6.94 

4460 

.621 

7.82 

6.13 

5460 

.525 

7.40 

5.78 

6460 

.525 

7.30 

5.78 

Density  = 107.6  lbm/ft3  = 1.72  g/cc 
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TABLE  B-7.  THERMAL  PROPERTIES  G-90  (Reference  B-12) 


Temperature 

°R 

Specific  Heat 
Btu/1bm-°P 

Thermal  Conductivity 
Btu/ft-sec-°R  x 10' 3 

W/G 

A/G 

460 

0.283 

30.60 

21.00 

960 

.340 

22.68 

16.92 

1460 

.390 

16.68 

13.44 

1960 

.430 

12.96 

10.76 

2460 

.462 

10.88 

9.14 

3460 

.505 

8.10 

7.16 

4460 

.521 

7.16 

6.36 

5460 

.525 

6.95 

6.13 

6460 

.525 

6.95 

6.13 

Density  * 118.5  Ibm/ft’  =1.9  g/cc 

8-10. 


THERMAL  PROPERTIES  P-03- ISOTROPIC  (Reference  B-13) 


Temperature 

°R 

Specific  Heat 
Btu/lbm-°R 

Thermal  Conductivity 
1 Btu/ft-sec-°R  x 10"3 

1110 

0.350 

27.768 

1460 

.390 

21.672 

1936 

.428 

14.580 

2460 

.462 

10.692 

3460 

.505 

7.908 

4460 

.521 

7.224 

5460 

.525 

6.948 

Density  * 114  lbm/ft3  = 1.83  g/cc 

B-1 1 . THERMAL  PROPERT! 


Temperature 

°R 

Specific  Heat 
Btu/lbm-°R 

460 

0.340 

1460 

.390 

1960 

.430 

2460 

.462 

3460 

.505 

4460 

.521 

6460 

.525 

Density  = 109.9  lbm/ft3  = 1.76  g/cc 
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SPEER  8882- E (Reference  B-14) 


Thermal  Conductivity 
Btu/ft-sec-°R  X 10“3 

W/G 

A/G 

53.16 

44.04 

36.96 

25.44 

21.00 

14.04 

12.96 

9.80 

10.50 

9.06 

10.50 

9.06 

10.50 

9.06 

Figure  B-l.  Thermal  conductivity  of  bulk  graphites  against  grain  direction 
(P-03  is  isotropic). 
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ivity,  Btu/ft-sec- 


Specific  heat  - 8tu/lbm 


k-,4 


Figure  B-3.  Specific  heat  of  bulk  graphites. 
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TABLE  B-12.  EFFECT  OF  PROPERTY  VARIATION  ON  THERMAL 
RESPONSE  FOR  BULK  GRAPHITES 


Property 

Variation* 

tbw 

(2320°R) 

Ts 

(5490°R) 

Recess ion+ 
(.47  in) 

£ 

+30% 

+4% 

+0.5% 

+2% 

p 

+10% 

— 

— 

+10% 

k 

+14% 

+7% 

+0.1% 

+2% 

CP 

+10% 

+6% 

+2% 

+1% 

Worst  combination 
of  Cp  and  k 

+16% 

+2% 

+3% 

* 

ATJ-S  selected  for  nominal  values;  with-grain  results  shown  in 
parenthesis  above. 

+ 

Evaluated  at  end  of  40  seconds. 
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were  the  same  for  both  the  with  grain  and  against  grain  values  of  k.  Also,  a severe  change  in  the 
surface  kinetic  model  (from  the  bulk  graphite  model  to  the  layer  PG  model)  had  very  little  effect 
on  the  sensitivity  to  property  variations.  (The  sensitivity  to  the  kinetic  model  itself  was  not  con- 
sidered in  this  investigation.) 

From  the  results  of  the  sensitivity  analysis  and  the  availability  and  scatter  of  materials 
properties,  the  following  property  values  are  recommended: 

• e =0.9 

• hQ  = 0.0  Btu/lbm 

• Cp  - see  Figure  B-3 

• k,p  = a.  Use  manufacturers  recommended  values 

b.  Use  values  for  G-90  or  ATJ-S  as  representative  values 

B.3  CARBON/CARBON  COMPOSITES 

This  section  contains  a brief  description  of  the  manufacturing  process,  a list  of  carbon/ 
carbon  manufacturers  and  thermal  properties  of  their  products,  and  the  results  of  a sensitivity 
study  of  the  thermal  response  to  possible  property  variations  for  carbon/carbon  composite  materials. 

B.3.1  Manufacturing  Summary 

The  class  of  materials  known  as  carbon/carbon  composites  is  a very  broad  category  and  there 
are  many  variations  in  material  processing.  Also,  the  precise  details  of  these  processes  are  pro- 
prietary; so,  only  a brief  and  very  general  description  can  be  given  for  the  processing  of  carbon 
composites.  (See  Reference  B-30  for  more  details  on  the  processing  of  carbon  composites.) 

A typical  2-D  carbon/carbon  composite  starts  out  as  a carbon  or  graphite  phenolic,  which  is 
densified  through  a combination  of  liquid  phase  Impregnation  and/or  gas  phase  (Carbon  Vapor  Deposi- 
tion, CVD)  densification.  This  material  Is  then  carbonized  and  graphitized  or  annealed  at  tempera- 
tures in  excess  of  5000°R.  It  is  not  unusual  for  the  carbon  composite  materials  to  go  through  sev- 
eral impregnation  or  CVD  cycles  with  graphitization  or  annealing  at  the  appropriate  time.  The  manu- 
facture of  3-D  composites  is  somewhat  different,  as  data  from  AVCO  shows.  A preform  is  first  con- 
structed consisting  of  graphite  fabrics  pierced  with  graphite  fibers.  This  preform  is  impregnated 
with  a phenolic  resin;  carbonized  at  ~2460°R  and  then  graphitized  at  ~5460°R.  The  great  diversity 
among  the  carbon  composites  results  from  all  the  possible  variations  of  process  cycles  and  condi- 
tions that  are  possible. 
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B.3.2  Materials  Survey 

Table  B-13  lists  the  companies  and  the  products  for  which  thermal  properties  were  made  avail- 
able. Values  for  the  density,  specific  heat,  and  thermal  conductivity  are  given  in  Tables  B-14 
through  B-19  for  these  materials.  This  same  information  is  presented  graphically  in  Figures  B-4 
through  B-6  where  the  differences  in  thermal  conductivity  are  apparent.  The  heat  of  formation  can 
be  assumed  to  be  zero  and  the  surface  emissivity  can  be  taken  as  0.90.  It  should  be  noted  that 
these  composite  materials  are  not  as  well  characterized  as  the  bulk  graphites  and  that  the  proper- 
ties are  strongly  dependent  on  the  particular  processing.  The  density  may  even  vary  significantly 
for  the  same  material  specification,  for  example,  Pyrocarb  901  is  given  here  with  p = 1.6  g/cm3 
whereas  a sample  of  Pyrocarb  901  tested  by  Aerotherm  had  p = 1.8  g/cm3. 

Billet  to  billet  and  within  billet  variations  of  properties  are  discussed  in  Reference  B-31 
and  B-32.  The  density  both  within  a billet  and  between  billets  was  found  to  vary  as  much  as  ±5  per- 
cent. Thermal  conductivity  variations  can  be  expected  to  be  ±7  percent  due  to  material  variations 
and  experimental  error  (Reference  B-31).  The  specific  heat  is  not  always  reported  and  in  some  cases 
must  be  approximated  from  values  for  other  carbon/carbon  materials. 

B.3.3  Sensitivity  Analysis  and  Recommended  Values 

The  significant  results  of  the  sensitivity  study  for  the  carbon/carbon  composites  are  given 
in  Table  B-20.  The  sensitivity  to  emissivity  and  heat  of  formation  are  not  shown  since  thermal  re- 
sponse was  found  to  be  insensitive  to  either.  The  property  values  for  Haveg  SP8040  were  used  as 
nominal  values.  Variations  in  Cp  and  k for  this  sensitivity  study  were  larger  than  those  used  in 
the  graphite  studies.  This  was  done  tor  two  reasons: 

t Property  values  vary  widely  for  the  various  materials  classed  as  carbon/carbon  composites 

t Qualitative  results  for  small  variations  were  evaluated  in  the  bulk  graphite  section 
(Section  B.2) 

The  bulk  graphite  surface  kinetic  model  in  GASKET  was  used  since  at  the  time  of  these  studies  there 
were  no  models  for  carbon/carbon  materials.  This  is  acceptable  since  the  results  of  the  bulk  graph- 
ite sensitivity  study  showed  that  the  sensitivity  to  property  variations  was  not  dependent  cn  the 
surface  kinetic  model.  The  results  shown  in  Table  B-20  support  the  general  conclusion  of  Section  B.l. 

For  detailed  calculations  the  manufacturers  recommended  values  should  be  used.  The  MOD  3 ma- 
terial is  a three-dimensional  weave  and  should  not  be  used  as  representative  of  two-dimensional  ma- 
terials. If  property  values  are  not  known,  the  SP8040  values  can  be  used  as  representative.  For 
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♦Estimated 

4* 

The  thermal  conductivity  in  the  45°  x-y  direction  was  also  reported 
in  Reference  B-30  and  is  different  from  the  above  values. 


% 
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TABLE  B-16.  THERMAL  PROPERTIES  C/C  700  (Reference  B-30) 


Temperature 

°R 

Specific  Heat 

R4>.  | / 1 km 

Thermal  Conductivity 
Btu/ft-sec  °R  x 10'3 

DlU/  1 DIN 

W/P 

A/P 

530 

0.18 

15.0 

5.62 

960 

0.31 

9.5 

3.70 

1960 

0.43 

5.9 

2.67 

2460 

0.46 

5.66 

2.08 

2960 

0.48 

5.55 

1.86 

3460 

0.49 

5.7s 

1.97 

3960 

0.51 

5.75 

2.92 

4460 

0.51 

6.50 

2.77 

5460 

0.53 

7.30 

3.54 

5960 

0.53 

7.85 

3.85 

Density  = 93.5  lbm/ft3  = 1.5  g/cc 
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TABLE  B-17.  THERMAL  PROPERTIES  SP-8040  (Reference  B-32) 


Temperature 

°R 

Specific  Heat 
Btu/lbm-°R 

Thermal  Com 
Btu/ft-sec- 

ductivity 
3R  x 10"’ 

W/P 

A/P 

530 

0.245 

5.30 

2.50 

800 

0.260 

5.75 

2.95 

1100 

0.335 

6.30 

3.20 

1500 

0.450 

7.60 

3.40 

2000 

0.570 

7.95 

3.20 

2500 

0.580 

7.85 

2.80 

3000 

0.580 

7.60 

2.60 

3500 

0.580 

7.40 

2.70 

4000 

0.580 

7.40 

2.90 

4500 

0.580 

7.70 

3.30 

5500 

0.580 

9.30 

4.90 

Density  = 87.0  lbm/fts  * 1.4  g/cc 
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TABLE  B-18.  THERMAL  PROPERTIES  PYROCARB  502-1  (Reference  B-32) 


Temperature 

°R 

Specific  Heat 
Btu/lbm-°R 

Thermal  Conductivity 
Btu/ft-sec-°R  x 10" 3 

W/P 

A/P 

530 

0.240 

3.10 

1.50 

800 

0.275 

3.70 

1.52 

1100 

0.340 

4.30 

1.55 

1500 

P.415 

6.10 

1.90 

2000 

0.540 

6.35 

2.20 

2500 

0.54C 

5.75 

2.00 

3000 

0.540 

4.90 

1.90 

3500 

0.540 

5.40 

2.40 

4000 

C.540 

6.30 

3.00 

4500 

0.540 

7.35 

4.15 

5500 

0.540 

9.75 

7.60 

Density  = 74.5  lbm/ft}  = 1.2  g/cc 

H 
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TABLE  B-19.  THERMAL  PROPERTIES  FM-5228  (Reference  B-32) 


Temperature 

°R 

Specific  Heat 
8tu/1bm-°R 

Thermal  Con 
Btu/ft-sec 

ductivity 
°R  x 10-3 

W/P 

A/P 

530 

0.215 

4.15 

3.0 

800 

0.240 

4.70 

3.5 

1100 

0.315 

5.35 

3.8 

1500 

0. '.55 

8.30 

4.0 

2000 

0.590 

8.80 

3.9 

2500 

0.600 

8.70 

3.6 

3000 

0.600 

8.20 

3.4 

3500 

0.600 

8.10 

3.9 

4000 

0.600 

8.15 

5.0 

4500 

0.600 

8.50 

6.0 

5500 

0.600 

9.60 

8.0 

Density  = 94  lbm/ft3  = 1.5  g/cc 
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Temperature  ~ °R  x 103 


Figure  B-4.  Thermal  conductivities  of  carbon  composites,  against  ply. 
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Figure  B- 


0 2.0  3.0  4.0  5.0  6.0 

Temperature,  °R  x 10’ 


. Thermal  conductivities  of  carbon  composites,  with  ply. 
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TABLE  B-20.  EFFECTS  OF  PROPERTY  VARIATION  ON  THERMAL 
RESPONSE  FOR  CARBON-CARBON  COMPOSITES 


Property 

Variation* 

(2370  R) 

Ts+ 

(5460°R) 

Recession 
(0.625  in) 

P 

+10% 

— 

— 

+10% 

k 

+50% 

+40% 

+1% 

+7% 

+25% 

+25% 

+1% 

+2% 

*Haveq,  SP-8040  used  for  nominal  values,  with  ply  results 
shown  in  parenthesis  above* 

'i'Bulk  Graphite  kinetic  model  used  In  GASKET. 

Effect  of  property  variations  on  results 
are  the  same  for  wlth-ply  and  agalnst-ply 
values.- 
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simplicity,  constant  values  of  k,  both  with  and  against  ply,  and  Cp  could  be  used.  Although  there 
are  larger  differences  at  low  temperatures  (T  < 2500°R)  these  can  usually  be  neglected  since  the 
carbon/carbon  materials  would  probably  be  used  in  high  temperature  applications  and  the  low  tempera- 
ture transient  response  would  not  greatly  effect  the  results.  However,  when  tailored  properly  pro- 
perty carbon/carbon  materials  are  developed  to  the  point  where  they  are  used  as  insulators  as  well 
as  flame  surfaces,  then  low  temperature  properties  become  significant. 

B.4  PYROLYTIC  GRAPHITE 

Pyrolytic  graphite  is  a polycrystalline,  highly  anisotropic  form  of  graphite  produced  by  the 
thermal  decomposition  of  a hydrocarbon  gas.  PG  plates  (c  plane)  have  been  established  as  state-of- 
the-art,  however,  a-b  plane  coatings  have  yet  to  be  successfully  developed  for  ICBM  size  rocket 
nozzle  throat  inserts.  A brief  description  of  the  manufacturing  process,  the  reported  thermal  pro- 
perties and  the  results  of  a sensitivity  analysis  are  presented  in  this  section. 

B.4. 1 Manufacturing  Summary 

Pyrolytic  graphite  is  formed  by  the  vapor  deposition  of  a hydrocarbon  gas,  typically  methane; 
however,  other  hydrocarbons  such  as  acetylene  or  propane  are  often  used  either  by  themselves  or 
mixed  with  methane.  The  deposition  takes  place  on  a preformed  graphitic  substate  in  an  induction 
furnace  operating  at  temperatures  between  3500°R  and  4500°R. 

B.4. 2 Material  Survey 

Table  B-21  lists  some  pyrolytic  graphite  manufacturers.  Thermal  conductivity  c.nd  specific 
heat  values  used  by  the  rocket  motor  companies  are  shown  in  Figures  B-7  through  B-9. 

The  density  of  pyrolytic  graphite  is  2.2  gr/cc.  The  heat  of  formation  is  assumed  to  be  zero 
(the  value  for  elemental  carbon).  The  surface  emissivity  (total  hemispherical)  is  initially  a 
function  of  the  surface  fininsh,  however,  once  the  surface  starts  ablating,  all  of  the  graphites 
have  essentially  the  same  value,  approximately  0.90. 

B.4. 3 Sensitivity  Analysis 

Table  B-22  shows  only  the  results  of  the  sensitivity  analysis  for  variations  in  thermal  con- 
ductivity. All  other  results  are  consistent  with  those  reported  in  Section  B.l.  The  large  varia- 
tions in  thermal  conductivity  were  arbitrarily  selected  to  show  the  influence  of  very  large  property 
variations. 
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TABLE  B-21.  PYROLYTIC  GRAPHITE  MANUFACTURERS 


MANUFACTURER 

PRODUCT(S) 

Atlantic  Research  Corporation 

PURE  PG 

General  Atomic  Company 

ISOTROPIC  CARBON 

General  Electric  Company 

PURE  PG 

Hltco  Inc. 

PURE  PG 

Materials  Technology  Corporation 

PURE  PG 

Pfizer  Corporation 

PURE  PG 

Raytheon  Corporation 

PURE  PG 

Rocket  Propulsion  Establishment 

PURE  PG 

Super  Temp.  Corporation 

PURE  PG 

Union  Carbide  Corporation 

PURE  PG 
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Figure  B-7.  Thermal  conductivity,  a-b  plane. 


Temperature,  °R  x 103 


Figure  8-8.  Theraal  conductivity,  c direction. 


TABLE  B-22.  THERMAL  RESPONSE  FOR  LAYER  PG  FOR 
VARIATIONS  IN  THERMAL  CONDUCTIVITY 


Variation* 

Tpu 

(860°R) 

V 

(6200°R) 

Recession 
(0.0173  in) 

-80S 

-38S 

+1S 

+7S 

-50S 

-m 

+1S 

+AS 

+50S 

+13% 

-IS 

-3S 

+500S 

+55% 

-3S 

-17S 

★ 

Layer  PG  kinetic  model 

used  in  GASKET. 

*Aerotherm  PG  property  values  taken  as  nominal. 
100  mil  PG  on  2-inch  ATJ  substrate. 

244 


B.5  MODIFIED  PYROLYTIC  GRAPHITES 


8.5.1  Codeposited  Silicon  Carbide/ Pyrolytic  Graphite 

Codeposited  S1C/PG  is  a relatively  new  material  for  rocket  nozzle  inserts.  The  primary  manu- 
facturer is  Atlantic  Research  Corporation.  A brief  description  of  the  manufacturing  process,  the 
reported  thermal  properties,  and  the  results  of  the  sensitivity  analysis  for  SiC/PG  are  presented 
in  this  section. 

B.5. 1.1  Manufacturing  Summary 

Silicon  carbide  codeposited  pyrolytic  graphite  (Si f./PG)  is  formed  by  the  vapor  deposition  of 
a silicon  and  carbon  carrying  gas  (e.g.,  2 percent  methane,  25  percent  methyl  tri chi orosi lane  in 
nitrogen)  on  a preformed  graphite  substrate.  The  deposition  process  takes  place  in  an  induction 
furnace  operating  at  temperatures  between  3500°R  and  4500°R.  After  being  coated  the  material  is 
annealed  in  a nitrogen  environment  at  temperatures  of  42Q0°R  to  4700°R. 

B.5. 1.2  Material  Survey 

The  only  information  collected  on  the  thermal  properties  of  SIC/PG  is  that  provided  by 
Atlantic  Research  Corporation  and  Southern  Research  Institute  (References  B-50  to  B-52).  The  mate- 
rial tested  was  cut  from  one  billet  of  SIC/PG  which  was  nominally  20  percent  SiC  (by  weight)  with 
density  of  2.29  g/cc.  The  actual  SiC  content  varied  from  16  to  24  percent.  Reference  B-52  indi- 
cates that  the  a-b  plane  thermal  coroucti vity  is  only  slightly  affected  by  minor  variations  in  SiC 
content;  however,  the  c di>  ction  values  may  be  altered  by  as  much  as  50  percent  for  a 5 percent 
change  in  SiC  content.  In  addition  to  these  variations  the  measured  values  of  c direction  thermal 
conductivity  from  two  SoRI  programs  and  the  values  calculated  from  TRW  thermal  stress  tests  are 
greatly  different  (Figure  B-10).  Figure  B-ll  shows  the  values  of  thermal  conductivity  and  Figure 
B-12  shows  the  specific  heat.*  All  of  these  are  for  a nominal  20  percent  SiC  content. 

8. 5. 1.3  Sensitivity  Analysis  and  Recommended  Values 

The  recommended  values  for  thermal  conductivity  (Reference  B-50)  and  specific  heat  are  given 
in  Table  B-23. 

The  heat  of  formation  of  SiC/PG  can  be  calculated  from: 

Ho  = -7-86  PSiC 

* ”” 

Values  supplied  b.,  Atlantic  Research  Corporation. 
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TEMPERATURE  f'F) 


Figure  B-10.  Thermal  conductivity  of  PG/20X  SiC  in  the  "c1 
direction  (Reference  B-50). 

a.  Old  SoRI  data  generated  under  AIR  LAUNCH 

b.  New  SoRI  data  (Reference  B-50) 

c.  Data  backed  out  TRW  thermal  stress  tests 
and  recommended  for  nozzle  design 


Thermal  conduct!  vity,  Btu/ft-sec-°R 


Temperature,  °R  x 103 

Figure  B-ll.  Thermal  conductivity  of  20%  SiC/PG. 
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TABLE  B-23. 


THERMAL  PROPERTIES  SiC/PG  ATLANTIC 
RESEARCH  - 20%  SiC  (References  B-50 
TO  B-52 


Temperature 

Specific  Heat 
Btu/lbm-°R 

Thermal  Conductivity 
Btu/ft-sec°R  x 10"3 

a-b  plane 

c-direction 

525 

0.17 

33.6 

— 

960 

0.31 

28.8 

4.5 

1460 

0.37 

24.0 

3.46 

2460 

0.43 

15.6 

1,80 

2960 

0.45 

13.2 

1.23 

3460 

U.467 

11.3 

1.26 

4460 

0.475 

9.0 

1.45 

4960 

0.480 

9.0 

1.58 

5460 

0.480 

9.0 

— 

density  = 

142.96  1 bit)/ ft 3 = 

2.29  g/cc 

where  Is  the  weight  percent  of  SIC  and  HQ  Is  given  in  Btu/lbm.  It  should  be  noted  that  the 
thermal  response  is  very  Insensitive  to  the  heat  of  formation. 


TABLE  B-24.  EFFECTS  OF  PROPERTY  VARIATION  ON 
THERMAL  RESPONSE  FOR  SiC/PG  (20 
PERCENT  BY  WEIGHT  SiC) 


Property 

Variation 

tbw 

(2560°R) 

V 

(5400°R) 

Recession 
(0.33  in)+ 

P 

+10? 

— 

— 

+10% 

k 

+50% 

+50% 

+2% 

+2% 

c 

+25% 

+10% 

+1% 

+1% 

P 

•k 

Bulk  graphite  kinetic  model  used  in  GASKET. 
+a-b  plane  results  shown  in  parenthesis. 

Effect  of  property  variations  on  results  same 
for  a-b  plane  and  c direction  values. 
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